Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,495	(A)5→6	coding (798/1413 nt)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,490	1	.	A	96.8%	236.3 / 0.8	62	coding (793/1413 nt)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (30/30);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-01

CCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTG  >  NC_000913/1756342‑1756630                                                                                                                                                      |                                                                                                                                              ccACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑a                                                                                                                                              <  1:369158/150‑1 (MQ=255)     atcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aaa‑aaa                                                                                                                                        <  2:460029/151‑2 (MQ=255)          ctcCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTAt                                                                                                                                    >  1:461637/1‑151 (MQ=255)           tcCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATc                                                                                                                                   >  2:241801/1‑151 (MQ=255)            ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTAt                                                                                                                                    <  2:64992/149‑1 (MQ=255)            ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATcc                                                                                                                                  <  1:359260/151‑1 (MQ=255)            ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATc                                                                                                                                   <  2:896068/150‑1 (MQ=255)            ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aaa‑aaatgta                                                                                                                                    <  1:599799/148‑6 (MQ=255)             cAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCg                                                                                                                                 >  1:301230/1‑151 (MQ=255)              aaaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGt                                                                                                                                <  1:16002/151‑1 (MQ=255)                aaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGc                                                                                                                              >  2:842264/1‑151 (MQ=255)                 aTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCa                                                                                                                             <  2:306067/151‑1 (MQ=255)                    gAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGt                                                                                                                          >  2:578447/1‑151 (MQ=255)                     aaaaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTa                                                                                                                         <  2:128576/151‑1 (MQ=255)                       aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTaa                                                                                                                        >  2:777757/1‑150 (MQ=255)                          cAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTc                                                                                                                    <  1:846985/151‑1 (MQ=255)                               aGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAAGAATGTATCCGTGCACGTAAAGTCGTTAt                                                                                                               <  1:887048/151‑1 (MQ=255)                                 gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCa                                                                                                             <  1:487067/151‑1 (MQ=255)                                    tCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTg                                                                                                          >  2:331170/1‑151 (MQ=255)                                     caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACt                                                                                                           <  1:737170/149‑1 (MQ=255)                                     caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGc                                                                                                         <  1:166960/151‑1 (MQ=255)                                     caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGc                                                                                                         <  2:683374/151‑1 (MQ=255)                                     caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGc                                                                                                         <  2:742567/151‑1 (MQ=255)                                      aacaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCg                                                                                                        <  2:894715/151‑1 (MQ=255)                                             tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGAccc                                                                                                    >  1:496870/1‑148 (MQ=255)                                             tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCa                                                                                                   <  1:483566/149‑1 (MQ=255)                                              cgacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGAccc                                                                                                    >  1:279519/1‑147 (MQ=255)                                              cgacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGa                                                                                                 <  1:332397/150‑1 (MQ=255)                                               gacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATg                                                                                               >  1:897436/1‑151 (MQ=255)                                                    aaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTgg                                                                                           >  1:684287/1‑150 (MQ=255)                                                           gAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATg                                                                                   <  2:540426/151‑1 (MQ=255)                                                           gAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATg                                                                                   >  1:391198/1‑151 (MQ=255)                                                             aGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGAt                                                                                 >  1:817352/1‑151 (MQ=255)                                                               ccTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCa                                                                               <  1:456183/151‑1 (MQ=255)                                                                 tctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaa                                                                              >  2:857163/1‑150 (MQ=255)                                                                  ctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaaa                                                                             >  1:699558/1‑150 (MQ=255)                                                                    gACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAc                                                                          >  2:952763/1‑151 (MQ=255)                                                                     aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGATAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGTCTTCCATGATCaaaaa                                                                           <  2:385028/149‑1 (MQ=255)                                                                          atcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGc                                                                    >  1:477286/1‑151 (MQ=255)                                                                                     gcgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCgcgc                                                         <  2:68006/151‑1 (MQ=255)                                                                                     gcgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCgcgc                                                         >  1:284835/1‑151 (MQ=255)                                                                                      cgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCa                                                        <  2:830646/151‑1 (MQ=255)                                                                                         ggcggcCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAg                                                    <  1:924049/151‑1 (MQ=255)                                                                                            gACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGcc                                                  <  2:478140/151‑1 (MQ=255)                                                                                            gACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGcc                                                  >  1:483597/1‑151 (MQ=255)                                                                                             aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                 <  1:233070/151‑1 (MQ=255)                                                                                              ccTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                >  1:724494/1‑151 (MQ=255)                                                                                                       aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                        <  2:39985/150‑1 (MQ=255)                                                                                                           aTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTACGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                                   <  2:369965/151‑1 (MQ=255)                                                                                                                   agaagaAGTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCt                           <  2:718330/151‑1 (MQ=255)                                                                                                                     aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTc                          <  1:630889/150‑1 (MQ=255)                                                                                                                     aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                         <  1:129051/151‑1 (MQ=255)                                                                                                                      agaagTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                         >  1:687584/1‑150 (MQ=255)                                                                                                                        aagTTATCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACg                      >  2:581616/1‑151 (MQ=255)                                                                                                                             aTCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                 >  2:127160/1‑151 (MQ=255)                                                                                                                             aTCTTCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                 >  2:47383/1‑151 (MQ=255)                                                                                                                                 tCGCCCAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc              <  1:342383/150‑1 (MQ=255)                                                                                                                                    cccAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc            >  1:774303/1‑149 (MQ=255)                                                                                                                                    cccAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc            >  1:31480/1‑149 (MQ=255)                                                                                                                                      cAGAAGATGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg        >  1:476583/1‑151 (MQ=255)                                                                                                                                        gaagaTGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGAt      >  2:764689/1‑151 (MQ=255)                                                                                                                                          agaTGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc    <  2:774057/151‑1 (MQ=255)                                                                                                                                          agaTGATGATCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc    <  2:483515/151‑1 (MQ=255)                                                                                                                                           gatgatgatCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCt   >  2:350081/1‑151 (MQ=255)                                                                                                                                            atgatgatCGAAAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTg  >  1:198585/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                              CCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AAA‑AATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTG  >  NC_000913/1756342‑1756630

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