Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	370,874	A→G	S200G (AGC→GGC) ‡	mhpB →	3‑carboxyethylcatechol 2,3‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	370,874	0	A	G	100.0%	132.2 / NA	40	S200G (AGC→GGC) ‡	mhpB	3‑carboxyethylcatechol 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/16);  total (24/16)

TGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGC  >  NC_000913/370737‑371009                                                                                                                                              |                                                                                                                                         tGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGc                                                                                                                                        >  1:20515/1‑140 (MQ=255)    gTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCggg                                                                                                                                     <  2:863557/140‑1 (MQ=255)     tGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGa                                                                                                                                    <  2:691556/140‑1 (MQ=255)             ttGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑TTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑g                                                                                                                                <  2:27672/135‑1 (MQ=255)             ttGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑g                                                                                                                                >  1:27672/1‑135 (MQ=255)              tGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCTTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTAc                                                                                                                          >  1:1642361/1‑140 (MQ=255)                      tCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGt                                                                                                                  >  1:1283987/1‑140 (MQ=255)                               ctctCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGaa                                                                                                         >  1:1584920/1‑140 (MQ=255)                               ctctCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGaa                                                                                                         >  1:899085/1‑140 (MQ=255)                                      tAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTc                                                                                                  >  1:463187/1‑140 (MQ=255)                                         ccgggtgtTGTCCCTGGGTCCGGGGGGG‑CTTC‑CCCACCGCCCGCGGGCCCCCGAC‑CGGGAAAAA‑GCCGATCCCCATAGGCGCGCCTGTGTGTGGGGGGGGGGAAAA‑AATTTACCCGCCAGTGAGCGCGAATTGCGTCAGc                                                                                               <  2:1450982/136‑1 (MQ=255)                                          cgcgTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCa                                                                                              >  2:1240420/1‑140 (MQ=255)                                          cgcgTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCa                                                                                              >  1:900756/1‑140 (MQ=255)                                            cgTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAAc                                                                                            >  1:1227250/1‑140 (MQ=255)                                              tgtttttCCGGGGTCCGGGTGGT‑CTTC‑CCCACCGCCCGCGGGCGCCGGAA‑CGGGCAAAA‑GCCGACGCCCATATGCGCGCCCGTTTGTGGGGGGGGGGGAAA‑GTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACgg                                                                                          <  2:1283987/135‑1 (MQ=255)                                              tGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACgg                                                                                          >  1:85095/1‑140 (MQ=255)                                              tGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACgg                                                                                          >  1:997886/1‑140 (MQ=255)                                              tGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACgg                                                                                          >  2:305457/1‑140 (MQ=255)                                               ggctgtCCCGGGGTCCGGGGGG‑TCTT‑CCCCACCGCCGGCGGGGGCCGGA‑GCGGGGAAAGGCCGTCGCCAATAGGCGCGCCCTTTGGTGGGGGGGGGGAAAA‑GTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACgg                                                                                          <  2:138867/135‑1 (MQ=255)                                                   ttCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGAt                                                                                     >  1:1450982/1‑140 (MQ=255)                                                      gggggTGCGGGGGGG‑TTTC‑CCCACCGGCCGCGGGTGCCCGAG‑CGGGCGAAA‑GCCGATCCCAATAGGCGCGACCGTTTGTGGGGGGGGGGGAAA‑GTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAg                                                                                  <  2:301080/138‑1 (MQ=255)                                                       gggggTCCGGGGGGCCTTC‑CCCACCGGCCGCGGGTGCCCGAG‑GGGGAGAAA‑GCCGATGCCCATATGCGCGCCCGTTTGTGGGGGGGGGGGAAA‑TATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAg                                                                                  <  2:463187/139‑1 (MQ=255)                                                         gggtCCGGGGGG‑CTTTCCCCACCGCCCGCGGGGGCCCGAA‑CGGGGAAAG‑GCGGATGCCCATAGGCGCGCCCGTCTGTGGGGGGGGGGAAAA‑GTTTTCCCCGCCAGTGAGCGCGAATTGCATCAGCAACGGGTGATTAgcg                                                                                <  2:900756/137‑1 (MQ=255)                                                           ttCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGcgccg                                                                             >  1:723153/1‑140 (MQ=255)                                                             ccGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCt                                                                           >  1:301080/1‑140 (MQ=255)                                                              ggggggg‑CTTTCCCCACCCGCCGCGGGTGCCCGAA‑GGGGCAAAA‑GCCGCTCCCAATATGCGCGCCCGTGTGGGGGGGGGGGGAAAA‑TTTTCACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCt                                                                           <  2:506581/136‑1 (MQ=255)                                                                            ccaccggccgcgGGTGCCCGAA‑CGGGAGAAG‑GCCGACGCCAATATGCGGGCGGGTGTGTGGGGGGGGGGAAAA‑TTTTTGCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTttgttg                                                              <  2:1642361/134‑1 (MQ=255)                                                                                cggccgccGG‑GCCCGAAGCTGGCAAAG‑GCGGTGGCCTATGGGCGCGGCCGTTTGTGGGGGGGGGGAAAA‑TTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGa                                                          <  1:305457/138‑1 (MQ=255)                                                                                      ccgGTGCCCGAA‑GGGGCAAAA‑GCCGATGCCCATAGGCGCGACCGTTTGTGGGGGGGGGGAAAA‑TTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCagag                                                    <  2:899085/140‑1 (MQ=255)                                                                                                aa‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCAt                                          >  2:977174/1‑140 (MQ=255)                                                                                                aa‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCAt                                          >  2:48960/1‑140 (MQ=255)                                                                                                    gggCAAAA‑GCCGAGCCCCATAGGCGCGACTGTTTGTGGGGGGGGGGAAAA‑GTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCg                                       <  2:1227250/139‑1 (MQ=255)                                                                                                           aggcggATCCCCATGGGCGCGCCCGTGTGTGGGGGGGGGGAAAA‑TTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAAc                                 <  2:714998/135‑1 (MQ=255)                                                                                                                gATGCCAATAGGCGCGCCCGTGTGTGGGGGGGGGGAAAA‑GTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAt                            <  2:1584920/140‑1 (MQ=255)                                                                                                                 aTGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAtt                           >  1:714998/1‑140 (MQ=255)                                                                                                                 aTGCCCATATGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGAtt                           >  2:2495/1‑140 (MQ=255)                                                                                                                        tatGCGCGACCGTCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATa                    >  2:1556160/1‑140 (MQ=255)                                                                                                                          tGCGCGCCCGTTTGTGGGGGGGGGGGAAA‑GTTTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAAc                  <  1:1038212/140‑1 (MQ=255)                                                                                                                                    tCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGa        >  1:524476/1‑140 (MQ=255)                                                                                                                                    tCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGa        >  1:1290967/1‑140 (MQ=255)                                                                                                                                    tCTGTTGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGa        >  2:556350/1‑140 (MQ=255)                                                                                                                                      ttgtGGGGGGGGGGAAATTATTTCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGAc       <  2:742714/137‑1 (MQ=255)                                                                                                                                         tGGGGGGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTg   >  1:742714/1‑140 (MQ=255)                                                                                                                                          ggggggggGAAAA‑GTTTCCCCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGc  <  2:1607155/140‑1 (MQ=255)                                                                                                                                          ggggggggGAAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGc  <  1:1556160/140‑1 (MQ=255)                                                                                                                                              |                                                                                                                                         TGGGTGAAGCCATTGGACGTTTCACCAGCACTCTCAATAAACGCGTGCTGTTCCTGGGTTCCGGTGGG‑CTTT‑CCCATCAGCCGCCGGTGCCCGAA‑CTGGCGAAA‑GCCGATGCCCATATGCGCGACCGTCTGTTGGGGAGCGGGAAA‑GATTTACCCGCCAGTGAGCGCGAATTGCGTCAGCAACGGGTGATTAGCGCCGCTGAGAAGTTTGTTGAGGATCAGAGAACGCTGCATCCGCTCAACCCGATTTGGGATAACCAGTTCATGACTTTGC  >  NC_000913/370737‑371009

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