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	98.7%	291.6 / ‑3.2	75	coding (793/1413 nt)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (32/42);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTG  >  NC_000913/1756341‑1756634                                                                                                                                                       |                                                                                                                                                  tCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                                                                                                                                    >  1:1059622/1‑150 (MQ=255) tCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                                                                                                                                    >  2:106679/1‑150 (MQ=255)  ccACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa                                                                                                                                                 <  2:217352/151‑1 (MQ=255)  ccACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑a                                                                                                                                                  >  2:684631/1‑150 (MQ=255)   cacaTCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa‑a                                                                                                                                               <  2:713868/151‑1 (MQ=255)    acaTCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa‑a                                                                                                                                               >  2:240736/1‑150 (MQ=255)     catcatCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa‑aaa                                                                                                                                             <  1:1095812/151‑1 (MQ=255)     catcatCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa‑aa                                                                                                                                              >  2:971902/1‑150 (MQ=255)     catcatCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa‑aa                                                                                                                                              >  1:579676/1‑150 (MQ=255)      atcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa‑aaa                                                                                                                                             >  1:406396/1‑150 (MQ=255)         atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑aa‑aaaatgt                                                                                                                                         >  1:642265/1‑147 (MQ=255)              cAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCg                                                                                                                                     <  2:441384/151‑1 (MQ=255)               aaaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGt                                                                                                                                    <  1:154486/151‑1 (MQ=255)               aaaaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGt                                                                                                                                    <  2:760397/151‑1 (MQ=255)                 aaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGc                                                                                                                                  <  1:291996/151‑1 (MQ=255)                    cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCAc                                                                                                                                <  1:47760/150‑1 (MQ=255)                    cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCAc                                                                                                                                <  1:97212/150‑1 (MQ=255)                    cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACg                                                                                                                               <  1:228240/151‑1 (MQ=255)                        aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTaaa                                                                                                                           <  1:315501/151‑1 (MQ=255)                        aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTaaa                                                                                                                           <  2:1009429/151‑1 (MQ=255)                        aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTaa                                                                                                                            >  1:249541/1‑150 (MQ=255)                        aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTaa                                                                                                                            >  1:245176/1‑150 (MQ=255)                            aGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCg                                                                                                                       <  2:510083/151‑1 (MQ=255)                             ggAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGt                                                                                                                      >  2:952622/1‑151 (MQ=255)                                aGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTAt                                                                                                                   >  2:1037160/1‑151 (MQ=255)                                  gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATc                                                                                                                  >  2:778291/1‑150 (MQ=255)                                   ccTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCa                                                                                                                 >  2:790929/1‑150 (MQ=255)                                      caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACt                                                                                                               <  2:250560/149‑1 (MQ=255)                                      caacaaCTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACt                                                                                                               <  2:504003/149‑1 (MQ=255)                                      caacaaCTTCGACGAAAGCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCAc                                                                                                                <  2:555785/148‑1 (MQ=255)                                       aacaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCg                                                                                                            <  2:866068/151‑1 (MQ=255)                                        acaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCg                                                                                                            >  1:866936/1‑150 (MQ=255)                                        acaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCg                                                                                                            >  1:426093/1‑150 (MQ=255)                                           acTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGAccc                                                                                                        <  1:1144319/151‑1 (MQ=255)                                              tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCa                                                                                                       <  2:582676/149‑1 (MQ=255)                                               cgacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGa                                                                                                     <  2:896520/150‑1 (MQ=255)                                               cgacgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGAt                                                                                                    <  2:653254/151‑1 (MQ=255)                                                 acgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGc                                                                                                  <  2:911759/151‑1 (MQ=255)                                                      aTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGAt                                                                                             >  2:88943/1‑151 (MQ=255)                                                         cTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTcc                                                                                          >  2:742654/1‑151 (MQ=255)                                                           cGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCAt                                                                                        <  2:238532/151‑1 (MQ=255)                                                            gAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATg                                                                                       >  2:1080362/1‑151 (MQ=255)                                                             aaGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGa                                                                                      <  2:1120613/151‑1 (MQ=255)                                                             aaGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGa                                                                                      <  2:6175/151‑1 (MQ=255)                                                                  tctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaa                                                                                  >  2:798603/1‑150 (MQ=255)                                                                   ctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaaa                                                                                 >  2:308311/1‑150 (MQ=255)                                                                     gACGGCATCCTGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAACTTATCTTCGCCCAGAAGATGATGATCG‑AAGAAATGTATCCGTGCACGCAACGTCGTTATCACTGCGAGCCAGATGCTGGATTCGATGATCaaca                                                                                >  2:716082/1‑147 (MQ=255)                                                                     gACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaaaaa                                                                               <  2:815292/150‑1 (MQ=255)                                                                        ggCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCa                                                                           <  1:780636/151‑1 (MQ=255)                                                                          catcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACg                                                                         >  1:458595/1‑151 (MQ=255)                                                                           atcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGc                                                                        <  2:184974/151‑1 (MQ=255)                                                                              atGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                                     >  2:276239/1‑151 (MQ=255)                                                                                       cgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCa                                                            <  1:847034/151‑1 (MQ=255)                                                                                          ggCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGa                                                          >  2:501817/1‑150 (MQ=255)                                                                                                  gggTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                                   <  2:388811/149‑1 (MQ=255)                                                                                                     tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                              >  2:44813/1‑151 (MQ=255)                                                                                                     tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                              >  2:1093537/1‑151 (MQ=255)                                                                                                      gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGtt                                             <  2:399099/151‑1 (MQ=255)                                                                                                      gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                              <  2:975994/150‑1 (MQ=255)                                                                                                           aaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCaaa                                        <  2:231746/151‑1 (MQ=255)                                                                                                            aTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                                       >  1:584443/1‑151 (MQ=255)                                                                                                              cccGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                                     >  2:51704/1‑151 (MQ=255)                                                                                                              cccGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                                     <  2:907110/151‑1 (MQ=255)                                                                                                                     gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTc                              >  2:651359/1‑151 (MQ=255)                                                                                                                     gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTc                              <  1:1003357/151‑1 (MQ=255)                                                                                                                     gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTc                              <  2:687819/151‑1 (MQ=255)                                                                                                                        gaagTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                           <  1:81506/151‑1 (MQ=255)                                                                                                                          agTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACgg                         <  2:848902/151‑1 (MQ=255)                                                                                                                           gTTATCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGt                        >  1:888312/1‑151 (MQ=255)                                                                                                                               tCTTCGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGGGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTAc                      >  1:1104207/1‑149 (MQ=255)                                                                                                                                   cGCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                >  1:217109/1‑151 (MQ=255)                                                                                                                                    gCCCAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa               <  2:444420/151‑1 (MQ=255)                                                                                                                                     cccAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                >  2:179885/1‑149 (MQ=255)                                                                                                                                     cccAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa               <  1:1119752/150‑1 (MQ=255)                                                                                                                                      ccAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt             <  1:157177/151‑1 (MQ=255)                                                                                                                                      ccAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt             <  2:227333/151‑1 (MQ=255)                                                                                                                                       cAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt             >  1:1030397/1‑150 (MQ=255)                                                                                                                                       cAGAAGATGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg            <  1:972034/151‑1 (MQ=255)                                                                                                                                        agaagaTGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGa           <  2:843301/151‑1 (MQ=255)                                                                                                                                         gaagaTGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTACAAACGCCATCCTCGACGGTACTGACGCAGTGAt          >  2:368863/1‑151 (MQ=255)                                                                                                                                           agaTGATGATCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc        >  2:172455/1‑151 (MQ=255)                                                                                                                                            gatgatgatCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCt       >  2:59963/1‑151 (MQ=255)                                                                                                                                             atgatgatCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTg      >  2:1080856/1‑151 (MQ=255)                                                                                                                                                atgatCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt   >  1:467570/1‑151 (MQ=255)                                                                                                                                                atgatCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt   >  2:52338/1‑151 (MQ=255)                                                                                                                                                  gatCGAAA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg  <  1:106698/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  TCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCG‑AA‑AAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTG  >  NC_000913/1756341‑1756634

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