Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	51.1%	68.0 / 142.5	90	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (38/6);  new base A (0/46);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 2.00e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-16

AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5009‑5307                                                                                                                                                       |                                                                                                                                                     aaTCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgag                                                                                                                                                      <  2:824391/151‑1 (MQ=255)       caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                                >  1:243953/1‑151 (MQ=255)       caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                                >  1:1029353/1‑151 (MQ=255)        aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                                >  1:1242406/1‑150 (MQ=255)           aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                             <  2:1154533/150‑1 (MQ=255)           aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                            <  2:854120/151‑1 (MQ=255)            aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                           <  1:261367/151‑1 (MQ=255)            aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                           <  1:391102/151‑1 (MQ=255)            aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                           <  2:428891/151‑1 (MQ=255)             aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                          <  1:393189/151‑1 (MQ=255)                cTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGt                                                                                                                                       >  1:543355/1‑151 (MQ=255)                     cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                   >  1:155645/1‑150 (MQ=255)                           ctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                            >  2:688877/1‑151 (MQ=255)                               aTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                        >  1:841245/1‑151 (MQ=255)                                 cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                      >  1:329687/1‑151 (MQ=255)                                       ggttttttttttAAGCACCCCGGTTTTTTTATAAAAAAATGATGAAAAAAAATAACGGGAAAAAAGGAAAATTTCTAAAAAAATGCGTTAATTTAGAGTTTAGGTTTGGGGAAAATAACACCCGCCGTTCTCATCGAGTAAtct                                                                                                                       <  1:228668/136‑1 (MQ=255)                                       gggTTTGCTTTTATGCAGCCCGTTTTTTTTATAAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                                <  1:212952/151‑1 (MQ=255)                                          tttGTTTTTATGCGGCCCGGTTTTTTTATAAAAAAATTAGGAAAAAAAAGAAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAATTTAGAAATGGGGATTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                              <  1:686205/150‑1 (MQ=255)                                          tttGTTTTTATGCAGCCGGTTTTTTTTATAAAAAAATTAGGAAAAAAAATGAAGGGAAAAAAGAAAAAATTCTAAAAAAAGGGGGTAACTAAGATATTAGGTTGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                               <  2:926996/149‑1 (MQ=255)                                          tttGTTTTTAGGCAGCCCGGTTTTTTTAGGAAAAAATTAGGGAAAAAAATGACGGGAAAAAAGGAAAAATTCTAAATAAATGGGTAAACTAAGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                               <  2:799558/149‑1 (MQ=255)                                          tttGTTTTTAGGCAGCCCGGTTTTTTAATAAAAAAATTAGGAAAAAAAATGAAAGGAAAAAAGAAAAAATTCTCAAAAAATGGGTTACTTTAGAGATTAGGAGTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                              <  1:2254/150‑1 (MQ=255)                                            tgtttttAGGCAGCCGGGTTTTTTTAGAAAAAAATTAGGGAAAAAAAGGAAAGGAAAAAGGGAAAAATTCTAAAAAAAGGCGGTACCTTAGAGATAGGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGa                                                                                                            <  2:568550/147‑1 (MQ=255)                                             gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                          >  2:709180/1‑151 (MQ=255)                                                   aggCACCCCGGTTTTTTTATAAAGAAATTAGGGAAAAAAATGACAGGAAAAAAGAAGAAATTCTAAAAAAAGGCGAAAACTAAGAGATTGGATTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      <  1:1201049/147‑1 (MQ=255)                                                      cggCCGGGTTTTTTTATGAAAAAATTAGGAAAAAAAAGGACAGGAAAAAAGAAAAAATTCTAAATAAGTGTGGAATCTTAGAGATTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      <  2:893585/144‑1 (MQ=255)                                                      cggCCGGGTTTTTTAATAAAGAAATTAGGGAAAAAAAAGAAAGGGAAAAAGAAAAAATACTAAAAAAACGGGGAACCTTAGAGATTAGGATGGGGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                     <  2:97802/145‑1 (MQ=255)                                                         ccctgtttttttAAAAAAAAATTAGGAAAAAAAAAGAAAGAAAAAAAGAAAAAATTCAAAAAAATTGCGGTAACTAAGAGATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                              <  1:275245/145‑1 (MQ=255)                                                         cccgttttttttATAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAATTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                    <  2:351703/139‑1 (MQ=255)                                                         cccgttttttttAGAAAAAATTTAGGGAAAAAAATGACAGGAAAAAAGGAAAATTTCACAATAAATGCGGTAATTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                 <  2:593525/142‑1 (MQ=255)                                                         cccgttttttttAAAAAAAAATTAGGGAAAAAAATAAAAGGAAAAAAGGAAAAATTCTAAAAAAAGGGGGTAATTAAGAGATTGGGTTTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                   <  2:566655/140‑1 (MQ=255)                                                         cccGGTTTTTTTATAAAAAAATAAGGAAAAAAAATGAAAGGAAAAAAGAAAAAATTCAAAAAAAAGGCGATAATTAAGAGATGAGGTTGGCGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                <  1:209460/149‑1 (MQ=255)                                                         cccGGTTTTTTTAAAAAAAAATTAGGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCTCAAAAAATGCGTTAACTAAGAGTTTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                   <  1:1007715/146‑1 (MQ=255)                                                          ccgttttttttAAAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGAAAAAATTCTAAAAAAGGGGGGTAACTGAGAAATTAGGATGCCGGAAATAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      <  1:874458/137‑1 (MQ=255)                                                          ccGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                             >  1:488819/1‑151 (MQ=255)                                                           cGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                            >  2:21775/1‑151 (MQ=255)                                                               tttttaATAAAGAAATTAGGAAAAAAAATGAAAGGAAAAAAGGAAAATTTCTAAAAAAAGGCGGAAACTTAAAAATTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          <  2:385092/149‑1 (MQ=255)                                                                tttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                        >  1:67154/1‑150 (MQ=255)                                                                ttttaAAGAAAAAATTATGAAAAAAAAGAACGGGAAAAAAGGAAAAATTCTCAATAAATGCGTTAACTAAGAAATTAGAATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                        <  1:1019029/150‑1 (MQ=255)                                                                      tGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                 >  2:953045/1‑151 (MQ=255)                                                                         agaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                                                                              >  2:1206463/1‑151 (MQ=255)                                                                           aaaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                            >  1:876899/1‑151 (MQ=255)                                                                            aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           >  2:433147/1‑151 (MQ=255)                                                                            aaTTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           <  2:349698/151‑1 (MQ=255)                                                                               taggaaaaaaagTGACAGGAAAAAAGAAAAAATACTAAAAAAAGGAGTTAACTAAGAGATGAGGAGGGCGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           <  1:1006617/141‑1 (MQ=255)                                                                                aTGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                        >  1:545699/1‑150 (MQ=255)                                                                                 tGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                       >  2:830113/1‑150 (MQ=255)                                                                                    aaaaaaagaaaaGGAAAAAAGAAAAAATTCAAAAAAAAGGCGGTAATTTAAAAATTGGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                          <  1:1117666/133‑1 (MQ=255)                                                                                    aaaaaaaTGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAACTAAGAGATTAGTTTGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                                   <  2:201645/149‑1 (MQ=255)                                                                                      aaaaaTGAGAGGAAAAAAGGAAAAATTCTCAAAAAATGGGGTAATTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                 <  2:210188/151‑1 (MQ=255)                                                                                       aaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                >  1:270630/1‑151 (MQ=255)                                                                                       aaaaTGAAAGGAAAAAAGGAAAATTTCTAAAAAAATGGGATAACTTAGATATTAGGATGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                <  1:855894/151‑1 (MQ=255)                                                                                         aaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               >  2:124313/1‑150 (MQ=255)                                                                                              cAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCAGATATCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                          >  1:466443/1‑150 (MQ=255)                                                                                               aGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                        >  2:1104508/1‑151 (MQ=255)                                                                                                gggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                       >  2:34714/1‑151 (MQ=255)                                                                                                ggaaaaaaGGAAAAATTCTCAATAAAGGGGGTAACTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                         <  2:624459/146‑1 (MQ=255)                                                                                                ggaaaaaaGGAAAAATTCTAAAAAAATGCGTTACTTAAGAGATTGGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                        <  2:905691/147‑1 (MQ=255)                                                                                                ggaaaaaaGAAAAAATTCTCAAAAAGTGCGAAAATTAAGAGATTGGGTTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                        <  1:632652/147‑1 (MQ=255)                                                                                                gaaaaaaaGAAAAAATCCAAAAAAAATGCGGAAACTAAGAGATTGGGAGTGGGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                            <  1:532847/143‑1 (MQ=255)                                                                                                   aaaagGGAAAATTTCTAAAAAAAGGCGGTAATTAAGAAATTAGGTTGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTa                                                    <  1:682647/151‑1 (MQ=255)                                                                                                   aaaaaGGAGAATTTCTAAATAAATGCGGTAACTTAGAGATTGGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                       <  2:370592/148‑1 (MQ=255)                                                                                                   aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                     >  1:1190608/1‑150 (MQ=255)                                                                                                   aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                     >  1:243387/1‑150 (MQ=255)                                                                                                   aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                     >  1:102248/1‑150 (MQ=255)                                                                                                   aaaaaGGAAAAATTCAAAAAAAATGGGGTAATTTAGAGATGAGGATGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                             <  1:540082/142‑1 (MQ=255)                                                                                                   aaaaaGAAAAAATTCTAAAAAAATGGGTTAATTTAGAGTTAAGGTTGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTTTGAAAAAGATGCAATCTAt                                                        <  2:819060/147‑1 (MQ=255)                                                                                                     aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                  >  1:581584/1‑151 (MQ=255)                                                                                                     aaaGGAAAAATTCTCAATAATGGCGGTACCTTAGAGATTGGGAGCGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                  <  2:719184/151‑1 (MQ=255)                                                                                                      aaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                 >  2:695523/1‑151 (MQ=255)                                                                                                      aaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                 >  1:952907/1‑151 (MQ=255)                                                                                                          aaaaatccaaaaaaaaTGCGTAAACTAAGAGATGAGGATGGCGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                             <  1:191013/138‑1 (MQ=255)                                                                                                           gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                            >  1:1097171/1‑151 (MQ=255)                                                                                                            aatttcaaaaTAAAGGCGATAATTTAGATATTAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                           <  2:452157/143‑1 (MQ=255)                                                                                                            aaaTTCTAAAAAAATGCGGTAACTTAGAGATTGGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAGAAAGATGCAATCTATCGTACTCGCACtt                                           <  2:641472/151‑1 (MQ=255)                                                                                                               ttCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                        >  1:120307/1‑151 (MQ=255)                                                                                                                    aaaaaaTGAGATAAGTAAGAGATTGGAATCGCGAAAAAAACCAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                    <  2:354651/147‑1 (MQ=255)                                                                                                                    aaaaaaGGCGGTAATTAAGAGATTAGGATTGGGAAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                       <  1:1207004/144‑1 (MQ=255)                                                                                                                    aaTAAATGCGGTAATTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                   <  2:261148/151‑1 (MQ=255)                                                                                                                        aaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               >  1:762864/1‑151 (MQ=255)                                                                                                                          tGCGGTAACTTAGAAATAAGGTTTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                              <  2:272159/150‑1 (MQ=255)                                                                                                                           gCGGTAAATTAGAGATTGGGATTGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                             <  1:514008/150‑1 (MQ=255)                                                                                                                            gggTAACTAAGAAATAGGGATGGGGGAAAAAAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAGGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               <  2:1098812/146‑1 (MQ=255)                                                                                                                             gttAACTTAGAGATTGGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                           <  2:35394/148‑1 (MQ=255)                                                                                                                                    tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                   >  1:649730/1‑151 (MQ=255)                                                                                                                                    tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                   >  2:773292/1‑151 (MQ=255)                                                                                                                                     agagATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg                  >  2:31281/1‑151 (MQ=255)                                                                                                                                      gagaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  2:630827/1‑151 (MQ=255)                                                                                                                                          ttAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc              >  2:573091/1‑150 (MQ=255)                                                                                                                                             ggATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca             <  2:969316/148‑1 (MQ=255)                                                                                                                                                 tGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg      >  2:1043121/1‑151 (MQ=255)                                                                                                                                                  ggggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc        <  2:1242431/146‑1 (MQ=255)                                                                                                                                                      agaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  >  2:455426/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                     AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5009‑5307

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