Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	791,688	G→A	S123F (TCC→TTC)	galE ←	UDP‑glucose 4‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	791,688	0	G	A	96.0%	145.0 / ‑1.8	51	S123F (TCC→TTC)	galE	UDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/31);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

CGATGTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGATAACAG  >  NC_000913/791443‑791931                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cGATGTATGGCATCAGGTTATGCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTa                                                                                                                                                                                                                                                 <  2:402250/251‑1 (MQ=255)   aTGTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTaaa                                                                                                                                                                                                                                               <  1:165883/251‑1 (MQ=255)   aTGTATGGCATCAGGTTATTCGGAATGCCTTGCGGATATTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTaa                                                                                                                                                                                                                                                >  1:132998/1‑250 (MQ=255)     gTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTACAGTCCGGCTGGGCTTTTTGAAGATCGGTGAGGATCTGTTCCACCATCAGCTTGATTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTaaaa                                                                                                                                                                                                                                              >  1:170904/1‑250 (MQ=255)     gTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCAAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTaaaa                                                                                                                                                                                                                                              >  1:413246/1‑250 (MQ=255)      tATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAAt                                                                                                                                                                                                                                            >  1:43301/1‑251 (MQ=255)              cAGGTTATTGGGAATGCCTTGCGGATCTTCGCCCATAGCGCCCGACGGATGCGCGCTATCCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTGCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGTTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAAGAAAGttt                                                                                                                                                                                                                                     <  2:395724/250‑1 (MQ=255)              cAGGGTATTCGGAATGCCCTGCGGATCTTCGCCCAGATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGttt                                                                                                                                                                                                                                     <  2:26649/250‑1 (MQ=255)                 gTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGa                                                                                                                                                                                                                                 <  1:342386/251‑1 (MQ=255)                    aTTCTGAATGCCTTTCGTATCGTCGCCCATCTCGCCCGACGTATGCGCGCCTACCGGGTTCAAGTAGCGCAGCAGGGCATGGCTCCAGGCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCGTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGt                                                                                                                                                                                                                              <  2:340449/251‑1 (MQ=255)                       cGGAAGGCCGTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTGCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAg                                                                                                                                                                                                                           <  2:445803/251‑1 (MQ=255)                         gAATGCCTTGCGGAGCGTCGCCCATATCGCCCGACGGATGCGCGCCATCCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTTCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTTATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAgcg                                                                                                                                                                                                                         <  2:341997/251‑1 (MQ=255)                             gCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAgcggcg                                                                                                                                                                                                                      <  1:122081/250‑1 (MQ=255)                               cTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGgcgcgc                                                                                                                                                                                                                   >  1:82801/1‑251 (MQ=255)                                      aTCTTCGCCCATATCGCCCGACGGCTGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGgcg                                                                                                                                                                                                             <  2:100485/250‑1 (MQ=255)                                                        cGACGGATGCGCGCCAACCGGGTTAAAGTAGCGCAGCAGGGCACTGCTCCAGTCCGGCTGGGCTTTTTGCAGAACGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTAGGGAAGCTTTCAACGTATGGGATTTTGGGCTGATCGCCATAAACGGTGGCGGAGCAGCTAAAAATAAAGATTTTGACGTTAGCGGCGCGAAGGGCGCCAAACAggagcggaaggc                                                                                                                                                                                          >  2:262892/1‑241 (MQ=255)                                                           cGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGCGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGt                                                                                                                                                                                       >  2:267375/1‑251 (MQ=255)                                                            ggATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCACTGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGCGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATACGGGGCAGGGTGCCGtt                                                                                                                                                                                      >  2:426581/1‑251 (MQ=255)                                                                     ccAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACAttgt‑t                                                                                                                                                                            >  2:397115/1‑251 (MQ=255)                                                                        aCCGGGTTGAAGTAGCGCAGCAGGGCAATGCGCCAGTCCGGATGGGCTTTTTGCAGAGCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGgc                                                                                                                                                                                                   <  1:128140/226‑1 (MQ=255)                                                                             gTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAat                                                                                                                                                                    <  1:321766/251‑1 (MQ=255)                                                                               tGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGAGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAata                                                                                                                                                                   <  1:291884/250‑1 (MQ=255)                                                                                 aaGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGTATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTATCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTc                                                                                                                                                                <  1:149604/251‑1 (MQ=255)                                                                                       cgcAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCgg                                                                                                                                                          <  1:434616/251‑1 (MQ=255)                                                                                         cagcagGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGTTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGAGAAGAGTGCCTTTGCCATTTT‑TGTCGTAATTTTCCAGCGGtt                                                                                                                                                        >  2:408194/1‑251 (MQ=255)                                                                                          agcagGGCAGTGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCATCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTATCTGCGCGCATGGCGATAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGttt                                                                                                                                                       <  2:199023/251‑1 (MQ=255)                                                                                          agcagGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTg                                                                                                                                                                                                                                  >  1:345625/1‑177 (MQ=255)                                                                                          agcagGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTg                                                                                                                                                                                                                                  <  2:345625/177‑1 (MQ=255)                                                                                           gcagGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGtttt                                                                                                                                                      <  2:446509/251‑1 (MQ=255)                                                                                               ggCAATGATCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACAT‑GTGTGGCGGAATATTCCAACGGTTTTTGt                                                                                                                                                   >  2:326319/1‑250 (MQ=255)                                                                                                         cAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGCGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGcc                                                                                                                                        <  1:267375/251‑1 (MQ=255)                                                                                                           gTCCGGCTGGGCTTTTTGCAGATCGGGGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCa                                                                                                                                      <  2:250828/251‑1 (MQ=255)                                                                                                                gCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCt                                                                                                                                 >  2:248092/1‑251 (MQ=255)                                                                                                                     gCTTTTTGCAGATCGGTGAGGATCTGTTCCGCCATCAGCGTGCTTTTGACGTAAGGGCTTGGCGGTGGGCCGGTCGGGAAGGTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATGGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAg                                                                                                                            <  2:394405/251‑1 (MQ=255)                                                                                                                         tttGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCg                                                                                                                        >  2:133146/1‑251 (MQ=255)                                                                                                                                      gAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATc                                                                                                            >  1:457181/1‑250 (MQ=255)                                                                                                                                              gTTCCACCATCAGCTAGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATAAGGCGCCGAGTGCCGTTGAAATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAGGTGGATCCAGGTGTCg                                                                                                   >  2:421946/1‑251 (MQ=255)                                                                                                                                                 ccaccaTCAGCTTGCTGTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATa                                                                                                <  1:421946/251‑1 (MQ=255)                                                                                                                                                  caccaTCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAg                                                                                               <  1:209731/251‑1 (MQ=255)                                                                                                                                                                     ccGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGAtctc                                                                            <  2:409219/251‑1 (MQ=255)                                                                                                                                                                            ggCTTTGCGGTTTGCCGGTCGGGAAGCTGTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGtcatc                                                                     <  2:411281/251‑1 (MQ=255)                                                                                                                                                                              cTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAAGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCAc                                                                                                                                     >  1:70334/1‑185 (MQ=255)                                                                                                                                                                              cTTTGCGGTGTGCCGGTCGGGAAGCTTGCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAAGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATGGT‑TGTCGTAATATTCCAGCGGTGTTTGTACCGATTCGCCCAc                                                                                                                                     <  2:70334/185‑1 (MQ=255)                                                                                                                                                                                                 ggAAGGTTTCAACGCATGGAAATTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGGTTTTGACGGTAGCGGCGCGCATGGCGCTAATCAGGCGCAGGTTCCCGTTTACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTTCCGATTCGCCCTCGGCTTTCATCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTTCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCg                                                <  2:345105/251‑1 (MQ=255)                                                                                                                                                                                                            aCGTATTGAATTTTGGGCTGAGCGCCATTAACGGGGGCGGAGAAGCTAAAAATAAAGTTTTTGCCGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCTCCGGCAACaaa                                     <  2:420654/251‑1 (MQ=255)                                                                                                                                                                                                                    aaTTTTGGGCTGATCGCCATAAACGGTAGCGGAGAAGCTAAAAATAAAGTTTTAGACGGTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGGGCCGTTGAAATAGT‑TGTCGTAATCTACAAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCACGGCGCAGTGGATCACGGTGGCGCTAGCGTGCTCGTGCAGGATCTCGGGCATCAACGCTTCGATAAGGATATCGCCTTCATCAAACGTTGGAt                             >  1:290911/1‑251 (MQ=255)                                                                                                                                                                                                                         tGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGttt                         <  1:40037/250‑1 (MQ=255)                                                                                                                                                                                                                            gCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTgccg                     >  1:156529/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTaaa               <  1:248092/251‑1 (MQ=255)                                                                                                                                                                                                                                          acgggtgCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGAt       <  2:457181/245‑1 (MQ=255)                                                                                                                                                                                                                                               ggCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCATCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGATAACAg  <  2:410005/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CGATGTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGT‑TGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGATAACAG  >  NC_000913/791443‑791931

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