Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,951,452	A→G	T20A (ACC→GCC)	yecE →	DUF72 domain‑containing protein YecE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,951,452	0	A	G	100.0%	117.6 / NA	35	T20A (ACC→GCC)	yecE	DUF72 domain‑containing protein YecE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/3);  total (32/3)

AGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGTT  >  NC_000913/1951316‑1951585                                                                                                                                         |                                                                                                                                      agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCa                                                                                                                                    <  1:959814/140‑1 (MQ=255)    agcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACGTCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGcc                                                                                                                                 <  2:411798/140‑1 (MQ=255)                   tCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGccc                                                                                                                  >  2:259745/1‑140 (MQ=255)                   tCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGccc                                                                                                                  >  1:956329/1‑140 (MQ=255)                   tCAACCATATCCCCCCGCGCACCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTCCCGCAATGGTCGCATCCTAAATGGGGGCGGTTGGGGAGCGCCAGCCTTGAAGAGTATGccc                                                                                                                  <  2:985040/140‑1 (MQ=255)                      aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGcc                                                                                                               >  1:713170/1‑140 (MQ=255)                           atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAaga                                                                                                                          <  2:548043/124‑1 (MQ=255)                           atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAaga                                                                                                                          >  1:548043/1‑124 (MQ=255)                           atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACttt                                                                                                          >  1:528110/1‑140 (MQ=255)                           atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACttt                                                                                                          >  1:57303/1‑140 (MQ=255)                           atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACttt                                                                                                          >  1:638192/1‑140 (MQ=255)                           atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACttt                                                                                                          >  1:84909/1‑140 (MQ=255)                                cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTg                                                                                                     >  1:60859/1‑140 (MQ=255)                                cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTg                                                                                                     >  1:447726/1‑140 (MQ=255)                                 ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGc                                                                                                    >  2:400290/1‑140 (MQ=255)                                  cgcgCACCGCCGGTGTGCGTGGGGGGGAAGAGCCCCCCAACGTTTTATGTTTTACATCGGTCTCCCGCAATGGTCGCATCAAAAAGGGGTGGGGTGGGGGACCGCCAGCCTTAAAGAGTATGCCCGCCACTTTAACTGCg                                                                                                   <  2:60859/140‑1 (MQ=255)                                  cgcgCACCGCCCGTGGGCGTGGGGGGAAAGAGATCCCCAACGTTTTATGATTTCCAGCGGTCTCCCGAAATGGCCGCACCCAAAGGGGGGGGGGTGGGGGACCGCCAGCTTGGAAGAGTAGGCCCGCCACTTTAACTGCg                                                                                                   <  2:447726/140‑1 (MQ=255)                                    cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTg                                                                                                 >  1:302880/1‑140 (MQ=255)                                               tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAacac                                                                                      >  1:638746/1‑140 (MQ=255)                                               tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAacac                                                                                      >  1:111281/1‑140 (MQ=255)                                               tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAacac                                                                                      >  2:298046/1‑140 (MQ=255)                                                          tGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACg                                                                           >  1:569593/1‑140 (MQ=255)                                                          tGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACg                                                                           >  2:504261/1‑140 (MQ=255)                                                          tGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACg                                                                           >  2:656395/1‑140 (MQ=255)                                                            gAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGcc                                                                         >  1:238413/1‑140 (MQ=255)                                                                       caccgctTTATTATTTACATCGGTCCCCCGCAGTGGCCGCATCAAAAAGGGGTGGGGTGGGGGATCGCCAGCTTTAAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAAcc                                                              <  2:638746/137‑1 (MQ=255)                                                                       cAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAAcc                                                              >  2:656183/1‑140 (MQ=255)                                                                           gcgTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAg                                                          >  2:883305/1‑140 (MQ=255)                                                                            cgTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGTGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAgg                                                         >  1:249664/1‑140 (MQ=255)                                                                                                    gCAAGGGTCGCATCAAAAGGGGGTGGGGTGGGGGCCCGCCAGCCTGGAAGAGTATCCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGa                                 <  2:956329/140‑1 (MQ=255)                                                                                                                tCCTAAATGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACt                     >  2:148139/1‑140 (MQ=255)                                                                                                                    aaaTGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCg                 >  2:115506/1‑140 (MQ=255)                                                                                                                    aaaTGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCg                 >  2:359496/1‑140 (MQ=255)                                                                                                                    aaaTGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCg                 >  2:183616/1‑140 (MQ=255)                                                                                                                    aaaTGGGTGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCg                 >  2:1034452/1‑140 (MQ=255)                                                                                                                           tGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGt          >  2:529694/1‑140 (MQ=255)                                                                                                                           tGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGt          >  1:11817/1‑140 (MQ=255)                                                                                                                           tGCGGTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGt          >  1:1150613/1‑140 (MQ=255)                                                                                                                               gTTGGGGATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTa      >  1:657179/1‑140 (MQ=255)                                                                                                                                   gggATCGCCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGtt  >  2:356525/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      AGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGTT  >  NC_000913/1951316‑1951585

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