Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,377	T→C	L227P (CTC→CCC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,377	0	T	C	96.0%	259.1 / ‑3.3	75	L227P (CTC→CCC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (29/43);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCAC  >  NC_000913/1756231‑1756524                                                                                                                                                   |                                                                                                                                                    gATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCca                                                                                                                                                   <  1:796491/149‑1 (MQ=255)     ttttGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACt                                                                                                                                             >  2:538497/1‑151 (MQ=255)      ttGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCcaaca                                                                                                                                                <  2:177789/147‑1 (MQ=255)      ttGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCcaac                                                                                                                                                 <  1:920397/146‑1 (MQ=255)       tGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCcaaca                                                                                                                                                <  2:181330/146‑1 (MQ=255)           tGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTcgacg                                                                                                                                       >  1:991079/1‑151 (MQ=255)            gCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTcgacga                                                                                                                                      <  1:286163/151‑1 (MQ=255)             cGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTcgacga                                                                                                                                      >  2:601755/1‑150 (MQ=255)              gAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGaaa                                                                                                                                    <  1:221174/151‑1 (MQ=255)              gAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGaa                                                                                                                                     >  1:1018084/1‑150 (MQ=255)                aCAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATc                                                                                                                                  >  1:58627/1‑151 (MQ=255)                 cAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATcc                                                                                                                                 <  1:1001091/151‑1 (MQ=255)                  aaGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCt                                                                                                                                <  2:23528/151‑1 (MQ=255)                  aaGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCt                                                                                                                                <  1:894669/151‑1 (MQ=255)                  aaGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCt                                                                                                                                <  2:1058993/151‑1 (MQ=255)                   aGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTc                                                                                                                               <  1:369740/151‑1 (MQ=255)                              ttgttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTg                                                                                                                     <  2:156720/150‑1 (MQ=255)                              ttgttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTg                                                                                                                     >  1:347561/1‑150 (MQ=255)                                gttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACg                                                                                                                  >  1:334345/1‑151 (MQ=255)                                gttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACg                                                                                                                  <  2:643831/151‑1 (MQ=255)                                 ttgCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAACCCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCAGa                                                                                                                    >  2:398238/1‑148 (MQ=255)                                     tgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatc                                                                                                             >  2:932875/1‑151 (MQ=255)                                       ctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatcat                                                                                                           >  1:388105/1‑151 (MQ=255)                                        ttCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATg                                                                                                          >  1:314027/1‑151 (MQ=255)                                         tCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGcatcat                                                                                                           <  1:113802/149‑1 (MQ=255)                                         tCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATgg                                                                                                         >  1:469856/1‑151 (MQ=255)                                          ccTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGt                                                                                                        <  1:559302/151‑1 (MQ=255)                                             ttATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGtt                                                                                                       <  2:768085/149‑1 (MQ=255)                                                  cGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGt                                                                                                 <  2:226144/150‑1 (MQ=255)                                                  cGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTg                                                                                                <  1:160615/151‑1 (MQ=255)                                                     aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCg                                                                                             <  2:679390/151‑1 (MQ=255)                                                      aGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCg                                                                                             >  2:420029/1‑150 (MQ=255)                                                         gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                                                                                         >  2:31480/1‑151 (MQ=255)                                                         gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                                                                                         >  1:805674/1‑151 (MQ=255)                                                         gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                                                                                         <  1:646355/151‑1 (MQ=255)                                                             tGACGTTATCGAAATCCGTAAGCACCTGTCAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGCCGAAATCCTCGAAGCCTCTGACGGCATCAGGGCTGCGCGTGGCGACCTgg                                                                                       >  2:697616/1‑149 (MQ=255)                                                               aCGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgtgt                                                                                   >  2:322979/1‑151 (MQ=255)                                                                cGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTa                                                                                  <  1:94389/151‑1 (MQ=255)                                                                   tATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGa                                                                                >  2:786181/1‑150 (MQ=255)                                                                    aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCTAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaa                                                                               >  1:499839/1‑150 (MQ=255)                                                                    aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaa                                                                               <  2:930575/150‑1 (MQ=255)                                                                    aTCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGaa                                                                               <  1:57683/150‑1 (MQ=255)                                                                      cGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAAt                                                                             <  2:966162/150‑1 (MQ=255)                                                                         aaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCg                                                                         <  2:805628/151‑1 (MQ=255)                                                                            ccGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTa                                                                      <  2:573332/151‑1 (MQ=255)                                                                             cGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGt                                                                       >  1:715107/1‑149 (MQ=255)                                                                             cGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTag                                                                     <  2:284140/151‑1 (MQ=255)                                                                                 aGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTagaaga                                                                 <  2:608415/151‑1 (MQ=255)                                                                                  gCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAgaagaa                                                                <  1:16723/151‑1 (MQ=255)                                                                                     ccTGAATGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAACCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGTCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGtt                                                             <  2:113066/151‑1 (MQ=255)                                                                                      cTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGtt                                                             >  1:216274/1‑150 (MQ=255)                                                                                                   gcggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCagaa                                               <  2:873308/151‑1 (MQ=255)                                                                                                    cggcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCaga                                                >  1:753087/1‑149 (MQ=255)                                                                                                      gcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGGGACTGGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTAGACCa                                                  >  1:965507/1‑145 (MQ=255)                                                                                                      gcgAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgat                                            <  2:453021/151‑1 (MQ=255)                                                                                                           aaCATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatga                                       <  1:385559/151‑1 (MQ=255)                                                                                                            aCATCCACATCAACTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatgat                                      <  1:286288/151‑1 (MQ=255)                                                                                                             cATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATc                                     <  2:176700/151‑1 (MQ=255)                                                                                                              aTCCCCATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                    <  2:68310/151‑1 (MQ=255)                                                                                                              aTCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                    >  1:14986/1‑151 (MQ=255)                                                                                                              aTCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                    <  1:689294/151‑1 (MQ=255)                                                                                                               tCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                    >  2:641250/1‑150 (MQ=255)                                                                                                                  acaTCATCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaa                                 >  1:1023152/1‑150 (MQ=255)                                                                                                                   catcatCTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaa                                 <  2:901102/149‑1 (MQ=255)                                                                                                                    atcatcTCCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAt                              >  1:42304/1‑151 (MQ=255)                                                                                                                         ctcCAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATc                         >  2:277708/1‑151 (MQ=255)                                                                                                                           ccAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTAt                          <  1:311636/148‑1 (MQ=255)                                                                                                                           ccAAAATCGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCAGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATcc                        <  2:616355/150‑1 (MQ=255)                                                                                                                                  cGAAAACCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGt                >  1:760753/1‑151 (MQ=255)                                                                                                                                    aaaaCCAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTaa              <  2:531077/151‑1 (MQ=255)                                                                                                                                        ccAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGt           <  2:496417/150‑1 (MQ=255)                                                                                                                                        ccAGGAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTc          <  1:1034214/151‑1 (MQ=255)                                                                                                                                            gAAGGCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTa      <  1:1065182/151‑1 (MQ=255)                                                                                                                                              aGGCCTCCAGCACTTCGACGGCATCGTCGAGGGCTCGGGCGGTATCATGGTGGCGCGTGGAGACTTGTGTGTAGAAATCTTGGTAGAATAAGTTATGTTTGCCCAGAAGATGATGATCGAAAAATGTAGCCGTGCACGTAAAGTCGTTAt     <  2:182147/150‑1 (MQ=255)                                                                                                                                                gCCCCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCAc  >  1:282612/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GATCTTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCAC  >  NC_000913/1756231‑1756524

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