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

TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGC  >  NC_000913/5011‑5304                                                                                                                                                        |                                                                                                                                                  tcatcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGgaga                                                                                                                                                  >  1:900658/1‑150 (MQ=255)     caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                             >  1:167712/1‑151 (MQ=255)     caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                             <  2:524967/151‑1 (MQ=255)        tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                          >  1:350805/1‑151 (MQ=255)           aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                       >  1:941386/1‑151 (MQ=255)               tATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGtt                                                                                                                                   >  2:1011362/1‑151 (MQ=255)               tATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGtt                                                                                                                                   >  2:783537/1‑151 (MQ=255)                       aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                            >  1:679184/1‑150 (MQ=255)                         ctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                         >  2:591868/1‑151 (MQ=255)                            aaTCAGGCGGGGTTTGTTTTTATGCACCCCGTTTTTTTTATAAAAAAATT‑AGGGAAAAAAATGAAAGGAAAAAAGGAAAAA‑TTCTAAAAAAAGGCGGTAA‑CTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTaa                                                                                                                       <  2:876321/150‑1 (MQ=255)                             accAGGCGGGTTTTTCTTTTAGGAGGCCGGGTTTTTTTATGAAAAAATT‑ATGAAAAAAAATGCAAGGAAAAAAGGAAAAA‑TTCAAAAAAAACGGGGTAC‑CTAAGAGTTTGGGTTGGGGGAGAAAAACAACCGCCGTTCTCATCGa                                                                                                                           <  1:585277/143‑1 (MQ=255)                              tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                    >  1:577326/1‑151 (MQ=255)                                   cggggTTTGTTTTTAGGCAGCCCGGTTTTTTTAGAAAGAAATT‑ATGGAGAAAAATGACAGGAAAAAAGGAGAAA‑TTCTAAATAAAGGCGGTAA‑TTTAGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                <  1:273850/148‑1 (MQ=255)                                     ggttttctttttAGACACCCCGGTTTTTTTGGAAAGAAATT‑AAGAAAAAAAAGGAAAGGAAAAAAGAAAAAA‑TTCAAATAAAATGGTGTAC‑TTTAGAGTTTGGAGTTGAGGAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                             <  1:166605/143‑1 (MQ=255)                                     gggTTTGCTTTTATGCAGCCCGTTTTTTTTATAAAAAAATT‑AGGAAAAAAAATGACAGGAAAAAAGGAAAAA‑TTCTAAATAAATGCGGTAA‑CTAAGAGATTAGGTTTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                <  2:653547/148‑1 (MQ=255)                                     gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                              >  2:313833/1‑150 (MQ=255)                                             ttttAAGAACCCCGGTTTTTTAATAAAAATATT‑AGGAAAAAAAAGGACAGAAAAAAAGAAAAAT‑TTCTCAATAAGTGGGCAAA‑CTAATAAATTTGTGTTGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                        <  1:281229/148‑1 (MQ=255)                                              tttATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                     >  1:751530/1‑150 (MQ=255)                                              tttAGGCAGCCGGGTTTTTTAATAAAAAAATT‑AGGGAAAAAAAGGAAAGGAAAAAGGGAAAAA‑TTCTCAAAAAATGCGGTAA‑TTTAGAGATTAGGTTGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                     <  2:819410/150‑1 (MQ=255)                                                tATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                   >  2:428419/1‑150 (MQ=255)                                                    cAGCCCGGTTTTTTTATAAAAAAATA‑AGGGAAAAAAAGGAAAGGAAAAAAGGAAAAA‑TTCTAAAAAAAGGCGGTAA‑TTTAGAGATTAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                                                                                               <  1:962399/150‑1 (MQ=255)                                                          gttttttttATAAAAAAATT‑ATGGAAAAAAATGACAGGGAAAAAGGAAAAA‑TTCACAATAAATGCGTTAA‑CTTAGAGATTAGGATTGGGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                         <  2:78167/147‑1 (MQ=255)                                                          gtttttttaATAAAAAATTT‑AGGGAAAAAAAGGAAAGGAAAAAAGGAAAAA‑TTCTAAAAAATTGCGTTAC‑CTTAGAAATAAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                         <  1:407990/147‑1 (MQ=255)                                                          ggCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                         >  2:810938/1‑150 (MQ=255)                                                            tttttttATGAAAAAATT‑AGGGAAAAAAATGACAGGAAAAAAGGAAAAA‑TTCTAAAAAAATGCGGTAA‑TTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                       <  2:560264/149‑1 (MQ=255)                                                             ttttttATAAAAAAATT‑AGGGAAAAAAAAGAAAGGAAAAAAGGAAAAA‑TTCTAAAAAAATGTGTTAA‑TTTAGAAATAAGGAGGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                      <  1:394092/150‑1 (MQ=255)                                                              tttttAAGAAGAAATA‑AGGGAAAAAAGGGAAAGAAAAAAAGGAAAAT‑TTCTAAAAAATGGCGTTAA‑CTAAGATATAGGGGGTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                     <  1:703025/150‑1 (MQ=255)                                                                tttaaagaaaaaTTAATGGAAAAAAATGACAGGAAAAAAGGAAAAA‑TTCCAAAAAAATCCGCTAT‑ATAAGAGTTTAGGAGCGGGGAAAATAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                     <  2:195617/140‑1 (MQ=255)                                                                   ataaaaaaaTT‑ATGGAAAAAAAGGACAGGAAAAAAGGAAAAA‑TTCTAAAAAAATGGGTTAA‑CTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                <  2:430345/144‑1 (MQ=255)                                                                        gaaATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                          >  2:895551/1‑151 (MQ=255)                                                                            tt‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                      >  2:728157/1‑151 (MQ=255)                                                                             taaggaagaAAAAGGACGGGGAAAAGGGAAAAT‑TTCAAAAAAAGTGCGTTAA‑CTAAGAGTTTGGGTTTGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                        <  2:753858/143‑1 (MQ=255)                                                                                   agaAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                  >  2:449308/1‑149 (MQ=255)                                                                                   agaAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                  >  2:530075/1‑149 (MQ=255)                                                                                    gaAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                  >  2:964568/1‑148 (MQ=255)                                                                                    aaaaaaTAACGGGAAAAAAGAAAAAA‑TTCTAAATAATGGCGTTAT‑CTAAGAAATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa                                                                   <  1:537700/146‑1 (MQ=255)                                                                                      aaaaAGACAGGAAAAAAGGAAAAA‑TTCTCAAAAAGTGCGTTAC‑CTAAGATAGAAGGATGGGGGAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                             <  1:1030364/151‑1 (MQ=255)                                                                                        aagaaCAGGAAAAAAGGAGAAT‑TTCTAAAAAAATGCGATAA‑CTTAGAAATTAGGGTTGCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                           <  2:1035099/147‑1 (MQ=255)                                                                                          tGAAAGGAAAAAAGGAAAAA‑TTCTAAAAAAATGCGGTAA‑TTAAGAATTTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                          <  2:751050/150‑1 (MQ=255)                                                                                              aGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                     >  2:469444/1‑151 (MQ=255)                                                                                               ggaaaaaaGGAAAAT‑TTCTAAAAAAGGGGGGAAA‑TTAAGAAATTGGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                    <  1:775481/148‑1 (MQ=255)                                                                                               ggaaaaaaGGAAAAA‑TTC‑CAAAAAAAGGGGAAACCTTAGAGATAAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                    <  2:905427/148‑1 (MQ=255)                                                                                                 gaaaaagaaaaaATTTCTAAAAAAATGCGGTAA‑CTAAGAAATTAGGAGTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                      <  1:1039346/148‑1 (MQ=255)                                                                                                  aaaaaGAAAAAA‑TTCACAATAAAGGCGATAA‑CTAAGAAATTAGGATTGCGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                    <  1:292383/148‑1 (MQ=255)                                                                                                     acggAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                              >  1:797218/3‑151 (MQ=255)                                                                                                      aGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                             >  1:188347/1‑151 (MQ=255)                                                                                                           aaa‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                        >  2:315133/1‑151 (MQ=255)                                                                                                           aaa‑TTCAAAATAATTGGGTTAA‑CTAAGAGATTGGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                        <  1:583892/151‑1 (MQ=255)                                                                                                             a‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTc                                      >  2:1030970/1‑151 (MQ=255)                                                                                                               ttCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                     >  2:883572/1‑151 (MQ=255)                                                                                                               ttCTCAAAAAAGGTGGAAA‑CTTAGATATTGGGATCGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTc                                      <  2:755695/150‑1 (MQ=255)                                                                                                                 ctaaaaaaaTGCGATAA‑CTAAGATATAGGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                     <  1:591746/143‑1 (MQ=255)                                                                                                                  tcAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                  >  1:532952/1‑151 (MQ=255)                                                                                                                    aaaaaaGGCGGTAA‑CTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                   <  2:661359/145‑1 (MQ=255)                                                                                                                      tAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                              >  2:77814/1‑151 (MQ=255)                                                                                                                      tAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                              >  2:823194/1‑151 (MQ=255)                                                                                                                       aaaTGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                             >  1:161394/1‑151 (MQ=255)                                                                                                                         aTGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                           >  1:65118/1‑151 (MQ=255)                                                                                                                          tGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAAGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                          >  1:232806/1‑151 (MQ=255)                                                                                                                                aa‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                    >  2:122132/1‑151 (MQ=255)                                                                                                                                       gagaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc              >  1:20147/1‑151 (MQ=255)                                                                                                                                       gaatttagagtggCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg               <  1:434587/138‑1 (MQ=255)                                                                                                                                         gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGCTAAAAGGAGTAACCTGTGAAAAAGCTGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc              >  1:73779/1‑149 (MQ=255)                                                                                                                                         gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc              >  1:327476/1‑149 (MQ=255)                                                                                                                                              ggATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg       <  1:889135/151‑1 (MQ=255)                                                                                                                                                   gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc  >  2:213964/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                  TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATT‑ATGGAGAAAAATGACAGGGAAAAAGGAGAAA‑TTCTCAATAAATGCGGTAA‑CTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGC  >  NC_000913/5011‑5304

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: ‑