Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,969,962	A→G	S341S (AGT→AGC)	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,962	0	A	G	100.0%	90.8 / NA	27	S341S (AGT→AGC)	tap	methyl‑accepting chemotaxis protein Tap
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/10);  total (17/10)

TCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGAT  >  NC_000913/1969827‑1970097                                                                                                                                        |                                                                                                                                        tcCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGa                                                                                                                                     >  1:666531/1‑140 (MQ=255)   cGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGAcc                                                                                                                                   >  2:520874/1‑140 (MQ=255)       gcgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGaa                                                                                                                               <  1:351419/140‑1 (MQ=255)       gcgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGaa                                                                                                                               <  1:758516/140‑1 (MQ=255)        cgAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAac                                                                                                                              >  1:347845/1‑140 (MQ=255)            ccGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACAccg                                                                                                                          <  1:559469/140‑1 (MQ=255)                 tCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACAccgcccgc                                                                                                                     <  2:154424/140‑1 (MQ=255)                       gCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTgcgc                                                                                                               >  1:28706/1‑140 (MQ=255)                       gCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTgc                                                                                                                 >  2:661120/1‑138 (MQ=255)                       gCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTgc                                                                                                                 <  1:661120/138‑1 (MQ=255)                           cATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGtt                                                                                                           >  2:35245/1‑140 (MQ=255)                                              ggTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCa                                                                                        <  2:432865/140‑1 (MQ=255)                                                  tGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTc                                                                                    >  1:211026/1‑140 (MQ=255)                                                    gAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCg                                                                                  >  2:69654/1‑140 (MQ=255)                                                     aaaGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTACATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGa                                                                                 <  2:525124/140‑1 (MQ=255)                                                                                aTGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTg                                                      <  2:69220/140‑1 (MQ=255)                                                                                     gCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTa                                                 >  2:125551/1‑140 (MQ=255)                                                                                     gCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGACCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTa                                                 >  2:730301/1‑140 (MQ=255)                                                                                       cAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTAcc                                               >  1:102825/1‑140 (MQ=255)                                                                                                  gAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACCCCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTaa                                    <  2:347845/140‑1 (MQ=255)                                                                                                         tggCGATCTCCTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTcc                             <  1:125551/140‑1 (MQ=255)                                                                                                                  ccTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTggcg                    >  2:261018/1‑140 (MQ=255)                                                                                                                  ccTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTggcg                    >  2:477901/1‑140 (MQ=255)                                                                                                                  ccTGCATGGTGTGAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTggcg                    <  1:69654/140‑1 (MQ=255)                                                                                                                           tgtgAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTtgt           >  2:642865/1‑140 (MQ=255)                                                                                                                           tgtgAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTtgt           >  1:429688/1‑140 (MQ=255)                                                                                                                           tgtgAGTCATGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTtgt           >  1:162948/1‑140 (MQ=255)                                                                                                                                    tGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGAt  <  1:595858/140‑1 (MQ=255)                                                                                                                                    tGGTGCTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGAt  >  2:595858/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        TCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGCATGGTGTGAGTCATGGTACTGACCTGAACACCGCCCGCCTGCGCCGTTGTCGCGGCATTTTTTGCCAGTTCCGACGCCTGTCGTGCGTTATCGGCGTTTTGCCCTACCGTGGCGGTTAATTGCTCCATACTGGCGGCCGTTTGTGCCAGCGAT  >  NC_000913/1969827‑1970097

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