Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	42.3%	83.0 / 50.3	78	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (42/3);  new base A (0/33);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.55e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5011‑5307                                                                                                                                                         |                                                                                                                                                     tcatcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGgaga                                                                                                                                                     >  1:903542/1‑150 (MQ=255)  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGaa                                                                                                                                                    <  2:484498/150‑1 (MQ=255)   atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGACA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                   >  2:1262769/1‑150 (MQ=255)          aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                           >  1:1219139/1‑151 (MQ=255)            aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                         >  1:291069/1‑151 (MQ=255)             tCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAAccgccg                                                                                                                                        >  1:960356/1‑151 (MQ=255)                aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                         <  1:575307/147‑1 (MQ=255)                   cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                   >  2:704929/1‑150 (MQ=255)                            aaTCAGGCGGTGTTTGTTTTTAGGCAGCCTGTTTTTTATAAGAAAAAATTAGGAAA‑AAAAATGAGAGGAAAAAAGGAAAAA‑TTGTA‑AAAAATTGGGG‑TAACTTAGAGATTAGGATTGGTGAGAATAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                         <  2:588363/151‑1 (MQ=255)                                aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                     >  2:1295735/1‑151 (MQ=255)                                 gcggggttttttttttATGCGGCCCGGTTTTTTAAAAAAAAATTTAGGAAA‑AAAAAGGAAAGGAAAAAAGGAAAAT‑TCCTA‑AAAAAGGGCGG‑TAACTTGGAATTTGGGTTGGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                    <  1:1137105/139‑1 (MQ=255)                                  ggcGGGTTTTTTTTTTGGCCCCCCCGGTTTTTTAATAAAAAATTTAGGAA‑AAAAAAGGAAGGAAAAAAGGGAAAA‑TTCCCAAAAAAAGGCTG‑TAACTAAGAGATTGGAATTGGGAAAAATAAAAACCGCCGTTCTCATCGAGTAAtct                                                                                                                       <  2:1239184/146‑1 (MQ=255)                                     ggttttttttttATGCAGCCCGTTTTTTTAAGAAAAAAATTAGGAAA‑AAAAAGAAGAGGAAAAAAGGAAAAT‑TTCAA‑AAAAAGTGCGG‑TAACTTAGAGATAAGGTTGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   <  2:107455/140‑1 (MQ=255)                                       gTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                              >  2:887434/1‑151 (MQ=255)                                        tttttttttATGCGGCCGGTTTTTTTTAAAAAAAAATTAGGAAA‑AAAAATAAAAGGAAAAAAGAAAAAA‑TTCTA‑AATAAATGCGT‑TAACTTAGAGATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                              <  2:596365/145‑1 (MQ=255)                                           gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                          >  2:595313/1‑151 (MQ=255)                                             ttttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                         >  2:904843/1‑150 (MQ=255)                                                  tGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                   >  2:865480/1‑151 (MQ=255)                                                  tGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                   >  1:728083/1‑151 (MQ=255)                                                    cAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                 >  1:577665/1‑151 (MQ=255)                                                    cAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                 >  2:539127/1‑151 (MQ=255)                                                       cctgttttttttATAAAAAAATTATGAAA‑AAAATGGAAAGGAAAAAAGAAAAAT‑TTCAA‑AAAAAGTGCGT‑TACCTTAGAAATTAGGTTGGAGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                     <  2:597502/138‑1 (MQ=255)                                                       cccgtttttttaATAAAAAAATTAGGGAA‑AAAAAGGACAGAAAAAAAGAAAAAT‑TTCTA‑AATAAGTGGGT‑TATTTTAAAAATAAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                <  1:1229087/143‑1 (MQ=255)                                                         cagtttttttAAAAAAAAAATAAGGAAAAAAAATGAAAGGAAAAAAGGAAAAA‑TTCTA‑AAAAAAGGCGT‑TATCTAAGAAATAAGGATGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                             <  2:231500/147‑1 (MQ=255)                                                            tttttttAAAAAAAAATTAGGAAA‑AAAAATGACAGGAAAAAAGGAAAAA‑TTCTC‑AAAAAGTGCGA‑TAATTTAGAGATTAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                 <  2:1111941/142‑1 (MQ=255)                                                            cTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGAGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGCTCTCATCGAGAAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          >  2:365066/1‑150 (MQ=255)                                                             ttttttATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                         >  2:538204/1‑150 (MQ=255)                                                                   aTGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                    >  2:870786/1‑149 (MQ=255)                                                                    tGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                                  >  2:314040/1‑150 (MQ=255)                                                                     aaagaaATTAAGGAA‑AAAAATAAGAGGAAAAAAGGAAAAA‑TTCTC‑AAAAAAGGCGG‑TAACTGAGAGTTTAGTATCGGGAAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                 <  1:940248/149‑1 (MQ=255)                                                                     aaaaaaaTTAGGGAA‑AAAAAGGACAGGAAAAAAGAAAAAA‑TTCTA‑AAAAAAGGCGT‑TAACTTAGAGTTTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                 <  1:332186/146‑1 (MQ=255)                                                                      aagaaATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                >  1:1170137/1‑150 (MQ=255)                                                                      aagaaATTAGGAAA‑AAAAATAAAAGGAAAAAGGAAGAAA‑TTCTA‑AATAAATGCGG‑TATTTTAGAAATTAGAATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                                  <  1:549050/148‑1 (MQ=255)                                                                        gaaATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                             >  1:844244/1‑151 (MQ=255)                                                                         aaaTTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGTCAATGATAAAAGGAGTAAcc                                                                             >  1:982797/1‑150 (MQ=255)                                                                          aataagggAA‑AAAAAGGAAGGGAAAAAAGGAAAAT‑TTCTA‑AAAAAAGGGGG‑AAACTAAGAGTTTAGGATTGGGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           <  1:910474/145‑1 (MQ=255)                                                                             taatgaaaaaaaaGGACAGGAAAAAAGAAAAAA‑TTCAA‑AAAAAAGGCGG‑TAACTAGGAATTTGGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                          <  2:1067/142‑1 (MQ=255)                                                                              aTGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                        >  2:407449/1‑150 (MQ=255)                                                                                   aaaaaaaTGAAAGGAAAAAAGGAGAAA‑TTCTA‑AAAAAATGCGG‑TAACTAAAAGATAAGGTTGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         <  2:367608/143‑1 (MQ=255)                                                                                   aaaaaaaGGACA‑GGAAAAAGGAAAAACTTCAA‑ATTAACGGGGT‑TACCTTGGAGATTGGGAGAGAGAAAAAAAAAATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                       <  2:1137595/145‑1 (MQ=255)                                                                                     aaaaaTGAGAGAAAAAAAGGAAAAA‑TTCTA‑AAAAAAGGCGG‑TAACTTAGAAATTGGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                        <  1:578889/144‑1 (MQ=255)                                                                                      aaaatgaaggaaaaaaaGGAAAAA‑TTCAA‑AAAAAAGGCGT‑TAATTAAGAGATTGGGATTGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                <  1:200911/139‑1 (MQ=255)                                                                                      aaaagaaaaGGAAAAAAGAAAAAT‑TTCAC‑AAAAAAGGCGT‑TACTTTAGAGTTTAGGATGGGAAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                <  2:731674/143‑1 (MQ=255)                                                                                      aaaaagaaagaaaaaaaGGAAAAT‑TTCAA‑AAAAAATGGGT‑TAACTAAGAAATAGGGATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                <  2:1147219/139‑1 (MQ=255)                                                                                        aaggaaaGGAAAAAAGAAAAAA‑TTCTA‑AATAAGTGGGT‑TAACTTAGAAATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                              <  2:460882/145‑1 (MQ=255)                                                                                             cAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATCAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                          >  2:202780/1‑150 (MQ=255)                                                                                                 aaaaaaGGA‑AAA‑TTTCA‑AAAAAATGGGGTTAATTAAGAGTTTAGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                      <  2:1267605/149‑1 (MQ=255)                                                                                                 aaaaaaGAAAAAA‑TTCTA‑AATAAATGCGG‑TATCTGAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                      <  1:817627/149‑1 (MQ=255)                                                                                                  aaaaaGAAAAAA‑TACAA‑AAAAAAGGGGG‑TAACTAAGAGATGAGGAGTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                        <  2:77679/147‑1 (MQ=255)                                                                                                    aaaGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                  >  2:601110/1‑151 (MQ=255)                                                                                                      aGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                                >  2:56007/1‑151 (MQ=255)                                                                                                         agaAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                             >  1:1122832/1‑151 (MQ=255)                                                                                                         aaaaa‑TTCTC‑AAAAATTGCGG‑TAACTAAGAGATTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                                <  1:426107/146‑1 (MQ=255)                                                                                                          gaAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                            >  1:60029/1‑151 (MQ=255)                                                                                                            aa‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                          >  2:729273/1‑151 (MQ=255)                                                                                                                tctc‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                       >  2:415691/1‑151 (MQ=255)                                                                                                                    aaaaaaaGGCGG‑AAACTTAGAGATAAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                     <  2:31524/146‑1 (MQ=255)                                                                                                                     aaaaaaGGCGT‑TAATTAAAAGATTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                        <  1:539143/143‑1 (MQ=255)                                                                                                                     aaTAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                   >  2:367688/1‑151 (MQ=255)                                                                                                                     aaTAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                   >  2:1304815/1‑151 (MQ=255)                                                                                                                       tAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                 >  2:426815/1‑151 (MQ=255)                                                                                                                       tAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                 >  1:358504/1‑151 (MQ=255)                                                                                                                        aagtGCGT‑TAACTTAGAAATTAGTAGTGAAGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                  <  1:911447/146‑1 (MQ=255)                                                                                                                         aaTGCGG‑TAACTGGGATATTAGGTTTGCGGAAAAAAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                    <  1:218728/146‑1 (MQ=255)                                                                                                                          aTGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                              >  2:1133389/1‑151 (MQ=255)                                                                                                                           tGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                             >  2:446956/1‑151 (MQ=255)                                                                                                                              gg‑TATCTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                           <  1:87123/150‑1 (MQ=255)                                                                                                                                 taccTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCAGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTc                        >  2:1214568/4‑151 (MQ=255)                                                                                                                                    cTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                      >  2:790518/1‑150 (MQ=255)                                                                                                                                    cTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                      >  2:292697/1‑150 (MQ=255)                                                                                                                                        gagaTTAGGATTGCGTAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg                  >  1:224495/1‑150 (MQ=255)                                                                                                                                        gagaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  2:159821/1‑151 (MQ=255)                                                                                                                                         agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  2:785992/1‑150 (MQ=255)                                                                                                                                          gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  1:496704/1‑149 (MQ=255)                                                                                                                                            ttgggATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    <  2:318660/141‑1 (MQ=255)                                                                                                                                              aggtgggCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 <  1:1165716/139‑1 (MQ=255)                                                                                                                                               ggATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    <  2:1096246/141‑1 (MQ=255)                                                                                                                                                 atgggggAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg         <  1:684216/145‑1 (MQ=255)                                                                                                                                                     cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg    >  1:724314/1‑151 (MQ=255)                                                                                                                                                        agaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  >  1:163477/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                     TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAAGGAGAAA‑TTCTC‑AATAAATGCGG‑TAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5011‑5307

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