Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,862	0	C	A	40.7%	97.8 / 28.9	59	S55R (AGC→AGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base C (22/13);  new base A (0/24);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 1.96e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATG  >  NC_000913/1755716‑1756012                                                                                                                                                   |                                                                                                                                                        caaacatgatctcatatcagggacttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCaaa                                                                                                                                                      <  2:650506‑M1/111‑1 (MQ=255)   aacatgatctcatatcagggacttgttcgcaccttccTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACt                                                                                                                                                   <  1:953459‑M1/114‑1 (MQ=255)           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)                                                                                                                                                   |                                                                                                                                                        GTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAA‑CCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATG  >  NC_000913/1755716‑1756012

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