Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 625,777 A→G D143G (GAC→GGC)  cstA → pyruvate/proton symporter CstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273625,7770AG100.0% 116.3 / NA 33D143G (GAC→GGC) cstApyruvate/proton symporter CstA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/21);  total (12/21)

AGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTGCAGTGCTGGCGATGATCGTGGTGAAAGCCCT  >  NZ_CP009273/625643‑625906
                                                                                                                                      |                                                                                                                                 
aGCAGGTCCGCTGGTGGGTCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTCCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGt                                                                                                                               <  2:336476/139‑1 (MQ=255)
  cAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGtcg                                                                                                                             <  1:539698/139‑1 (MQ=255)
     gTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGtcgctc                                                                                                                          >  2:615115/1‑139 (MQ=255)
     gTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGgcgctc                                                                                                                          >  2:567203/1‑139 (MQ=255)
      tCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGtcgctcg                                                                                                                         >  2:490091/1‑139 (MQ=255)
                gggCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAg                                                                                                                <  1:405258/138‑1 (MQ=255)
                     ggTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTc                                                                                                          <  1:126007/139‑1 (MQ=255)
                            gcggcgCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAaag                                                                                                      >  2:125001/1‑136 (MQ=255)
                                      tGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGcc                                                                                         <  2:30718/139‑1 (MQ=255)
                                                               gctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCtg                                                                   <  2:276477/136‑1 (MQ=255)
                                                                ctgctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCtggtgg                                                               >  1:415063/1‑139 (MQ=255)
                                                                   ctCGCTGGGGTTGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                                                              <  1:346837/137‑1 (MQ=255)
                                                                   ctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                                                              <  1:490091/137‑1 (MQ=255)
                                                                   ctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                                                              <  2:393427/137‑1 (MQ=255)
                                                                   ctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                                                              <  1:561753/137‑1 (MQ=255)
                                                                   ctCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGc                                                              <  1:509797/137‑1 (MQ=255)
                                                                          gggTGGTTCTCGCCGGTGCGGTGCAGGATTTGATGTTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCATCCGGCGGGGTGATTGCGCTGGTGGCCTGCTTTa                                                      <  2:271404/138‑1 (MQ=255)
                                                                          gggTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTa                                                      <  1:541232/138‑1 (MQ=255)
                                                                                ttCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCAt                                               <  1:495321/139‑1 (MQ=255)
                                                                                     gCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGt                                            <  1:330203/137‑1 (MQ=255)
                                                                                          tgcggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATc                                     >  2:506339/1‑139 (MQ=255)
                                                                                          tgcggtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATc                                     >  2:413645/1‑139 (MQ=255)
                                                                                              gtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTg                                 <  1:174170/139‑1 (MQ=255)
                                                                                              gtgcAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTg                                 <  2:396045/139‑1 (MQ=255)
                                                                                                         cATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAg                                                                                                                >  1:126196/1‑49 (MQ=255)
                                                                                                         cATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAg                                                                                                                <  2:126196/49‑1 (MQ=255)
                                                                                                         cATGGTGCTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTGCAGTGCTGGc                       >  1:78767/1‑138 (MQ=255)
                                                                                                                cTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGAt                                                  <  1:578876/104‑1 (MQ=255)
                                                                                                                cTGTTTGTTTCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGAt                                                  >  2:578876/1‑104 (MQ=255)
                                                                                                                        ttCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGcc                                                                                         <  2:366401/57‑1 (MQ=255)
                                                                                                                        ttCTACGCGCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGcc                                                                                         >  1:366401/1‑57 (MQ=255)
                                                                                                                             cgcgCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTGCAGTGCTGGCGATGATCGTGGTGAAAGCCCt  >  2:214264/1‑139 (MQ=255)
                                                                                                                             cgcgCCGTGGCGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTGCAGTGCTGGCGATGATCGTGGTGAAAGCCCt  <  1:506339/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                 
AGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATGATCTGGCTGCTCGCTGGGGTGGTTCTCGCCGGTGCGGTGCAGGATTTCATGGTGCTGTTTGTTTCTACGCGCCGTGACGGTCGCTCGCTGGGTGAGCTGGTCAAAGAAGAGATGGGGCCAACCGCCGGGGTGATTGCGCTGGTGGCCTGCTTTATGATCATGGTCATTATCCTTGCAGTGCTGGCGATGATCGTGGTGAAAGCCCT  >  NZ_CP009273/625643‑625906

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