Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	371,179	Δ1 bp	coding (903/945 nt)	mhpB →	3‑carboxyethylcatechol 2,3‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	371,178	0	G	.	100.0%	169.8 / NA	40	coding (902/945 nt)	mhpB	3‑carboxyethylcatechol 2,3‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (16/24);  total (16/24)

CCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCT  >  NC_000913/371040‑371309                                                                                                                                           |                                                                                                                                    ccGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTga                                                                                                                                    >  1:180051/1‑139 (MQ=255)  cGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTgat                                                                                                                                   >  1:1148826/1‑138 (MQ=255)  cGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTgat                                                                                                                                   >  2:368790/1‑138 (MQ=255)        tAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCgg                                                                                                                            <  1:432518/140‑1 (MQ=255)           cGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGAtt                                                                                                                         >  1:673032/1‑140 (MQ=255)           cGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGAtt                                                                                                                         >  1:62809/1‑140 (MQ=255)           cGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGAtt                                                                                                                         >  2:224070/1‑140 (MQ=255)             aaGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTg                                                                                                                       >  2:234182/1‑140 (MQ=255)                agCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCt                                                                                                                    <  1:451933/140‑1 (MQ=255)                     tCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCNCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTa                                                                                                               <  1:267734/140‑1 (MQ=255)                     tCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTa                                                                                                               <  2:62809/140‑1 (MQ=255)                     tCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTa                                                                                                               <  2:696272/140‑1 (MQ=255)                     tCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTa                                                                                                               <  1:540590/140‑1 (MQ=255)                       cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAg                                                                                                             <  1:234182/140‑1 (MQ=255)                       cGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAg                                                                                                             <  1:963358/140‑1 (MQ=255)                            ttGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCa                                                                                                        >  2:997199/1‑140 (MQ=255)                            ttGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCa                                                                                                        <  1:881941/140‑1 (MQ=255)                               ccGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGaa                                                                                                     <  2:539086/140‑1 (MQ=255)                               ccGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGaa                                                                                                     >  1:882563/1‑140 (MQ=255)                                     aGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAgaga                                                                                               <  1:156707/140‑1 (MQ=255)                                        cGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACt                                                                                            <  2:1140008/140‑1 (MQ=255)                                          acacATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGa                                                                                          <  2:276570/140‑1 (MQ=255)                                           cacaTGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGaa                                                                                         <  1:224070/140‑1 (MQ=255)                                             caTGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAata                                                                                       >  2:274729/1‑140 (MQ=255)                                                     cAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGgag                                                                               <  2:631934/140‑1 (MQ=255)                                                              ggTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGatgatg                                                                      >  2:1121608/1‑140 (MQ=255)                                                              ggTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGatgatg                                                                      >  2:370202/1‑140 (MQ=255)                                                                                tATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAAcc                                                    <  2:180051/140‑1 (MQ=255)                                                                                 aTTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCg                                                   <  1:274729/140‑1 (MQ=255)                                                                                             ggCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGccacca                                       <  2:1148826/140‑1 (MQ=255)                                                                                              gCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAg                                      <  2:590906/140‑1 (MQ=255)                                                                                                     gCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGtttt                               <  2:739524/140‑1 (MQ=255)                                                                                                        tAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTg                            <  1:539138/140‑1 (MQ=255)                                                                                                             aaGGGCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGt                       >  1:1143025/1‑140 (MQ=255)                                                                                                               gggCGTTATTACCGCCCAATCCCGGAGT‑GATTGTCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTag                     >  1:597927/1‑140 (MQ=255)                                                                                                               gggCGTTATTACCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTag                     >  1:808607/1‑140 (MQ=255)                                                                                                                   gTTATTCCCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTagaaga                 <  2:70034/140‑1 (MQ=255)                                                                                                                    ttattaCCGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAgaagaa                >  1:899808/1‑140 (MQ=255)                                                                                                                           cGCCCAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTa         >  2:188066/1‑140 (MQ=255)                                                                                                                             cccAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTaaa       >  1:300678/1‑140 (MQ=255)                                                                                                                              ccAATCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTaaaa      >  1:668699/1‑140 (MQ=255)                                                                                                                                  tCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCt  <  1:370202/140‑1 (MQ=255)                                                                                                                                  tCCCGGAGT‑GATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCt  <  1:1121608/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                    CCGTCAGTAACGAAGAGCTTTCCGCCATTGCCGGAAAGTCGACACATGAAATCAAAACCTGGGTCGCCGCTTTTGCCGCTATTTCTGCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAACGCT  >  NC_000913/371040‑371309

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