Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 853,311 G→A E601K (GAA→AAA)  SAUSA300_RS04120 → ribonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793853,3110GA100.0% 99.2 / NA 37E601K (GAA→AAA) SAUSA300_RS04120ribonuclease R
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/2);  total (35/2)

TAATCCGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTGAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAATATATGATTCAACATATTGGTGAT  >  NC_007793/853231‑853406
                                                                                |                                                                                               
tAATCCGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGa                                                                              >  1:4909634/1‑100 (MQ=255)
tAATCCGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGa                                                                              >  1:6295550/1‑100 (MQ=255)
 aaTCCGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGaa                                                                             >  1:875263/1‑100 (MQ=255)
   gccGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAAcg                                                                           >  1:617311/2‑100 (MQ=255)
    ccGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAAcgt                                                                          >  1:3567157/1‑100 (MQ=255)
    ccGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAAcgt                                                                          >  1:2970366/1‑100 (MQ=255)
    ccGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAAcgt                                                                          >  1:3366384/1‑100 (MQ=255)
       tAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTGAACGTGAAcgtcgt                                                                       >  1:4980179/1‑100 (MQ=255)
       aaaGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAAcgtcgt                                                                       >  1:1235683/2‑100 (MQ=255)
              ttAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTg                                                                <  1:3692496/100‑1 (MQ=255)
              ttAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTg                                                                >  1:3539857/1‑100 (MQ=255)
                    gagaAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAg                                                          >  1:1218264/1‑100 (MQ=255)
                           cAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATCGAGGCAGAACGTGa                                                   >  1:6111637/1‑100 (MQ=255)
                           cAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCATGAGATAGCTAAACATACTTCTAAACGTGAACGTCGTGCGATTGAGGCAGAACGTGa                                                   >  1:1711797/1‑100 (MQ=255)
                            tctGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATGGAGGCAGATCggaag                                                  >  1:2492062/3‑96 (MQ=255)
                            aaTGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGCGGCAGATCggaag                                                  >  1:1891220/1‑96 (MQ=255)
                                   aaCAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACtgatg                                           >  1:5736491/1‑100 (MQ=255)
                                   aaCAAAGAAGTGAAGCGTTGGGAAGACAAAATGCCTGAGTTAGCTACACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACtgatg                                           >  1:6539986/1‑100 (MQ=255)
                                            ttGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGaaaa                                  >  1:4154510/2‑100 (MQ=255)
                                            gTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACCTACTTCTAAACGTGAACGCCGTGCTATTGCGGCAGAACGTGATACTGATGAATTGaaaa                                  >  1:6298797/1‑100 (MQ=255)
                                              gAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGaaaaaa                                >  1:5251242/1‑100 (MQ=255)
                                              ctaGCGTTGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGaaaaaa                                >  1:1346800/3‑100 (MQ=255)
                                                aGCGTTGGGAAGACAAAGTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCCATTGAGGCAAAACGTGATACTGATGAATTGAAAAAAGc                              <  1:6507777/100‑1 (MQ=255)
                                                   cttGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGa                           >  1:3988227/2‑100 (MQ=255)
                                                    tggggAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGaa                          >  1:4161994/3‑100 (MQ=255)
                                                     tGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGCTGAATTGAAAAAAGCAGAat                         >  1:1862105/1‑100 (MQ=255)
                                                     tGGGAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGaaa                         >  1:400444/1‑99 (MQ=255)
                                                      gggAAGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAata                        >  1:896422/1‑100 (MQ=255)
                                                         gggACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGCGAACGTCGTGCTATTGAGGCAGAACGCGATACTGATGAATTGTAAAAAGCAGAatagat                     >  1:6612442/3‑100 (MQ=255)
                                                         acgACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAGAAAGCAGAatatat                     >  1:4931708/3‑100 (MQ=255)
                                                         aaGACAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAACAAAGCAGAatatat                     >  1:4962088/1‑100 (MQ=255)
                                                            aCAAATTGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGatcggaagag                  >  1:2937621/1‑91 (MQ=255)
                                                                  tGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAATATATGATTCAACa            >  1:4709411/1‑100 (MQ=255)
                                                                  tGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAATATATGATTCAACa            >  1:1193798/1‑100 (MQ=255)
                                                                  tGCCTGAGTTAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAATATATGATTCAACa            >  1:1589624/1‑100 (MQ=255)
                                                                          gtAGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAATATATGATTCAACATATTGGtg    >  1:5268322/2‑100 (MQ=255)
                                                                            aGCTAAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAATATATGATTCAACATATTGGtgat  >  1:5330294/1‑100 (MQ=255)
                                                                                |                                                                                               
TAATCCGTAAGTATTTAATTGAGAAATCAATGGATAACAAAGAAGTGAAGCGTTGGGAAGACAAATTGCCTGAGTTAGCTGAACATACTTCTAAACGTGAACGTCGTGCTATTGAGGCAGAACGTGATACTGATGAATTGAAAAAAGCAGAATATATGATTCAACATATTGGTGAT  >  NC_007793/853231‑853406

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑