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%	56.5 / NA	19	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 (11/8);  total (11/8)

AAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGAT  >  NC_000913/2041056‑2041512                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           aagacgTTGCGTTCCATTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAAATCGCATCAGGTTATGTTTCGTTTGTCATAAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGGTAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAgtt                                                                                                                                                                                                                >  2:114323/1‑251 (MQ=255)                       aCCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTAc                                                                                                                                                                                          <  1:195214/251‑1 (MQ=255)                                          tCACAACAAATTATCTGTATAACTTCATATTGACGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGACCTGAAAAAGCCATGCTCTCATTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTCCAcatcc                                                                                                                                                                       >  2:23117/1‑250 (MQ=255)                                            ccaatataTTATCTGGATAACGTCATATTTCCGACAAAGTACGCGTTTTCCCTGCCGGATGTGTCGCGAGCCCCCTATCCGGCCTACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAAGCACCAAAATAATTTATTTTATTTCAATAAGATATATAAGTAGAACTTAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATATCATTACACATCata                                                                                                                                                                     <  2:59413/244‑1 (MQ=255)                                                     tATCTGTATAACTTCATATTGCCGTCAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCTCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATCCATTaa                                                                                                                                                             <  2:217478/250‑1 (MQ=255)                                                          gTATAACTCCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTTTTGCATTTTGCTATATTTTACAATATAACATTACACATCATATACATTAActct                                                                                                                                                         <  1:182127/249‑1 (MQ=255)                                                                                          tATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACaa                                                                                                                        >  1:74043/1‑250 (MQ=255)                                                                                            tCCCTGCCGGTTGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCAttt                                                                                                                                                                                                       >  1:134768/1‑169 (MQ=255)                                                                                            tCCCGGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCAttt                                                                                                                                                                                                       <  2:134768/169‑1 (MQ=255)                                                                                                       tgtgGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCttt                                                                                                                            >  1:29118/1‑233 (MQ=255)                                                                                                       tgtgGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCttt                                                                                                                            <  2:29118/233‑1 (MQ=255)                                                                                                              cgAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGAGTCGGCTTTACAAACTGGCTGTTGCTTTAGgtgt                                                                                                   <  1:231732/251‑1 (MQ=255)                                                                                                                aGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCt                                                                                                 >  1:181780/1‑251 (MQ=255)                                                                                                                                      gATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCtt                                                                           >  1:231076/1‑251 (MQ=255)                                                                                                                                                                          cAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAAc                                       <  1:14388/251‑1 (MQ=255)                                                                                                                                                                                                   tttCAATACGATAGATAAGTAGAACTGAAAAAGCCATGCTTTCTTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGTAAACTGTTTTGGCGATAAATCTTTACAAACTGGCTGTTTCATTAGGGGTCTTTATTTTTAACGCTCCTGCCGTTTGGCCTGGTCATCACTAATACGGGAAACCCCTAACAGAGTTCGAACAAAAAgctgctcat              >  2:60491/1‑248 (MQ=255)                                                                                                                                                                                                         tAAGATCGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTGTTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGATGTCGTACAAAACGCTGCTAATgttgtt        >  2:72267/1‑251 (MQ=255)                                                                                                                                                                                                             atagataAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTACGAGTTGTTGCATTTTGCAATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGAGCTCCTGCCTTTTCGCCTGGTCATCACTCACACGGCAAACCCTTAACAGTGGTCGAACCAAAAGATGCTACTGGGGTTTTtg    >  2:54589/1‑251 (MQ=255)                                                                                                                                                                                                               agataAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTtgat  >  1:8294/1‑251 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           AAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGAT  >  NC_000913/2041056‑2041512

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