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%	83.0 / NA	24	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 (9/15);  total (9/15)

TTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGC  >  NC_000913/2041041‑2041517                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                ttAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGctc                                                                                                                                                                                                                                     <  1:219099/250‑1 (MQ=255)        gaCGGTATCAGAAGTTGCGTTCCCTTCTTAACCGGTTACGCAATCAGGTTCACAACATATGATCTGTATAACTTCATATTGCCTACAAAGTACGCGTTATCCCTGCCGGATGTGTCGTGAGCGACCTATCCGGCATCCTCAGATCGCATCATGTTATGTTTCGTTTGACATCAACCTCAAATCACCAAAATAATTTATTTTATTGCAATAAGATAGTTAAGTAGAACTGAAACATCCATGCTCTCGCTTc                                                                                                                                                                                                                              <  2:104595/249‑1 (MQ=255)                          gTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCAtttt                                                                                                                                                                                                           <  1:3876/251‑1 (MQ=255)                             cccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCt                                                                                                                                                                                                        >  1:83743/1‑251 (MQ=255)                             cccTGTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGGTGCATTTTGc                                                                                                                                                                                                         <  1:206732/250‑1 (MQ=255)                              ccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAAGTATCGGTAGAACTTCATATTGCCGAAAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCta                                                                                                                                                                                                       <  2:205959/251‑1 (MQ=255)                                   ttAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGt                                                                                                                                                                                                  >  1:47544/1‑251 (MQ=255)                                                        ttCACAACAAATTCTCTGTATAACTTCATATTGCCGACAAAGTACGCTTTATCCCTGCCGGATGTGGCGCTTGCGCCCTATCCGGCATACACAGCTCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATCAGCTAGATAAGTAGCACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACACTATCACATTCCAcatc                                                                                                                                                                             >  2:221416/1‑251 (MQ=255)                                                             aacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCAtata                                                                                                                                                                        <  2:187497/251‑1 (MQ=255)                                                                                                    cgcgTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCttt                                                                                                                                 >  2:6084/1‑251 (MQ=255)                                                                                                         tATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGGTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCTATTCGTCTTTACaa                                                                                                                             <  2:108424/250‑1 (MQ=255)                                                                                                                 ccGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGCTCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAGAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGt                                                                                                                    <  2:162176/251‑1 (MQ=255)                                                                                                                              gAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTc                                                                                                       >  1:146965/1‑251 (MQ=255)                                                                                                                                                         gCATAAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCg                                                                            <  2:124588/251‑1 (MQ=255)                                                                                                                                                          cATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGc                                                                           <  2:104746/251‑1 (MQ=255)                                                                                                                                                                gttatggttccttcttcatcaACCTCCAATAACCAAAATACTTTATTTGATTTCGCTAAGCTGGATCAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGCTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGACACTGTTTCGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGtc                                                                     >  2:110038/16‑251 (MQ=255)                                                                                                                                                                                            aTCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCGCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACaga                                         <  1:96355/251‑1 (MQ=255)                                                                                                                                                                                                            attttatttCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAgct                         >  2:37535/1‑251 (MQ=255)                                                                                                                                                                                                            attttatttCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAgc                          <  2:177986/250‑1 (MQ=255)                                                                                                                                                                                                               ttatttCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAgctgc                       <  2:73397/250‑1 (MQ=255)                                                                                                                                                                                                                 atttCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTa                     >  2:90797/1‑250 (MQ=255)                                                                                                                                                                                                                           gatagataAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGttttt          <  1:6084/251‑1 (MQ=255)                                                                                                                                                                                                                                ataAGTAGAACTGAAAAAGACATGCTCTCGTTTCCTAAGCGATACTGAATTTTGCTATATGTTACAATCTAACTTTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGACGAACAAAAAGCTGATAATGGTGTTTTtgatgg     >  1:140144/1‑250 (MQ=255)                                                                                                                                                                                                                                   aGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGc  <  1:140279/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                TTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGC  >  NC_000913/2041041‑2041517

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