Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,351,156 T→C D71G (GAT→GGT)  sapA ← peptide ABC transporter substrate‑binding protein SapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,351,1560TC100.0% 76.9 / NA 22D71G (GAT→GGT) sapApeptide ABC transporter substrate‑binding protein SapA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/9);  total (13/9)

GGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTT  >  NZ_CP009273/1351070‑1351240
                                                                                      |                                                                                    
ggTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACACCca                                                                                   >  2:361234/1‑90 (MQ=255)
                 gcgcCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAAc                                                                  >  1:434046/1‑90 (MQ=255)
                            gAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAg                                                       >  2:503555/1‑90 (MQ=255)
                                   tCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCg                                                >  1:158394/1‑90 (MQ=255)
                                    cccAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGa                                               <  1:383421/90‑1 (MQ=255)
                                     ccAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGAc                                              >  2:224600/1‑90 (MQ=255)
                                                             cAGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTc                              <  1:346565/82‑1 (MQ=255)
                                                             cAGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTc                              >  2:346565/1‑82 (MQ=255)
                                                              aGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGAt                     <  1:224600/90‑1 (MQ=255)
                                                              aGGCGATAGGTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGAt                     <  2:243350/90‑1 (MQ=255)
                                                                      ggTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTaaa             <  2:423274/90‑1 (MQ=255)
                                                                      ggTATAGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTaaa             <  1:503555/90‑1 (MQ=255)
                                                                         ataGGGATCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGgt          >  2:142782/1‑90 (MQ=255)
                                                                               aTCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAAcc                                      >  2:60176/1‑56 (MQ=255)
                                                                               aTCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAAcc                                      <  1:60176/56‑1 (MQ=255)
                                                                                tCGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACt   >  2:150070/1‑90 (MQ=255)
                                                                                 cGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGAt                     <  2:241242/71‑1 (MQ=255)
                                                                                 cGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGAt                     >  2:330911/1‑71 (MQ=255)
                                                                                 cGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGAt                     <  1:330911/71‑1 (MQ=255)
                                                                                 cGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGAt                     >  1:241242/1‑71 (MQ=255)
                                                                                 cGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACtt  >  2:36800/1‑90 (MQ=255)
                                                                                 cGACACCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACtt  >  1:306390/1‑90 (MQ=255)
                                                                                      |                                                                                    
GGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTT  >  NZ_CP009273/1351070‑1351240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: