Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,862	0	C	A	49.2%	82.5 / 43.9	63	S55R (AGC→AGA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base C (29/3);  new base A (0/31);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 7.14e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGA  >  NC_000913/1755712‑1756008                                                                                                                                                        |                                                                                                                                                   aaTTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGc                                                                                                                                                    >  1:1209769/1‑151 (MQ=255)       ttGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACt                                                                                                                                              >  1:596715/1‑151 (MQ=255)       ttGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACt                                                                                                                                              >  2:1318259/1‑151 (MQ=255)                 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)                                                                                                                                                      accAAATCTGATAAAGCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGa  <  2:1415494/147‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   AATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTG‑ATGAGCAAAACTGGTAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGA  >  NC_000913/1755712‑1756008

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