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.4%	182.1 / ‑2.1	55	S123F (TCC→TTC)	galE	UDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (25/28);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

ATGTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGATA  >  NC_000913/791445‑791927                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                aTGTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGCGAAGCTaa                                                                                                                                                                                                                                           >  2:542719/1‑250 (MQ=255)      tGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATa                                                                                                                                                                                                                                      >  1:223970/1‑250 (MQ=255)                  aTTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGATTTTGACGt                                                                                                                                                                                                                         >  2:420764/1‑251 (MQ=255)                     cGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGAGCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTa                                                                                                                                                                                                                       <  1:498925/250‑1 (MQ=255)                      ggAAGGCCCTGCGGATCTTTCCCCATTACGGCCCACGGATGCGCGCCAACCGGGGTGAAGTAGCGCAGCAGGGGCATGCTCCAGGCCCGCTGGGCTTTTTTCAGATCGGGGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGGTAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAg                                                                                                                                                                                                                      <  2:40843/250‑1 (MQ=255)                       gAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAgcg                                                                                                                                                                                                                    <  2:83593/251‑1 (MQ=255)                        aaTGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAgcgg                                                                                                                                                                                                                   <  2:20666/251‑1 (MQ=255)                        aaTGCCTTGCGGATCTTCGCCCATATCGCACGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGCGAAGCTAAAAATAAAGTTTTTGACGTTAgcgg                                                                                                                                                                                                                   >  2:74697/1‑251 (MQ=255)                         aTGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAgcggc                                                                                                                                                                                                                  >  2:83174/1‑251 (MQ=255)                           gCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCACCGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAgcggcgg                                                                                                                                                                                                                >  2:127822/1‑250 (MQ=255)                                         cGCCCATATCGCCCGATGGATGCTCGCCAACCGGTTTGGAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAAt                                                                                                                                                                                                   <  1:305626/250‑1 (MQ=255)                                         cGCCCATATCGCCCGACGGATGCGCGCCAACCGGGGTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAAt                                                                                                                                                                                                   <  1:539278/250‑1 (MQ=255)                                          gCCCATATCGCCCGACGGCTGCGCGCCAACCGGGTTGAAGTAGCGCAGCACGGCAATGCTCCAGTCCGACTGGGCTTTTTGCAGATCGGAGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCCTGGCCCTAATc                                                                                                                                                                                                  >  2:370013/1‑250 (MQ=255)                                            ccATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGAGCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGa                                                                                                                                                                                                                            <  1:60046/222‑1 (MQ=255)                                                                   ccAACCGGGTTGAAGTAGCGCAGCAGCGCAATGCTCCAGTCCGGCTGGGCTTTCTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGACGCAGAGTGCCGTTGATAttgtt                                                                                                                                                                        >  1:214825/1‑251 (MQ=255)                                                                     aaCCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATtgttgt                                                                                                                                                                      <  1:475392/251‑1 (MQ=255)                                                                      aCCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGTCATTGTTGTc                                                                                                                                                                     >  1:310097/1‑251 (MQ=255)                                                                        cgtgTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGAGCGGTGAGGATCTGTTCCACCATCAGCGTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGt                                                                                                                                                                   <  2:55138/248‑1 (MQ=255)                                                                              gAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATAtt                                                                                                                                                             <  1:402246/251‑1 (MQ=255)                                                                                      gcagcaGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGt                                                                                                                                                     <  1:3866/251‑1 (MQ=255)                                                                                           agGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTg                                                                                                                                                >  2:104885/1‑251 (MQ=255)                                                                                              gCAATGCTCCGGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCTGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTT‑GTAcc                                                                                                                                            >  1:297829/1‑251 (MQ=255)                                                                                               cAATGCTCCAGGCCGGCTGGGCTTTTTGCAGAGCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTAcc                                                                                                                                            <  1:235373/251‑1 (MQ=255)                                                                                                                 gggCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCtt                                                                                                                            <  2:259059/249‑1 (MQ=255)                                                                                                                            aGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGa                                                                                                                >  1:175962/1‑250 (MQ=255)                                                                                                                              aTCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGGTTCGCCCACGGCTTTCAGCCCGGCGAAgg                                                                                                             >  2:23100/1‑250 (MQ=255)                                                                                                                              aTCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAgg                                                                                                             >  2:269515/1‑250 (MQ=255)                                                                                                                                          cggTTCCACCATCAGCTTGCGTTTGCAGTAAGGGCTTTGCGGGGGGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGATTGCCGTTGACATTGTTGTCGTAATAGTCCGGCGGTTTTTGTTCCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGgtgt                                                                                                 <  2:454238/249‑1 (MQ=255)                                                                                                                                            gTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGgtgt                                                                                                 <  1:33011/249‑1 (MQ=255)                                                                                                                                                            gcttTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGCAGTGGATc                                                                                                        >  2:285532/1‑226 (MQ=255)                                                                                                                                                            gcttTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATc                                                                                                        <  1:285532/226‑1 (MQ=255)                                                                                                                                                                     gTAAGGGCTTTGCGGAGTGCCGGTCGGGCAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCCTGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGGGTCGCTAGCTGGATCGTGCAGGATCTCgg                                                                      >  2:249246/1‑251 (MQ=255)                                                                                                                                                                      tAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGt                                                                     >  2:184641/1‑251 (MQ=255)                                                                                                                                                                               tGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACg                                                             <  1:139617/250‑1 (MQ=255)                                                                                                                                                                                  ggTGTGCCGTGCGGGAAGCTGTCAACGTATGGAATTTTGGGCTGATCGCCAGAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGGACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGGCGATAGCGTGATCGTGCAGGATCTCTGTCATCAACGCTTc                                                         <  2:223970/251‑1 (MQ=255)                                                                                                                                                                                     gtgCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGAGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGACGAAGTGGATCACGGTGTCGATAACGTGCTCGTGCAGGATCTCGGTCATCAACGCTTCGtt                                                      >  2:419776/1‑251 (MQ=255)                                                                                                                                                                                       gcctgTCGGGAAGCTTTCAACGTATGGAAGTTTGGGCTGATCGCCATAAACGGGGGCGGAGAAGCTAAAAATAAAGTTTGTGACGTTAGCGGCGCGCATGGCTCTAATCAGTCGCAGAGTGCCGCTGACATTGTTGTCGTAATATTCCAGCGGTGTTTGTACCGATTCGCCCACGGCTTGCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCGACGCTTCGTTa                                                     <  2:83720/246‑1 (MQ=255)                                                                                                                                                                                           gTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGAGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGCGGGCAGTTGGAATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGTTCACGGTGTCGATAGCGTGATCGTGCACGAGCTCGGTCATCAACGCTTCGTTCCGac                                                 >  2:547621/1‑249 (MQ=255)                                                                                                                                                                                             tgggAAGCTTTCAACGTATGGAATTTTGGGCTGATCGTCATTAACGGTGGCGGAGAAGCTAAAAATAAAGTATTTGACGTCTTCGTCGCGCAGGGGGCTAATCGGTCTCAGGGTGCCGTGGACATTGCTGTCGTAATATTCCTGCGGGTTGTGTGACGATTCGCCTTCGGCGTGCAGCCCGGCGAAGTGGATCACGGTATCGATAGCGTGAGCGTTCAGGATGTTGGCCTCCAACGCTGCGTTGCGAata                                               <  2:361459/249‑1 (MQ=255)                                                                                                                                                                                             cGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAatat                                              >  1:351837/1‑251 (MQ=255)                                                                                                                                                                                             cGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAatat                                              >  1:439009/1‑251 (MQ=255)                                                                                                                                                                                                 aaGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGCTAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTCCGAATATCGcc                                          >  2:80006/1‑251 (MQ=255)                                                                                                                                                                                                          aCGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACaaa                                 <  1:23100/251‑1 (MQ=255)                                                                                                                                                                                                                  aaTTTTTGGCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAATTTTTTGACGTTTGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTCGTACCGATTCGCCCACGGCTTTCAGCCCCGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTTCAGGATCTCGGGCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGAt                         <  2:515603/251‑1 (MQ=255)                                                                                                                                                                                                                        gggCCGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTgc                   <  2:112382/251‑1 (MQ=255)                                                                                                                                                                                                                         ggCTGATCGCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTgcc                  >  1:419096/1‑251 (MQ=255)                                                                                                                                                                                                                                 gCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTaaa           >  1:226612/1‑250 (MQ=255)                                                                                                                                                                                                                                 gCCATAAACGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAAc          <  1:184641/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaaCGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCg     <  1:542719/251‑1 (MQ=255)                                                                                                                                                                                                                                       aaCGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCg     >  1:236655/1‑250 (MQ=255)                                                                                                                                                                                                                                       aaCGGGGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGCCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGa    <  2:459375/251‑1 (MQ=255)                                                                                                                                                                                                                                        aCGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCTTTGACAt                                                                                                                                                                            <  2:498333/82‑1 (MQ=255)                                                                                                                                                                                                                                        aCGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACAt                                                                                                                                                                            >  1:498333/1‑82 (MQ=255)                                                                                                                                                                                                                                         cGGTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGATa  <  1:410653/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggTGGCGGAGAAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGATa  >  2:263928/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                ATGTATGGCATCAGGTTATTCGGAATGCCTTGCGGATCTTCGCCCATATCGCCCGACGGATGCGCGCCAACCGGGTTGAAGTAGCGCAGCAGGGCAATGCTCCAGTCCGGCTGGGCTTTTTGCAGATCGGTGAGGATCTGTTCCACCATCAGCTTGCTTTTGCCGTAAGGGCTTTGCGGTGTGCCGGTCGGGAAGCTTTCAACGTATGGAATTTTGGGCTGATCGCCATAAACGGTGGCGGAGGAGCTAAAAATAAAGTTTTTGACGTTAGCGGCGCGCATGGCGCTAATCAGGCGCAGAGTGCCGTTGACATTGTTGTCGTAATATTCCAGCGGTTTTTGTACCGATTCGCCCACGGCTTTCAGCCCGGCGAAGTGGATCACGGTGTCGATAGCGTGATCGTGCAGGATCTCGGTCATCAACGCTTCGTTACGAATATCGCCTTCAACAAACGTTGGATGTTTGCCGCCTAAACGCTCGATA  >  NC_000913/791445‑791927

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