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

GGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755728‑1756020                                                                                                                                                 |                                                                                                                                                      ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAAccg                                                                                                                                                <  2:1171368/151‑1 (MQ=255)   aCCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAAccgcc                                                                                                                                              <  1:658328/151‑1 (MQ=255)   aCCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAAccgcc                                                                                                                                              >  1:161067/1‑151 (MQ=255)            cGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTg                                                                                                                                     >  1:739089/1‑151 (MQ=255)             gAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGc                                                                                                                                    >  2:756040/1‑151 (MQ=255)                   gAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGAt                                                                                                                               >  1:389726/1‑150 (MQ=255)                       agaTTTTAGAAAAAATGCTGGACACGGGCTTAACCTTTTTGCTTCTGACCTCCTCTCATGGGAACTAGGAAAAACACGGCCAGCGCATTCAAAACCTGCGCAAGGTG‑ATGAGAAAACCTGGAAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                           <  2:1046602/150‑1 (MQ=255)                       agaTGTTAGAAAAAAGGCTGGACGCTGGAATGAACTTTATGCGTCTGACCTTCTCTCAGGGTTAATACGAAAAACACGATCAGCCCATTCAGAACCTGCGCAGCGGG‑ATGAAAAAACCTGTAAAAACCGCCGCTATCCTGCtt                                                                                                                                  <  1:598010/143‑1 (MQ=255)                          tGTTAGCAAAAAGGCTGGACGCGGGAATGAACGTTATGCGTCTAACCTTCTCTCATGGTGACTAGGCAAACCACGGCCAGCGCATTCAAAATCGGCGCAACGTG‑ATGAGAAAAATTGTAAAAACCGCCGCTATCCTGCTTGATAcc                                                                                                                            <  2:1240046/146‑1 (MQ=255)                          tGTTAGCAAAAAGGCGGGACGCGGGCATGAACTTTATGCGTCTAAACTTCTCTAATGGTGACTATGAAGAACCCGGCCGGCGCATTAAGAATCTGCGCACGGTG‑ATGAGAAAACCGGGAAAAACCGCCGCTATCCTGCTTGATAcc                                                                                                                            <  2:532847/146‑1 (MQ=255)                            ttAGCAAAAAGGCTGGACGCGGGAATGAACTTTATGCTTCTGACCTTCTCTCATGGTGACTAGGAAGAACACGGTCAGCGCATTCAAAATCTGCGCACGGTG‑ATGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCaa                                                                                                                          <  2:161498/146‑1 (MQ=255)                              ataaaaaaGGCTGGACGCGGGAATGAATGTTAGGCGTCTGACTTTCCCTCATGGTACTTAGGAAAAACCCGCCCGGCGCATTAAGATTCTGCCAAACGGG‑AGAAAAAAAATGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCg                                                                                                                   <  1:1229230/147‑1 (MQ=255)                              acaaaaatTGTTGGCCGCGGGAATGACTGTTGTGCGTCTGACCTTCTCTCAGGGTAATTAGGAAGACCAGGGCCGGCCCATTCAAATCCTGCG‑AACGGGAAGGAGAAAAACTGAAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCg                                                                                                                   <  1:475132/147‑1 (MQ=255)                              aGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCg                                                                                                                   >  2:153309/1‑151 (MQ=255)                               acaaaaaGGCTGGACGCGGGAATGAACTTTATGCTTCTGATCTTCTCTCATGGTGATTATGAAGAACACGGTCAGCGCATTAAAATTCTGCCCATCGTG‑ATGAGAAAAATTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCgg                                                                                                                  <  2:1202293/148‑1 (MQ=255)                                cTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCgg                                                                                                                  >  2:768195/1‑150 (MQ=255)                                 tAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGa                                                                                                                 >  1:662732/1‑150 (MQ=255)                                     aTGCTGGCCCCTGGAATTAGTGTTAGGCTTCTGAATTTCTCTCATGGTGACAATGCAAAACGCGGTCAGCACTTTAAGAATCTCCGCAGCTTG‑AGAAGAAAACCGGGAAAAACCGCCGCTATCCTGCTTGATGCCAAAGGTCCGGAAATcc                                                                                                            <  2:909173/151‑1 (MQ=255)                                      tGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATcc                                                                                                            >  2:670840/1‑150 (MQ=255)                                       gcgggccGCTGGAATGAACGTTGTGCGTCTGACTTTCCCTCAGGGTAACTATGAAGAACCCGGTCAGCGCATTAAGAATCTGCGCAAGGGG‑ATGAGAAAAAGTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATcc                                                                                                            <  2:657946/143‑1 (MQ=255)                                       gCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGc                                                                                                          >  2:389823/1‑151 (MQ=255)                                           gACGCTGGAATGAACGTTATGCGTCTGACCTTCTCTCGTGGTGACTATGCAAAACACGGCCAGCGAATTCAGAATCGGCGCACGGTG‑AGGAGAAAAACGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCa                                                                                                      <  2:448542/151‑1 (MQ=255)                                            aCGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCAGCTTGATACCAAAGGTCCGGAAATCCGCACCAt                                                                                                     >  2:80293/1‑151 (MQ=255)                                                gggCATGAACGTTTTGCTTCTAATCTTCTCTCAGGGTGACTAGGAAAACCCCGGCCAGCCCTTTCAAAATCTCCGCAACGTG‑AAGAGAAAAAGGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGaaa                                                                                                 <  1:1186799/150‑1 (MQ=255)                                                    aTGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTgg                                                                                             >  1:953524/1‑151 (MQ=255)                                                      gAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGa                                                                                            >  1:848182/1‑150 (MQ=255)                                                          gTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                       >  2:1429140/1‑151 (MQ=255)                                                           ttATGCGTCTGAACTTCTCTCATGGTGACTATGCAAAACACGGTCAGCGAATTCAAAATCTGCGCAACGTA‑ATGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                        <  2:1172096/149‑1 (MQ=255)                                                               gCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAac                                                                                  >  2:1351249/1‑151 (MQ=255)                                                               gCGTCACAACTTCTCTCATGGTGATTAGGCAGACCGCGGCCACCGGATTCAAAAT‑TGCGCAAAGGGAATGAGAAAACCGGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTaa                                                                                   <  1:522089/150‑1 (MQ=255)                                                                    tGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGt                                                                             >  2:619628/1‑151 (MQ=255)                                                                     gAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGtt                                                                            >  2:1424958/1‑151 (MQ=255)                                                                         ttCTCTCATGGTAACTATGAAAACCAGGGACAGCGCATTCAAAGTCGGCGAAAGGTG‑AGGAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtc                                                                          <  1:378704/149‑1 (MQ=255)                                                                         ttCTCTAATGGTGATTAGGCAGAACCGGGTCGGCCCATTCAAATTCTGCGAAAGGTG‑ATGAAAAAACTTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                               <  2:212070/144‑1 (MQ=255)                                                                            tctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTg                                                                      >  1:409134/1‑150 (MQ=255)                                                                             ctcATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCTCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGa                                                                     >  2:610245/1‑150 (MQ=255)                                                                              tcATGGTGACTATGAAAAACACGGTCAGCGCATTAAAAATCTGCGCAACGTG‑ATGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtc                                                                          <  2:1256110/144‑1 (MQ=255)                                                                                  gggcactaggaaaaaCACGGCCAGCGCATTCAAAATCGGCGAAACGGG‑AGGAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAg                                                                  <  1:727505/135‑1 (MQ=255)                                                                                  ggTGACTATGAAGAACACGGTCGGCGCATTAAAAACCGGCGCAACGTG‑AGGAGAAAACCTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTtctc                                                                        <  2:332688/142‑1 (MQ=255)                                                                                  ggTGACTATGAAAACCGCGGTCAGCCCATTCAGAATCTGCGAAAGGTA‑AGGAAAAAATCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCt                                                                <  1:394518/150‑1 (MQ=255)                                                                                    tGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                             >  1:917606/1‑151 (MQ=255)                                                                                     gACTATGCACACCCGGGCCAGCGCATTAAGAATCGGCAAAAGGTG‑AGGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                             <  1:351681/150‑1 (MQ=255)                                                                                        taggCAGAACCCGGTCAGCGCATTCAAATTCGGCGCAACGTG‑ATGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAg                                                          <  2:489732/147‑1 (MQ=255)                                                                                         aTGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAc                                                        >  1:688854/1‑151 (MQ=255)                                                                                           gCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCt                                                      >  2:844402/1‑151 (MQ=255)                                                                                             aGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                                    >  1:689641/1‑151 (MQ=255)                                                                                                acccGGTCAGCGCATAAAGAATCTGCGCAAGGTG‑AGAAAAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                                  <  1:77923/147‑1 (MQ=255)                                                                                                acacGGACGGCACATTCAGAAGCTGCGAAGAGTT‑ATATGAAAAATTGGAAAACCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACt                                                 <  1:242124/151‑1 (MQ=255)                                                                                                       cAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcacc                                           >  1:433864/1‑150 (MQ=255)                                                                                                         gcgcttaaaGATTCTGCGAAACGTA‑AGAAAAAAAACGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcacca                                          <  2:923196/141‑1 (MQ=255)                                                                                                                  aatcttGCGCACGGTA‑ATGAGAAAAGTTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                               <  2:372060/146‑1 (MQ=255)                                                                                                                   aaccggcgcAGCGGG‑ATGAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTg                                       <  2:1271658/137‑1 (MQ=255)                                                                                                                   aaTCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTg                              >  2:197338/1‑151 (MQ=255)                                                                                                                    aTCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGt                             >  2:547622/1‑151 (MQ=255)                                                                                                                        gcgcAAGGTG‑ATGAGCAAAATTGATAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATc                                <  2:1334935/144‑1 (MQ=255)                                                                                                                        gcgcAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATc                         >  1:1315922/1‑151 (MQ=255)                                                                                                                         cgaaaggTG‑ATAAGAAAAACTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCt                               <  2:1072438/138‑1 (MQ=255)                                                                                                                          gcACCGTG‑AGGAAAAAAACTGGAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                            <  1:12199/146‑1 (MQ=255)                                                                                                                            aatgTG‑ATGAGAAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                            <  2:476607/141‑1 (MQ=255)                                                                                                                            aaCGTG‑ATAAAAAAAAGTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGc                      <  1:881285/150‑1 (MQ=255)                                                                                                                               gTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAAc                   >  2:852050/1‑150 (MQ=255)                                                                                                                                   agaacaaaaaGTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCg                        <  2:1105928/136‑1 (MQ=255)                                                                                                                                      accAAATCTGATAAAGCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGa              <  2:1415494/147‑1 (MQ=255)                                                                                                                                         aaaaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg         >  2:1272846/1‑151 (MQ=255)                                                                                                                                         aaaaCTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg         >  1:1030440/1‑151 (MQ=255)                                                                                                                                         aaaaCTGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATg          <  2:462983/150‑1 (MQ=255)                                                                                                                                               ggTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCg    >  2:1138331/1‑150 (MQ=255)                                                                                                                                                 tAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGt  >  2:545187/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGT  >  NC_000913/1755728‑1756020

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