Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,041,278	G→A	intergenic (+160/‑97)	msrQ → / → zinT	membrane‑bound flavocytochrome MsrQ/metal‑binding protein ZinT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,041,278	0	G	A	100.0%	47.3 / NA	18	intergenic (+160/‑97)	msrQ/zinT	membrane‑bound flavocytochrome MsrQ/metal‑binding protein ZinT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/8);  total (10/8)

GCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGT  >  NC_000913/2041044‑2041503                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                  gCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGt                                                                                                                                                                                                                   >  1:337684/1‑251 (MQ=255)         gTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCAGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTa                                                                                                                                                                                                             >  1:323337/1‑249 (MQ=255)                       gTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCAtttt                                                                                                                                                                                             <  1:138292/251‑1 (MQ=255)                                                       nacaacaAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACAt                                                                                                                                                                       >  1:88257/2‑241 (MQ=255)                                                                                                    gTTATCCCTGCCGGATGTGGCTCGAGCGCCCTATCCGACCTACACAGATCGCATCAGGTTATGTTTCGTTTGTCATACACATCAAATCACCCCACTAATTTATTTTATTTCAATAAGATAGCTAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTAc                                                                                                                 >  1:347710/1‑250 (MQ=255)                                                                                                           cTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGc                                                                                                                              <  1:74699/230‑1 (MQ=255)                                                                                                           cTGCCGGATGTGGCGCGAGCGCCCTANCCGGCANNCANAGATCGCATCAGGTTATGTNTCNNNTNNNANCAACCNCANNTCACCAACANAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTNGTTTCCTAANNGTTGTTGCATTTTGCTATATGNTACNNTNNANCATNacn                                                                                                                                                                   >  2:74699/1‑192 (MQ=255)                                                                                                                                tccTATCCGGCATACACATATCGCATCAGGTTATGTTTCGTTTGTTATCAACCTCAAATCACCAAAATAAGTTAGTTTATTGGAATAATATAGAGAAGTAGTTCGGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACa                                                                                                                                                                                <  1:335784/158‑1 (MQ=255)                                                                                                                                nccTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAAAAACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCa                                                                                                                                                                                                                            >  1:72860/2‑115 (MQ=255)                                                                                                                                  cTATCCGGCATACACAGATCGCATCAGGTTATGTTTAGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTg                                                                                  >  1:301561/1‑251 (MQ=255)                                                                                                                                                     cGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTAGTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCttttn                                                               <  1:42102/251‑2 (MQ=255)                                                                                                                                                                      cgtttatcatgattcTCAAATCACCAGAAGAATTTAGTTTATTTCAATAAGAGATACAAGTAGATCTGAAAAATCCATTCTCTCGTTTCCTAAGAGTGGTTGCATGTTTCTATATGTTACAATATAACATTACACATCATATTCAt                                                                                                                                                       <  1:190019/132‑1 (MQ=255)                                                                                                                                                                            tcatcaACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATCAGTAGAACAGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATCACATTACACATCATATACATTAAc                                                                                                                                                   >  1:199767/1‑144 (MQ=255)                                                                                                                                                                            tcatcaACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACGCATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCtttt                                                                <  1:7493/227‑1 (MQ=255)                                                                                                                                                                                                                  cAATAAGATAGATAAGTAGAACTGAACAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATCTTAAAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGCGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGgt  >  1:201189/1‑251 (MQ=255)                                                                                                                                                                                                                              taAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCTCTCCTGCCtttt                                                                <  1:25316/177‑1 (MQ=255)                                                                                                                                                                                                                              taAGTAGAACTGAAAAAGCCATGCTCTCGATTCCTAAGAGGTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACNNTAGctcn                                                                                                                                                >  2:25316/1‑96 (MQ=255)                                                                                                                                                                                                                                 gTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCGGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTGTATTGTTAGCGCGCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACagagn                         <  1:64478/213‑2 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                  GCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGT  >  NC_000913/2041044‑2041503

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