Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,871	0	T	A	44.4%	68.3 / 45.5	63	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (26/9);  new base A (0/28);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 2.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGC  >  NC_000913/1755726‑1756017                                                                                                                                                  |                                                                                                                                                    ctcatatcagggacttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑c                                                                                                                                               <  2:518899‑M1/122‑1 (MQ=255)   catatcagggacttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTaaaa                                                                                                                                                 >  2:934824‑M1/28‑148 (MQ=255)   catatcagggacttgttcgcaccgtccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑cc                                                                                                                                              >  1:1000395‑M1/28‑150 (MQ=255)    atatcagggacttgttcgcgccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑ccgc                                                                                                                                            >  1:68307‑M1/27‑151 (MQ=255)      atcagggacttgttcgccccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCGCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑ccgcc                                                                                                                                           >  1:479123‑M1/25‑150 (MQ=255)        cagggacttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCt                                                                                                                                        <  1:198143‑M1/129‑1 (MQ=255)         agggacttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCt                                                                                                                                        <  2:634117‑M1/129‑1 (MQ=255)            gacttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATcc                                                                                                                                    <  1:908449‑M1/133‑1 (MQ=255)               ttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGc                                                                                                                                 <  2:837807‑M1/136‑1 (MQ=255)                  ttcgcaccttccttAGCAAAAATGCGGGGCGCTGGCATAAACGTTATGCGTCTGAACTTCTCTCAGGGTGAATATGCAGAACAGGGACAGCGCATCCAAATTCTGCGAAGGGTGATGAAAAAAACTGAAAAAA‑CCGCCGCTATCCTGCtt                                                                                                                               <  1:822033/138‑1 (MQ=255)                    cgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGAt                                                                                                                            <  2:974196‑M1/141‑1 (MQ=255)                     gcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGAt                                                                                                                            <  1:650521‑M1/141‑1 (MQ=255)                      ccccttccttAGCAAAAAGGCGGGACGCTGGCATAACCTTTATGCGTTTAACTTTCTCTCAGGGTGAATAGGAAAACCAGGGCCAGCCCTTTCAGAATCTGCGCAAGGTGATGAGCAAAACGGGAAAAA‑CCGCCGCTATCCTGCTTGATAc                                                                                                                          <  1:872470/143‑1 (MQ=255)                          ttccttAGCAAAAATGTGGGACGCTGACATGACCGTTATGCGTCTGACCTTCTCTCATGGTAACTAGGAAAAACACGGCCAGCGCATCCAAATTCTGCGAAAGGTGAGGAGAAAAAGGGGAAAAA‑CCGCCGCTATCCTGCTTGATAcc                                                                                                                         <  1:665055/144‑1 (MQ=255)                          ttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCaa                                                                                                                       <  2:1008931‑M1/146‑1 (MQ=255)                              ttAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                                  >  2:1025369/1‑151 (MQ=255)                                 gtaaaaaGGCGGGACGCTGGCATGAACTTTATGCGTCTAACCTTCTCTCATGTTGACTAGAAAAACCACGGTCGGCGCATTCAAATCCTGCGCAAGGTGAGGAGAAAACCTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCg                                                                                                                <  1:1004974/147‑1 (MQ=255)                                 gCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCgg                                                                                                               >  2:323062/1‑151 (MQ=255)                                    aaaaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGaa                                                                                                             >  1:944907/1‑150 (MQ=255)                                     aaaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCGACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAAt                                                                                                           >  1:682802/1‑151 (MQ=255)                                      aaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATc                                                                                                          >  1:772932/1‑151 (MQ=255)                                                  tGGCATGAACTTTATGCGTCTGACCTCCTCTCATGGTGACTATGCAGAACCCGGTCAGCCCTTTCAGAATCTCCCCAACGTGAGGAGAAAACCTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCAt                                                                                                  <  2:362267/147‑1 (MQ=255)                                                   ggaataactgTTTTGCGTCTGAACTTCTCCCATGGTAAATATGCAAAACAGGGTCAGCGCATTAAAATCCGGCGAAAGGTGATGAGAAAACCTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAAc                                                                                             <  1:801225/142‑1 (MQ=255)                                                   ggaataaaCTTTGTGCGTCTGATTTCCTCCCGTGGTAACTAGGAAGAACACGGCCGGCCCATTCAAAACCGCCGAAAGTTAATAAGAAAAAGGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGaa                                                                                               <  1:42277/143‑1 (MQ=255)                                                            gTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                    >  1:638045/1‑151 (MQ=255)                                                             tttggcttCTGACCTTCTCTCATGGTGACTAGGCAAACCAGGGCCAGCGCATTCAAAATCGGCGCAAGGGGAGGAGAAAAATTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                     <  2:496584/142‑1 (MQ=255)                                                              tATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                  >  1:250552/1‑151 (MQ=255)                                                               aggcttttaactttCTCTCATGGTGACTATGCAGAACACGGTCAGCGAATTCAGAACCGGCGCATCGTTATGAAAAAAAGGGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                  <  2:607718/138‑1 (MQ=255)                                                                        aaCTCCTCCCTTGGTGATTAGGAAAAACACGGCCAGCCCATTAAAATTCTGCAAAACTTAAGGAGAAAACTTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgacg                                                                           <  1:1032532/148‑1 (MQ=255)                                                                           ttctctaggggTGACTATGAAGACCGCGGTCAGCGAATCCAAAATCTGCGCAACGTGAGGAGAAAAACTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                            <  2:1018035/135‑1 (MQ=255)                                                                           ttCTCTCTGGTTGACTAACCAGAACACGGCCAGCTTATACATACTACGCGAAACGAAAAAAGAAAAACGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtc                                                                       <  1:447217/149‑1 (MQ=255)                                                                            tctctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                    >  2:317940/1‑151 (MQ=255)                                                                               ctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGa                                                                  >  1:536193/1‑150 (MQ=255)                                                                               ctcAGGGTGACTATGAAAAACGCGGCCGGCCCATTCAAAATCGGCGCACCGTGATGAGAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                     <  2:576409/147‑1 (MQ=255)                                                                                  aTGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                              >  1:180239/1‑151 (MQ=255)                                                                                    ggtacttcggCAAAACACGGCCAGATCATTAAGAACCCGCCCAACGTAATGA‑AAAAACCGGTAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                              <  1:969505/140‑1 (MQ=255)                                                                                    ggggACTATGAAAAACGCGGTCGGCGCATTAAAAATCGGCGCAACGTGATGAAAAAAATTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCt                                                             <  2:726578/147‑1 (MQ=255)                                                                                    ggTGACTATGCAAAACCCGGTCGGCCCATACAGAACCTGCGAAACGTAAGAAAAAAAACGGGAAAAC‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                     <  1:94222/142‑1 (MQ=255)                                                                                        aCTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTc                                                         >  2:986319/1‑150 (MQ=255)                                                                                          tATGCAAAACCCGGCCAGCGCATTCAAAATCTGCGCAAGGTGATGAGAAAAACTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCa                                                        <  2:679515/149‑1 (MQ=255)                                                                                           aTGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAc                                                     >  1:78597/1‑151 (MQ=255)                                                                                               aGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                                 >  1:795577/1‑151 (MQ=255)                                                                                                    acGGACAGCGCATTCAAAATCTGCGCAGGGTAAGGAAAAAAACGGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                   <  1:982718/144‑1 (MQ=255)                                                                                                      ggccAGCACATTCAGAATCGCCCCAACTTGATAAGAAAAACTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                           <  2:281783/147‑1 (MQ=255)                                                                                                          aGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                      >  1:921518/1‑151 (MQ=255)                                                                                                           gcgcACACAGAATTTGCCCAACAGTATGAACAAAACTGGAAAAC‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATTAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                                 <  1:62403/139‑1 (MQ=255)                                                                                                            ccctttCAAAACCGGCCCAAAGTGATGACAAAACCTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACt                                     <  2:839709/146‑1 (MQ=255)                                                                                                              cATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                  >  2:903326/1‑151 (MQ=255)                                                                                                                ttCAGAATCTCCGCAACGTGATGAGAAAAACGGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                  <  1:149773/149‑1 (MQ=255)                                                                                                                    aaccccgcgaAATGTGAGGAACAAAAGTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                            <  1:650393/145‑1 (MQ=255)                                                                                                                     aaTCTGCGCAACGGAATGAGAAAAAGTGTAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTg                           <  2:982139/151‑1 (MQ=255)                                                                                                                          gccccccGTGATGAGAAAAGCGGGAAAAT‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt                       <  2:278039/144‑1 (MQ=255)                                                                                                                           cgcAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGGTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCg                     >  1:1022597/1‑151 (MQ=255)                                                                                                                              aaCGTAAAGAGAAAAAGTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                   <  2:886160/150‑1 (MQ=255)                                                                                                                                gtagatgaGAAAAAGTGAAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc                <  1:217551/148‑1 (MQ=255)                                                                                                                                       aGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGaa          >  2:671284/1‑150 (MQ=255)                                                                                                                                       aGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGaa          >  1:9517/1‑150 (MQ=255)                                                                                                                                         cAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg       >  1:506988/1‑151 (MQ=255)                                                                                                                                         cAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg       >  1:556572/1‑151 (MQ=255)                                                                                                                                          aaaaCTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAAt        <  2:945043/149‑1 (MQ=255)                                                                                                                                          aaaaCTGGAAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg      <  2:215036/151‑1 (MQ=255)                                                                                                                                           aaaCTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGt     >  2:482477/1‑151 (MQ=255)                                                                                                                                              cTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGc  >  1:514497/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGC  >  NC_000913/1755726‑1756017

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