Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,040,801 T→C L107L (TTA→CTA)  msrQ → periplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,8010TC100.0% 69.0 / NA 19L107L (TTA→CTA) msrQperiplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/7);  total (12/7)

TTTAACAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGTTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGG  >  NC_000913/2040755‑2040844
                                              |                                           
tttAACAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGctat                                          <  1:2319340/50‑1 (MQ=255)
  tAACAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGctatt                                         <  1:1874128/49‑1 (MQ=255)
   aaCAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGctatta                                        <  1:1091938/49‑1 (MQ=255)
    aCAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAgg                                      <  1:1272816/50‑1 (MQ=255)
    aCAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAgg                                      <  1:698976/50‑1 (MQ=255)
    aCAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAg                                       <  1:1977449/49‑1 (MQ=255)
     cAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAgg                                      >  1:72329/1‑49 (MQ=255)
       aGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGt                                     >  1:1361166/1‑48 (MQ=255)
          tACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAg                                  >  1:341580/1‑48 (MQ=255)
            cGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGt                              <  1:1240246/50‑1 (MQ=255)
             gCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGt                              >  1:655037/1‑49 (MQ=255)
               aTTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTaa                           >  1:1708648/1‑50 (MQ=255)
                   cTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTAc                       <  1:1979571/50‑1 (MQ=255)
                   cTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTAc                       <  1:620929/50‑1 (MQ=255)
                        gCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGc                  >  1:679079/1‑50 (MQ=255)
                          tGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCt                >  1:562209/1‑50 (MQ=255)
                               gTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTAtt            >  1:558112/1‑49 (MQ=255)
                                  aataatCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAAc        <  1:888316/50‑1 (MQ=255)
                                     aatCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGt      >  1:1865987/1‑49 (MQ=255)
                                     aatCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGt      >  1:709249/1‑49 (MQ=255)
                                       tCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAg   >  1:1073872/1‑50 (MQ=255)
                                        cTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAgg  >  1:605331/1‑50 (MQ=255)
                                              |                                           
TTTAACAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGTTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGG  >  NC_000913/2040755‑2040844

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