Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,894,654	G→A	M1M (ATG→ATA) †	chrR →	chromate reductase, Class I, flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,894,654	0	G	A	100.0%	116.9 / NA	35	M1M (ATG→ATA) †	chrR	chromate reductase, Class I, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

ACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAAGTCAATGCGTTACCATCCATTGCCGACA  >  NC_000913/3894522‑3894787                                                                                                                                     |                                                                                                                                      acacGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGaa                                                                                                                                 >  1:810653/1‑139 (MQ=255)          gAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAg                                                                                                                        <  1:708760/139‑1 (MQ=255)            aaaGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAggt                                                                                                                      <  2:698808/139‑1 (MQ=255)                 cTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGgtta                                                                                                                 <  2:810653/139‑1 (MQ=255)                 cTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGgtta                                                                                                                 <  1:485314/139‑1 (MQ=255)                 cTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGgtta                                                                                                                 <  1:840653/139‑1 (MQ=255)                  tGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGC‑TTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGgttacg                                                                                                               <  2:339386/139‑1 (MQ=255)                            tCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATTTGAAAAATTGCAGGTGGTTACGTTACTgggg                                                                                                      <  1:205064/139‑1 (MQ=255)                                     cgcCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTgcgc                                                                                             >  2:66462/1‑139 (MQ=255)                                     cgcCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTgcgc                                                                                             >  2:129051/1‑139 (MQ=255)                                         aacaaTCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAg                                                                                         <  1:129051/139‑1 (MQ=255)                                           caaTCCTACCAGCCTGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGc                                                                                       <  1:66462/139‑1 (MQ=255)                                           caaTCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGc                                                                                       <  1:162060/139‑1 (MQ=255)                                                  aCCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTa                                                                                >  1:133866/1‑139 (MQ=255)                                                        cGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCa                                                                          <  2:518800/139‑1 (MQ=255)                                                                 cGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCAc                                                                 >  1:110864/1‑139 (MQ=255)                                                                           aaCTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGcc                                                       >  2:994446/1‑139 (MQ=255)                                                                                tctACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGaaaa                                                  >  2:1054389/1‑139 (MQ=255)                                                                                  tACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAAtt                                                >  2:295021/1‑139 (MQ=255)                                                                                       tGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTcc                                           >  1:292669/1‑139 (MQ=255)                                                                                       tGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTcc                                           >  1:919427/1‑139 (MQ=255)                                                                                       tGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTcc                                           >  1:138230/1‑139 (MQ=255)                                                                                        gAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCg                                          <  2:743093/139‑1 (MQ=255)                                                                                         aaGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTc                                                                                     <  2:882977/95‑1 (MQ=255)                                                                                         aaGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTc                                                                                     >  1:882977/1‑95 (MQ=255)                                                                                         aaGGCTTTCTCGCCAAATTGATAGACAACCACAGGAGTCATTATATCTGAAAAATTGCAGGTGGTTACGTAAGTGTGGAGCCTGAGCAAAGGCTCATTTAATGGCGTGGTTGCATCTACCCTGCCGAAAATTGCTCCga                                         >  1:130397/1‑138 (MQ=255)                                                                                          aGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGc                                        <  2:133866/139‑1 (MQ=255)                                                                                                             aTAGACAACCACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAAGTCAATGCg                     <  1:1054389/139‑1 (MQ=255)                                                                                                                     ccACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCa                                                                          <  2:169415/78‑1 (MQ=255)                                                                                                                     ccACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCa                                                                          >  1:169415/1‑78 (MQ=255)                                                                                                                     ccACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATgg                                <  2:723856/120‑1 (MQ=255)                                                                                                                     ccACAGGAGTCATCATATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATgg                                >  1:723856/1‑120 (MQ=255)                                                                                                                                atcatATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGTCTGCGCAAAGGCTCATTTAATGGCATGGt                                                                      <  2:265299/71‑1 (MQ=255)                                                                                                                                atcatATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGTCTGCGCAAAGGCTCATTTAATGGCATGGt                                                                      >  1:265299/1‑71 (MQ=255)                                                                                                                                atcatATCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAAGTCAATGCGTTACCATCCATTGCCGACa  >  2:998364/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                      ACACGGCTGGAAAAGCCTGCCGGATATTCACAGCCGCGCCAACAATCCTACCAGCCGGATGATGCGTTTTGCCCAACTCTCTACCGTGAAGGCTTTTTCGCCAAATTGATAGACAACCACAGGAGTCATCATGTCTGAAAAATTGCAGGTGGTTACGTTACTGGGGAGCCTGCGCAAAGGCTCATTTAATGGCATGGTTGCACGTACCCTGCCGAAAATTGCTCCGGCGAGCATGGAAGTCAATGCGTTACCATCCATTGCCGACA  >  NC_000913/3894522‑3894787

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