Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,969,994	T→C	T331A (ACA→GCA)	tap ←	methyl‑accepting chemotaxis protein Tap

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,969,994	0	T	C	100.0%	100.8 / NA	30	T331A (ACA→GCA)	tap	methyl‑accepting chemotaxis protein Tap
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/13);  total (17/13)

GGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGA  >  NC_000913/1969860‑1970127                                                                                                                                       |                                                                                                                                      ggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCgg                                                                                                                                  >  2:485861/1‑140 (MQ=255) ggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCgg                                                                                                                                  <  1:87074/140‑1 (MQ=255) ggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCgg                                                                                                                                  <  1:485377/140‑1 (MQ=255) ggCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCgg                                                                                                                                  >  1:572719/1‑140 (MQ=255)           ttGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTTAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGc                                                                                                                        >  1:630868/1‑140 (MQ=255)                    aaaGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGa                                                                                                               <  1:819676/140‑1 (MQ=255)                                          tGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCg                                                                                         >  2:378450/1‑140 (MQ=255)                                            aTAATGTCGCCATTTTTTTGCGAGCTGGGGGCGATCTCCTGCGTGGTGTGAGTCAGGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGtt                                                                                       <  2:430085/140‑1 (MQ=255)                                            aTAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGtt                                                                                       <  2:572719/140‑1 (MQ=255)                                                    gCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTa                                                                               <  1:213315/140‑1 (MQ=255)                                                    gCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTa                                                                               >  2:52839/1‑140 (MQ=255)                                                          tttttGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTgg                                                                         <  1:290933/140‑1 (MQ=255)                                                              tGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGt                                                                     >  2:652602/1‑140 (MQ=255)                                                                    cTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTg                                                               <  1:652602/140‑1 (MQ=255)                                                                          gCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCAt                                                         >  2:363444/1‑140 (MQ=255)                                                                             aTCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACt                                                      >  1:384753/1‑140 (MQ=255)                                                                                      aTGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGt                                             <  1:378450/140‑1 (MQ=255)                                                                                          tgtgAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTtgt                                         >  2:79609/1‑140 (MQ=255)                                                                                          tgtgAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTtgt                                         >  1:420226/1‑140 (MQ=255)                                                                                          tgtgAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTtgt                                         <  1:178420/140‑1 (MQ=255)                                                                                                      tACTGACCTGAACACCGCCCACCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGGTCATTCCTCCATAATGGGGGCCGCTTTGGCCCGCGATgcc                             >  2:367362/1‑140 (MQ=255)                                                                                                         tGACCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATgccgcc                          >  1:588437/1‑140 (MQ=255)                                                                                                           aCCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCctg                        >  2:291724/1‑140 (MQ=255)                                                                                                           aCCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCctg                        <  1:456899/140‑1 (MQ=255)                                                                                                           aCCTGAACACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCctg                        >  2:73431/1‑140 (MQ=255)                                                                                                                aaCACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGtt                   >  2:437122/1‑140 (MQ=255)                                                                                                                aaCACCGCCCGCCTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGtt                   >  2:523597/1‑140 (MQ=255)                                                                                                                     cgcccgccTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTa              >  1:276656/1‑140 (MQ=255)                                                                                                                     cgcccgccTGCGCCGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTa              <  1:252740/140‑1 (MQ=255)                                                                                                                                 ccGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTgagaga  <  1:52839/140‑1 (MQ=255)                                                                                                                                 ccGTTGCCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTgagaga  >  1:26164/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                      GGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGATGCCGCCTGCTGTTCGGTACGACTTGAGAGA  >  NC_000913/1969860‑1970127

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