Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,871	0	T	A	41.5%	90.2 / 38.7	65	G58G (GGT→GGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (35/3);  new base A (0/27);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 2.80e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755721‑1756020                                                                                                                                                         |                                                                                                                                                      cACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGt                                                                                                                                                       <  1:649034/151‑1 (MQ=255)  aCCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGt                                                                                                                                                       >  2:148872/1‑150 (MQ=255)  aCCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGt                                                                                                                                                       >  1:244460/1‑150 (MQ=255)       cGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAAc                                                                                                                                                  <  1:200503/150‑1 (MQ=255)        ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAAccg                                                                                                                                                >  1:39629/1‑151 (MQ=255)             gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCt                                                                                                                                           >  1:898959/1‑151 (MQ=255)                   cGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTg                                                                                                                                     >  1:1002378/1‑151 (MQ=255)                    gAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGc                                                                                                                                    >  1:228695/1‑151 (MQ=255)                      aTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCtt                                                                                                                                  >  1:959428/1‑151 (MQ=255)                           aaGAGTTGTTAGAAAAAAGGCTGGCCGCTGGCATGAACGTTTGGCTTTTGACCTTCTCCCAGGGTGACTAGGAAGAACACGGT‑AAGCCCATCAAAACTCTGCGAAAGGTGATGAGAAAACCG‑GGAAAACCCGCCGCTATCCTGCTTGATAc                                                                                                                             <  2:965209/151‑1 (MQ=255)                           aaGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATAc                                                                                                                             >  2:348113/1‑151 (MQ=255)                              agaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                           >  1:410722/1‑150 (MQ=255)                               gaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCaa                                                                                                                          >  1:638934/1‑150 (MQ=255)                                   ttAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTc                                                                                                                     >  1:423170/1‑151 (MQ=255)                                     agaaaaaaTGCGGGACGCGGGCATGAGCGTAAGGCGCGTGACCTTCTCTGAGGTTGACTAGGGAAAACACGGC‑CGGCGCATTCGGAATCTGCGAAACGGGAGAAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAgg                                                                                                                       <  1:543592/143‑1 (MQ=255)                                             tGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCg                                                                                                           >  1:781143/1‑151 (MQ=255)                                                 ggACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCAcc                                                                                                       >  1:70627/1‑151 (MQ=255)                                                     gCTGGCAGAATCGTTAGGCTTCTAAACTCCTCTCATGGTGACTATGAAGAACCCGTT‑CCGCGCATTAAGATTCTGCGCAACGTGATAAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATg                                                                                                    <  2:186451/150‑1 (MQ=255)                                                          aaTGAACGTTAGGCGTCTGAACTTCTCTAATGGTGACTATGAAGAACACGGT‑CAGCGCATTAAAAATCTGCACAACGTGATGAAAAAAACG‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTg                                                                                              <  1:741380/150‑1 (MQ=255)                                                           aTGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                                             >  1:611224/1‑151 (MQ=255)                                                              aaCTTTATGCGTCTGACCTTCTCTCATGGTGACTATGAAAAACACGGT‑CAGCGAATTAAGAATCTGCGAAACGTGATGAAAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                                            <  1:947117/149‑1 (MQ=255)                                                              aaCGTTATGCGTCTGAACTTCTCTCATGGTGACTATGAAGAACACGGT‑CGGCGCATTCAAAATCGGCGCAACGTGATGAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                                            <  1:264112/149‑1 (MQ=255)                                                                cGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                        >  1:944346/1‑151 (MQ=255)                                                                  ttAGGCTTCTGACCTTCTCTCATGGTGACTATGCAGAACGCGGT‑CAGCGCATTAAAATTCGGCGCAAGGTGATGAGAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                        <  1:694065/149‑1 (MQ=255)                                                                   tATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                     >  1:973934/1‑151 (MQ=255)                                                                   tATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                     >  1:286251/1‑151 (MQ=255)                                                                    aggCGTCGAAACTTCCCTCATGGTGACAATGAAAAACCCGGC‑CACCGAATTAAAAATCTGCCAAAGGTGAGGAGAAAACCT‑GGTAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                     <  1:37131/148‑1 (MQ=255)                                                                    aTGCGTCTGACTTTCTCTCATGGTGACTATGAAAAACACGGC‑GAGCGCATAAAAAATCTGCCCAACGTGAGGA‑AAAAACTGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                     <  1:456322/150‑1 (MQ=255)                                                                    aTGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                     >  1:487094/1‑150 (MQ=255)                                                                    aTGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGt                                                                                     >  1:503023/1‑150 (MQ=255)                                                                         tCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                               >  2:23547/1‑151 (MQ=255)                                                                         tCTGAACTTCTCTCATGGTGACTATGCAAAACACGGC‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAac                                                                                  <  1:225365/148‑1 (MQ=255)                                                                                ttCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                        >  1:712033/1‑151 (MQ=255)                                                                                 tctctcATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                       >  2:499118/1‑151 (MQ=255)                                                                                   tctttgggTGACTAGGCAGACCGCGTC‑GAGCGCATTGAGAACTGGCCCAACGGGAGGAAAAAAATT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                           <  1:328293/139‑1 (MQ=255)                                                                                      cATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaa                                                                    >  2:812857/1‑149 (MQ=255)                                                                                         ggTGACTATGCAGAACACGGT‑CAGCGCATTCAAAATCTGCGCACCGTGATGAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                               <  1:617847/151‑1 (MQ=255)                                                                                         ggTGACTAAGCAAAACACGGA‑CAGCGCATTCAGAATCGGCGCAAGGTGAGGAAAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                                 <  1:267830/149‑1 (MQ=255)                                                                                          gTGACTATCCAGAAC‑CGGTCCAGCGCATTCAAATTCTGCGCAACGTGATAAAAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                               <  2:802853/150‑1 (MQ=255)                                                                                           tGACTATGCAAAACCCGGT‑CGGCGCATTCAGAATCGGCGAAGCGTGATGAAAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTgg                                                              <  1:841224/150‑1 (MQ=255)                                                                                             ccTATGCAGAACACGGT‑CAGTGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTc                                                            >  2:672869/2‑150 (MQ=255)                                                                                                atgaaaaaCCCGGT‑CAGCGCATACAAATCCTGCGCAACGTGATGAAAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGa                                                         <  2:248418/144‑1 (MQ=255)                                                                                                    aGAACACGGT‑CAGCGCATTAAGAATCTGCGCAACGTAATGAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                      <  2:468048/149‑1 (MQ=255)                                                                                                      aaCACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                                  >  2:640144/1‑151 (MQ=255)                                                                                                      aaCACGGT‑CAGCGATTTAAAAATCGGCCCAACGTAAGAAGAAAAACT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTa                                                   <  1:165130/150‑1 (MQ=255)                                                                                                        cacGGT‑CAGCGAATTCAGAATCTGCGCAAGGTGATGAAAAAACCT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTa                                                   <  1:504819/148‑1 (MQ=255)                                                                                                                aGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                         >  2:333920/1‑151 (MQ=255)                                                                                                                    cATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                     >  2:1043765/1‑151 (MQ=255)                                                                                                                      ttCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATa                                    >  1:4907/1‑150 (MQ=255)                                                                                                                      ttCAAATTCGGCGAACCGTGAGGAGAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                     <  1:448626/149‑1 (MQ=255)                                                                                                                         atattccgcgcAACGCGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCAGCACCATTAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGCCAGACCTTTACTTTCCACACTGATaa                                   >  2:570249/8‑148 (MQ=255)                                                                                                                         agattttgaaaagcGTGAAGAAAAAATCT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGa                                      <  2:418070/132‑1 (MQ=255)                                                                                                                           aagcTGCCCAATGTAAAGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                               <  1:910967/147‑1 (MQ=255)                                                                                                                            aTCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGt                             >  2:320373/1‑151 (MQ=255)                                                                                                                               tGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTAt                          >  2:820352/1‑151 (MQ=255)                                                                                                                                  gcAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCgg                       >  1:466192/1‑151 (MQ=255)                                                                                                                                    aaggTGATGAAAAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                      <  1:398994/147‑1 (MQ=255)                                                                                                                                    aaCGAGATGAAAAAAACT‑GAAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                      <  2:1068007/150‑1 (MQ=255)                                                                                                                                     aCGTGATAAGCAAACAT‑GTAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCgg                       <  2:979754/148‑1 (MQ=255)                                                                                                                                      cGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc                   >  2:520798/1‑151 (MQ=255)                                                                                                                                          aggaGCAAAACT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACa                  <  2:787407/146‑1 (MQ=255)                                                                                                                                           gaaacAAAGCT‑GAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGc                <  2:4907/145‑1 (MQ=255)                                                                                                                                             agaaaacCT‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGc                <  1:1055262/144‑1 (MQ=255)                                                                                                                                               cAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg          >  1:168268/1‑151 (MQ=255)                                                                                                                                                aaaaCT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg         >  2:21907/1‑151 (MQ=255)                                                                                                                                                aaaaCG‑GGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg          <  1:940533/150‑1 (MQ=255)                                                                                                                                                 aaaCT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg         >  2:254851/1‑150 (MQ=255)                                                                                                                                                  aaCT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt       >  1:347276/1‑151 (MQ=255)                                                                                                                                                         tAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGt  >  1:687841/1‑150 (MQ=255)                                                                                                                                                         tAAAAACGCTGATATAATGTTTGATACCAAAGGTCCGGAAATCAGAACAATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTAAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGAAACAGCGAAATGGTTGCGGt  >  2:101866/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                      CACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGT‑CAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACT‑GGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755721‑1756020

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