Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,730,927	0	T	G	55.1%	2.9 / 108.3	69	intergenic (‑162/‑204)	xylA/xylF	D‑xylose isomerase/D‑xylose transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (14/17);  new base G (18/20);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

AATCAGGTAATGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTAAAAATTGGTTACGTT  >  NC_000913/3730778‑3731073                                                                                                                                                      |                                                                                                                                                   aaTCAGGTAATGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAAtt                                                                                                                                                    <  2:853241/150‑1 (MQ=255)    cAGGTAATGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTc                                                                                                                                                <  1:905043/151‑1 (MQ=255)    cAGGTAATGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTc                                                                                                                                                >  1:442840/1‑151 (MQ=255)      ggTAATGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCC‑TTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGaa                                                                                                                                             <  1:885322/151‑1 (MQ=255)         aaTGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAag                                                                                                                                            <  1:740242/150‑1 (MQ=255)          aTGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACTCTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAagag                                                                                                                                          <  1:419808/151‑1 (MQ=255)          aTGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAaga                                                                                                                                           <  1:732068/150‑1 (MQ=255)            gCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAgaga                                                                                                                                         <  1:168472/150‑1 (MQ=255)              cgcgGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGaaaa                                                                                                                                      <  2:718818/151‑1 (MQ=255)              cgcgGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGaaa                                                                                                                                       >  1:881819/1‑150 (MQ=255)                cgGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAAt                                                                                                                                    >  2:714606/1‑151 (MQ=255)                    tGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCa                                                                                                                                 >  1:468792/1‑150 (MQ=255)                     gatggatgATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCaa                                                                                                                                <  2:491770/150‑1 (MQ=255)                          atgatgTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTa                                                                                                                          <  1:113449/151‑1 (MQ=255)                           tgatgTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTAc                                                                                                                         >  1:187189/1‑151 (MQ=255)                                       atTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGcaac                                                                                                              >  1:648150/1‑150 (MQ=255)                                        ttGGGCCCTCCCTTTCATTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGCGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGcaaca                                                                                                             >  1:855189/1‑150 (MQ=255)                                         tGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGcaacaa                                                                                                            <  1:284536/150‑1 (MQ=255)                                         tGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGcaacaa                                                                                                            >  1:815233/1‑150 (MQ=255)                                              aCTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTa                                                                                                       >  1:751903/1‑150 (MQ=255)                                                  ccTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGAtc                                                                                                  <  2:558931/151‑1 (MQ=255)                                                    tttCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGAtctc                                                                                                >  1:155438/1‑151 (MQ=255)                                                      tCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCg                                                                                               >  2:759168/1‑150 (MQ=255)                                                       cAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGt                                                                                             <  1:847404/151‑1 (MQ=255)                                                       cAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGt                                                                                             <  2:631483/151‑1 (MQ=255)                                                       cAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGt                                                                                             <  2:349714/151‑1 (MQ=255)                                                          ttGCTCAATTATGTTATTTCACACTGCTATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCAt                                                                                          >  2:201275/1‑151 (MQ=255)                                                                   tATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCa                                                                                    >  2:566262/1‑148 (MQ=255)                                                                    aTGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGa                                                                                 >  2:733120/1‑150 (MQ=255)                                                                      gTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAAt                                                                              >  2:717643/1‑151 (MQ=255)                                                                      gTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAAt                                                                              >  2:113278/1‑151 (MQ=255)                                                                         aTTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATaa                                                                            <  2:416236/150‑1 (MQ=255)                                                                          tttCACACTGCTATTGAGATCATTCAAAAGTGTGCGCTCGCTCGCAAAATAAAATGGGATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAAt                                                                              >  2:188871/1‑147 (MQ=255)                                                                            tCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAAcc                                                                         <  1:802126/150‑1 (MQ=255)                                                                                acTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCaaaaa                                                                    <  1:703984/151‑1 (MQ=255)                                                                                acTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCaaaaa                                                                    <  2:46691/151‑1 (MQ=255)                                                                                   gCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAAt                                                                   <  1:226756/149‑1 (MQ=255)                                                                                   gCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATc                                                                  <  2:458544/150‑1 (MQ=255)                                                                                     tATTGAGATCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGt                                                                >  2:664046/1‑150 (MQ=255)                                                                                      aTTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGt                                                                <  2:218579/149‑1 (MQ=255)                                                                                      aTTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTaa                                                              <  1:894484/151‑1 (MQ=255)                                                                                        tGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATc                                                            >  2:664520/1‑151 (MQ=255)                                                                                         gagaTCATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATc                                                            <  2:31086/150‑1 (MQ=255)                                                                                         gagaTAATTCGCAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATc                                                            <  2:709830/150‑1 (MQ=255)                                                                                          agaTAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCg                                                           >  1:469551/1‑150 (MQ=255)                                                                                             tAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGaa                                                         >  2:442229/1‑149 (MQ=255)                                                                                              aaTTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGaa                                                         >  1:399860/1‑148 (MQ=255)                                                                                                   acaAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATaaaa                                                 <  1:714586/151‑1 (MQ=255)                                                                                                    caAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATaaaaa                                                <  1:632674/151‑1 (MQ=255)                                                                                                     aaGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAAt                                               >  2:214073/1‑151 (MQ=255)                                                                                                       gtgtGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCt                                             <  2:61149/151‑1 (MQ=255)                                                                                                               tcgctcgcAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGt                                     >  2:803336/1‑151 (MQ=255)                                                                                                                cgctcgcAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGtt                                    >  2:192431/1‑151 (MQ=255)                                                                                                                cgctcgcAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTg                                      <  1:759149/149‑1 (MQ=255)                                                                                                                 gctcgcAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGttt                                   >  2:497096/1‑151 (MQ=255)                                                                                                                 gctcgcAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGttt                                   <  1:272229/151‑1 (MQ=255)                                                                                                                 gctcgcAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGttt                                   >  1:234492/1‑151 (MQ=255)                                                                                                                    cgcAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTcc                                <  2:492184/151‑1 (MQ=255)                                                                                                                    cgcAAAATAAAATGGAATGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTcc                                <  2:151703/151‑1 (MQ=255)                                                                                                                     gcAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTccc                               >  1:598263/1‑151 (MQ=255)                                                                                                                              aatGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGt                      <  1:155042/151‑1 (MQ=255)                                                                                                                                 ggAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGc                   <  1:526127/151‑1 (MQ=255)                                                                                                                                   aaTGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCt                  >  1:63544/1‑150 (MQ=255)                                                                                                                                   aaTGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCt                  >  2:738120/1‑150 (MQ=255)                                                                                                                                   aaTGATGAAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTa                 <  1:423341/151‑1 (MQ=255)                                                                                                                                     tgatgaAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTaa                >  1:677215/1‑150 (MQ=255)                                                                                                                                      gatgaAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTaaa               <  1:76308/150‑1 (MQ=255)                                                                                                                                       atgaAACTGGGTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTaaaa              >  2:142446/1‑150 (MQ=255)                                                                                                                                               gggTAATGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTAAAAATTGGtt       <  1:665368/149‑1 (MQ=255)                                                                                                                                                   aaTGCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTAAAAATTGGTTACGtt  <  2:112326/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   AATCAGGTAATGCCGCGGGTGATGGATGATGTCGTAATATTGGGCACTCCCTTTCAGTTGCTCAATTATGTTATTTCACACTGCTATTGAGATAATTCACAAGTGTGCGCTCGCTCGCAAAATAAAATGGAATGATGAAACTGGGTAATTCCTCGAAGAGAAAAATGCAATAAGTACAATTGCGCAACAAAAGTAAGATCTCGGTCATAAATCAAGAAATAAACCAAAAATCGTAATCGAAAGATAAAAATCTGTAATTGTTTTCCCCTGTTTAGTTGCTAAAAATTGGTTACGTT  >  NC_000913/3730778‑3731073

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