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

ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5013‑5307                                                                                                                                                       |                                                                                                                                                      atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAACTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAAt                                                                                                                                                    >  1:807068/1‑150 (MQ=255) atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGATTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATa                                                                                                                                                   <  1:633019/151‑1 (MQ=255)   caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATaac                                                                                                                                                 <  2:498038/151‑1 (MQ=255)        aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AAccg                                                                                                                                           <  1:738304/151‑1 (MQ=255)            cTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGt                                                                                                                                       >  2:460271/1‑151 (MQ=255)             tATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGtt                                                                                                                                      >  2:740250/1‑151 (MQ=255)              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AAccgcc                                                                                                                                         <  2:71117/147‑1 (MQ=255)              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTtc                                                                                                                                     >  2:720240/1‑151 (MQ=255)                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTtctc                                                                                                                                   >  2:116616/1‑150 (MQ=255)                   ttATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATc                                                                                                                                >  2:363899/1‑151 (MQ=255)                    tATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCg                                                                                                                               >  1:662453/1‑151 (MQ=255)                            tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GCGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAAtc                                                                                                                        >  2:481225/1‑150 (MQ=255)                               ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                    >  1:676890/1‑151 (MQ=255)                                gCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   >  1:863138/1‑151 (MQ=255)                                 cgggttttGTTTTTATGCACCCCGTTTTTTTTATAAAGAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAATTTCTAAATAAAGGGGGTAACTA‑A‑GAGTTTA‑GGATT‑GGGGAAAATAAC‑ACCCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   <  2:705280/145‑1 (MQ=255)                                  cGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                 >  1:818037/1‑151 (MQ=255)                                   ggttttttttttATGCGCCCCGTTTTTTTTATAAAAAAATTATGAAAAAAAAAGAGAGGAAAAAAGGAAAATTTCTAAAAAAATGGGATAATTT‑A‑GAGTTTA‑GATTT‑GGGAAAAAAAAA‑AACCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                                <  1:966237/143‑1 (MQ=255)                                   ggttttGTTTTTATGCAGCCTGTCTTTTTAATGAAAAAATTAGGAAAAAAAAAGAGAGGAAAAAGGGAAAACTTCTCAAAAATGGCGTTAACAT‑A‑GAGATAA‑GTATG‑GAGGAAAACAAA‑AGCCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                                <  1:575054/148‑1 (MQ=255)                                   ggttttGTTTTTAGGCAGCCCGGCTTTTTTATGAAGAAATTATAGAAAAAAATGACAGGAAAAAAGAAGAAATTCTAAATAAAGGCGGTAACTT‑A‑GAAATAG‑GGATT‑GGGGAGAAAAAA‑AACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   <  1:581712/145‑1 (MQ=255)                                   gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGTAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                 >  2:580892/1‑150 (MQ=255)                                      tttGTTTTTACGCAGCCGGGTTTTTTTAAAAAAAAATTAGGAAAAAAAAAGACAGGAAAAAAGGAGAAATTCTAAAAAAGTGGGTTAACTA‑A‑GAGATTA‑GGATT‑GCGAAAAATAAA‑AACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                               <  2:632726/149‑1 (MQ=255)                                      tttGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                             >  2:922155/1‑151 (MQ=255)                                           ttttATGCAGCCCGGTTTTTTTATGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAATTTCTAAATAAAGGCGGTAATTT‑A‑GAGATTA‑GGTTG‑GCGGAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                          <  2:557129/149‑1 (MQ=255)                                             ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                       >  2:865206/1‑150 (MQ=255)                                              tATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      >  2:228535/1‑150 (MQ=255)                                               aTGCAGCCCGGTTTTTTTAAAAAAAAATAAGGGAAAAAAATGACAGGAAAAAAGAAAAAATTCTAAAAAAATGCGTTAACTA‑A‑GATATTG‑GGATG‑GCGAAAAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                     <  1:912732/150‑1 (MQ=255)                                                    gCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                >  1:248113/1‑150 (MQ=255)                                                     cccGGTTTTTTTATGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAGAAATTCTCAAAAAAGGCGGTAACTT‑A‑GAGTTTA‑GGATT‑GCGGAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                                                                                                  <  1:259015/147‑1 (MQ=255)                                                      ccgcttttttaAGGAAAAATTTATGGAAAAAAATGACAGGAAAAAAGAAAAAATTCTAAAAAAATG‑GGGAACTTAA‑GAGA‑TG‑GGATTGGGGGAAAAAAAA‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                             <  2:444302/146‑1 (MQ=255)                                                          tttttttATAAAGAAATTAAAGAAAAAAATGATAGAGAAAAGGGA‑AAATTTCCAAAAAAAGGGGCAAATT‑AGGCGTTTA‑GTTTG‑GCGGAAAAAAAG‑ACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          <  2:525772/149‑1 (MQ=255)                                                           ttttttAAAAAAAAATTAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAATTT‑A‑GAAATTA‑GGATT‑GCGGAAAAAAAC‑ACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                         <  1:720113/150‑1 (MQ=255)                                                            tttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                        >  2:808022/1‑150 (MQ=255)                                                                       aaaTTAGGGAGAAAAAGGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTG‑A‑GAGATTA‑GGATT‑GCGGAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                             <  2:623064/150‑1 (MQ=255)                                                                       aaaTAATGGAAAAAAATGAAAGGAAAAAGGGAAAAATCCTAAAAAAATGCGGTAATTT‑A‑GAGATTA‑GGATG‑GCGAAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                             <  2:738669/150‑1 (MQ=255)                                                                           tATGGAAAAAAACAACAGAAAAAAAGGAGAAACCCCAAACAGAGGCAAAAACTG‑A‑GAGATGA‑GATAT‑GCGGAAAAAACC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         <  2:940253/150‑1 (MQ=255)                                                                                  aaaatagaaaggaaaaaaggaaaaaTACTAAATAAATGGGGTAACTT‑A‑‑AAATTAGGGATA‑GCGAAAAAAAAA‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                 <  1:453200/142‑1 (MQ=255)                                                                                    aaaaaagaGGAAAAAATGAAAAATTCAAAAAAAAGGGAGTAGTTT‑A‑GAGTTAA‑GGATT‑GGGAAAAAAAAC‑ATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               <  1:942169/144‑1 (MQ=255)                                                                                      ataagggaaaaaaaGGAAAAATTCTAAAAAACGGCGTAATTTT‑A‑GAAATTA‑GAATG‑GCGAAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                             <  1:759149/142‑1 (MQ=255)                                                                                            ggaaaaaaGGAAAAATTCTCAAAAAATGGGGTATCTT‑A‑GAGATTA‑GGATT‑GCGGAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                       <  2:450265/148‑1 (MQ=255)                                                                                               aaaaaGAAGAAATTCTAAATAAACGGGTTACCTA‑A‑GATATCG‑GCATC‑G‑GGAAAAAAACAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTa                                                    <  2:446455/151‑1 (MQ=255)                                                                                                  aaGGAAAAATTCTAAATAAAGGCGGTAACTT‑A‑GAGATAA‑GGAGT‑GCGGAAAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                 <  1:738754/151‑1 (MQ=255)                                                                                                      agaAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGCTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                             >  1:688263/1‑151 (MQ=255)                                                                                                         aaTTCT‑AAAAAATGGCGTTACTTAA‑GAGATTA‑GGATG‑GCGAAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                          <  1:137558/151‑1 (MQ=255)                                                                                                                aataaggggggTAACTG‑A‑GAGATTA‑GTATT‑GCGGAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                      <  2:865632/139‑1 (MQ=255)                                                                                                                 aTAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                  >  1:35417/1‑151 (MQ=255)                                                                                                                  tAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                 >  1:435939/1‑151 (MQ=255)                                                                                                                    aaaggcggaaCTTAA‑GAGTTTA‑GGTTG‑GCGAAAAATAAA‑ACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                <  2:404473/143‑1 (MQ=255)                                                                                                                    aaTGCGGTAACTT‑A‑GAGATTA‑GGATT‑ACGGAGAATAAC‑TACCGCCGTTCTCATCGAGTAATCTAAGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAACGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                >  1:278474/1‑150 (MQ=255)                                                                                                                     aTGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                              >  2:258459/1‑151 (MQ=255)                                                                                                                     aTGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                              >  1:231707/1‑151 (MQ=255)                                                                                                                       gCGGTAATTA‑A‑GAGATTA‑GGATG‑GGGGAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                            <  2:653948/151‑1 (MQ=255)                                                                                                                          gTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCt                         >  2:933266/1‑151 (MQ=255)                                                                                                                           gataagtaaagaaTTA‑GGATT‑GCGGAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                            <  1:921987/137‑1 (MQ=255)                                                                                                                                     agaTTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCGCTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  1:27843/1‑150 (MQ=255)                                                                                                                                      gaTTA‑GGAGT‑GCGGAAAAAACC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgca                <  1:734068/150‑1 (MQ=255)                                                                                                                                       aTTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc              >  2:865839/1‑151 (MQ=255)                                                                                                                                          a‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac            >  2:297884/1‑150 (MQ=255)                                                                                                                                            ggATT‑GCGAAAAAAAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca             <  1:198953/148‑1 (MQ=255)                                                                                                                                            aggttgggggAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTCTCCCTGGTTCTGGTCGCTCCCa                     <  2:863295/133‑1 (MQ=255)                                                                                                                                             gATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg         >  2:703311/1‑151 (MQ=255)                                                                                                                                                t‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg      >  1:770422/1‑151 (MQ=255)                                                                                                                                                    ggAAAAAAAA‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt       <  2:630330/147‑1 (MQ=255)                                                                                                                                                      agaATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  >  2:472706/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTT‑A‑GAGATTA‑GGATT‑GCGGAGAATAAC‑AACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5013‑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: ‑