Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	372,529	G→A	V139M (GTG→ATG)	mhpD →	2‑hydroxypentadienoate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	372,529	0	G	A	100.0%	137.2 / NA	38	V139M (GTG→ATG)	mhpD	2‑hydroxypentadienoate hydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/16);  total (22/16)

CAACTGGGCGTTGATCAACCGGATTTTGGGACGTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGG  >  NC_000913/372325‑372769                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  caacTGGGCGTTGATCAACCGGATTTTGGGACGTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTGCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACTCCATTGAATGTGTAGTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATAc                                                                                                                                                                                                     <  1:194036/251‑1 (MQ=255)      gggCGTTGATCAACCGGATTTTGGGACGTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGt                                                                                                                                                                                                  <  2:281676/249‑1 (MQ=255)            tGATCAACCGGATTTTGGGACGTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGAt                                                                                                                                                                                           >  1:410690/1‑250 (MQ=255)             gATCAACCGGATTTTGGGACGTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATa                                                                                                                                                                                          <  1:250631/250‑1 (MQ=255)                taaCCGTATTTTGGGACGTTATTTACCAACATGTGTGATGGCGATAACTAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACGGATATCACCGTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTTGCAGATAAct                                                                                                                                                                                       <  1:391426/249‑2 (MQ=255)                 aaCCGGATTTTGGGACGTTATTTGCCGACATGTGTTATGGAGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGcc                                                                                                                                                                                     >  2:41588/1‑251 (MQ=255)                      gATTTTGGGACGTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCtg                                                                                                                                                                                >  1:57299/1‑251 (MQ=255)                                 gTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATg                                                                                                                                                                     >  1:229199/1‑251 (MQ=255)                                 gTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATGCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGGCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATg                                                                                                                                                                     <  2:177226/251‑1 (MQ=255)                                  ttATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGAGTATGt                                                                                                                                                                    >  2:251407/1‑251 (MQ=255)                                                    tATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCACCCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGCTTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTCCTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGACCGgcgc                                                                                                                                                   >  1:392295/1‑250 (MQ=255)                                                         cGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTc                                                                                                                                            <  2:410690/251‑1 (MQ=255)                                                                     cATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGACATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCCAAC‑TCCGGCgggggttg                                                                                                                                >  2:254928/1‑247 (MQ=255)                                                                         ccTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGAccc                                                                                                                            >  1:15401/1‑250 (MQ=255)                                                                           tttttCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGGCCGGCGC‑AAAGTCCGGCGGGGtgcgaacctg                                                                                                                          >  2:297822/1‑242 (MQ=255)                                                                               tCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGaaaa                                                                                                                       <  2:29770/250‑1 (MQ=255)                                                                                cccGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑ACCGTACGGGGGGGGTAGGGCTGaaaaa                                                                                                                      >  2:108815/1‑250 (MQ=255)                                                                                cccGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTAGACCTGaaaaa                                                                                                                      >  1:120608/1‑250 (MQ=255)                                                                                          gCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGGCCGGGGGGGGTAGGCCTGGAAAACTGCGccctg                                                                                                            >  2:158701/1‑247 (MQ=255)                                                                                                ccGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGTCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGa                                                                                                     <  1:108815/251‑1 (MQ=255)                                                                                                        aaGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGtaat                                                                                             <  2:120608/251‑1 (MQ=255)                                                                                                           cGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGGGC‑AACGTCCGGCGGGGTGGGAACTGAAAAAATGCGCCACTAAGATGAcgccgaagaa                                                                                           >  2:296154/1‑243 (MQ=255)                                                                                                            ggAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACg                                                                                         <  2:248749/251‑1 (MQ=255)                                                                                                            ggAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACg                                                                                         >  1:248749/1‑251 (MQ=255)                                                                                                                aTCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCGTCGACGAATTGTATAACGCCATTGAATGTGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACTCGTAATAACGAaga                                                                                     <  2:48063/251‑1 (MQ=255)                                                                                                                     aCTGGTGTTGAACAGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATACCGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTAGAACTGAAAAAACGCGCCACGAAGAAGACGCGttaaaaagacgagggtt                                                                                >  1:298507/1‑235 (MQ=255)                                                                                                                        ggTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCt                                                                              <  1:380856/250‑1 (MQ=255)                                                                                                                                                aaCCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGta                                                                                               <  1:43686/209‑1 (MQ=255)                                                                                                                                                    gATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGGCGCGTAATAAAGGAGAGGGTTCTAGCGGGGGCGGCAGCGGATGCCTggggg                                                  >  2:13606/1‑248 (MQ=255)                                                                                                                                                                gacgaATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGc                                     >  1:278276/1‑251 (MQ=255)                                                                                                                                                                gacgaATTGTATAACGCCATTGAATGGGTACTTCCGGCGATGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑CAAGGCCGGCGGGGGTAGAACTGCAAAAATGCGCCATGAAGATGACGCGGAATAAAGAAAAGGTTTTTAGCGGGGGCGGGAGCGAATGCCTGGGGCATCCGCTTCATgg                                     >  2:41823/1‑250 (MQ=255)                                                                                                                                                                 acgaATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTTAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCCTGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCCGGGACATCCGCTTAATGCg                                    >  1:47251/1‑251 (MQ=255)                                                                                                                                                                   gaATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGc                                  >  1:140551/1‑251 (MQ=255)                                                                                                                                                                                     gAATGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AAAGTCCGGCGGGGGTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGGTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGcc                 >  1:355034/1‑249 (MQ=255)                                                                                                                                                                                      aaTGGGTACTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGGTTGACCTGAAAAACTGCGCCCTGAAGAAGACGCGTAATAAAGAAGAGGGTTCTAGCGGGCGCGGCCGCGAAAGCCCGGGGCATCCGCTTAATTCGGCCGTctgggcggcccccca                >  2:325903/1‑239 (MQ=255)                                                                                                                                                                                             aCTTCCGGCGCTTGAAATGGTGGGGATCCGCATTCGCGACTGGTCGATTCAGTTTGTCGGTACCGTGGCAGATAACTCCTCCTGTGTGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGTCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGTAGATGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGTCATCCGCTTAATGCGGCCGTCTGTCTTGCACGCAAAATGGcc         <  1:275734/250‑1 (MQ=255)                                                                                                                                                                                             aCTTCCGGCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGcc         <  2:334973/250‑1 (MQ=255)                                                                                                                                                                                                   ggCGCTGGAAATGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTgg  <  1:108345/251‑1 (MQ=255)                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  CAACTGGGCGTTGATCAACCGGATTTTGGGACGTTATTTGCCGACATGTGTTATGGCGATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGC‑AACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGG  >  NC_000913/372325‑372769

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