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%	111.6 / NA	31	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 (18/13);  total (18/13)

TAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATG  >  NC_000913/2041042‑2041516                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               tAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAACTCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCtctc                                                                                                                                                                                                                                  >  1:200434/1‑251 (MQ=255)       tACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTc                                                                                                                                                                                                                             <  2:257550/250‑1 (MQ=255)                                        ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAat                                                                                                                                                                                           >  1:62287/1‑251 (MQ=255)                                                            aacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCatat                                                                                                                                                                        >  2:88899/1‑250 (MQ=255)                                                              caaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACa                                                                                                                                                                     <  2:170569/251‑1 (MQ=255)                                                                    aTCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAAct                                                                                                                                                               >  2:118742/1‑251 (MQ=255)                                                                       tGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTg                                                                                                                                                            >  2:283702/1‑251 (MQ=255)                                                                            aaCTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACTTTACACATCATATACATTAACTCTggagg                                                                                                                                                        <  1:103816/250‑1 (MQ=255)                                                                             aCTTCATATTGCCGCCAAAGTACGCGTTATCCCTGCCGAATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTggagga                                                                                                                                                       >  2:59834/1‑250 (MQ=255)                                                                                        ccgtcacagttcgcgTTACCCCTGCCGGATGTGGCGCGAGCCCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTGTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGATCTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTTCATGTTGCTAGATGTTACAATAGAACATTACTCATCATATACATTAACTCTGGAGGATGCTGTTATGGc                                                                                                                                           <  2:222455/240‑1 (MQ=255)                                                                                          gACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGa                                                                                                                                         <  1:296267/251‑1 (MQ=255)                                                                                             aaaGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTc                                                                                                                                      <  1:206068/251‑1 (MQ=255)                                                                                              aaGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCg                                                                                                                                     >  1:176860/1‑251 (MQ=255)                                                                                               aGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCACAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGt                                                                                                                                    >  1:216613/1‑251 (MQ=255)                                                                                                                                  cccTATCTGGCATACACAGATCGCATCAGGTTATGTGTCGGTAGTCATCAACCACAAAGCACCAAAATAATTTATTTTATTTCAAGAAGATAGATAAGTAGAAATGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTGTGCTATCTGTTACAATATAACATGATACATCATATACAGTAACTCTTGAGGATATTGTTTTGGCGAGTGGTCAGTACAAACTGGCTTTTGCTTTATGTGTCttttt                                                                                                 <  2:232888/251‑3 (MQ=255)                                                                                                                                         cGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAgc                                                                                          <  1:119492/251‑1 (MQ=255)                                                                                                                                                         cATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACAtt                                                                                                                                                                                   >  2:248406/1‑146 (MQ=255)                                                                                                                                                         cATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACAtt                                                                                                                                                                                   <  1:248406/146‑1 (MQ=255)                                                                                                                                                          aTCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCga                                                                          >  2:244637/1‑249 (MQ=255)                                                                                                                                                             aGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTAGTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTATCTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAATCTGGCTGTTGCTTTAGGTGTCTTGATTGTTAGCGCTCCTGCCTGTGCGCATgg                                                                      <  2:135821/251‑1 (MQ=255)                                                                                                                                                                                  caACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAAcc                                                 >  1:248349/1‑251 (MQ=255)                                                                                                                                                                                             cACCAAAATACTTTATTTTATTTCAATAAGATAGCTAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAgg                                      >  2:134263/1‑251 (MQ=255)                                                                                                                                                                                                         ttattttattTCAATAAGATAGATAAGTCGAACTGAAAAAGACATGCTCTCGTTTCCTAAGAGTTGTTGAATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTAGTCTTTACAAACTAGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTACCAGAGGTCGAACaaaa                            >  2:314502/1‑249 (MQ=255)                                                                                                                                                                                                         ttattttattTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTTTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCTCTCGGCAAACCCTTAACAGAGTTCGAACaaaaa                           <  1:244637/250‑1 (MQ=255)                                                                                                                                                                                                         ttattttattTCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACaaaa                            >  2:44783/1‑249 (MQ=255)                                                                                                                                                                                                          tattttatttCAATAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAgc                         >  2:139949/1‑251 (MQ=255)                                                                                                                                                                                                                      aTAAGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGgtgt             >  1:2134/1‑251 (MQ=255)                                                                                                                                                                                                                        aaGATAGATAAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGtt            <  2:248349/250‑1 (MQ=255)                                                                                                                                                                                                                         agatagatAAACAGAACTGAAAAAGACAGACTCTCGTCGCCTCAGAGTTGTTGCATTTTGCCATATGTTACCATATACCATTACACATCATCTAAATTAACTCTGGAGGAAACTGTTTTGGCGATTAGTCTTTACAAACTGGCTGTTGCTTTAGGGGTCTTTATTGTTAGCGCTCCTGCCTTTTAGCATGGTCATCACTCACACGGCAAAACCTTAACAGAGGCCGAACAAAAAGCTGCTACTGGTGtttt          >  2:44563/1‑251 (MQ=255)                                                                                                                                                                                                                             agataAGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAACGCTGCTAATGGTGTTTTtgat      >  2:148722/1‑251 (MQ=255)                                                                                                                                                                                                                                 aaGTAGAACTGAAAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCGGCTAATGGTGTTTTTGatgatg  <  1:283702/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               TAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATG  >  NC_000913/2041042‑2041516

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