Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,529	A→T	T278S (ACC→TCC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,529	0	A	T	98.7%	272.8 / ‑4.6	74	T278S (ACC→TCC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (36/37);  minor base C (0/1);  total (36/38)
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.75e-01

CAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTT  >  NC_000913/1756381‑1756673                                                                                                                                                     |                                                                                                                                                 caacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTcc                                                                                                                                                <  1:925020/151‑1 (MQ=255)      tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGAt                                                                                                                                           <  1:408966/151‑1 (MQ=255)      tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGAt                                                                                                                                           <  2:253960/151‑1 (MQ=255)             aaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAt                                                                                                                                    <  2:951183/151‑1 (MQ=255)                ccTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTcc                                                                                                                                 <  1:1145754/151‑1 (MQ=255)                   cGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCAt                                                                                                                               <  1:640407/150‑1 (MQ=255)                   cGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATg                                                                                                                              <  1:701544/151‑1 (MQ=255)                       gCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCa                                                                                                                          >  1:1131160/1‑151 (MQ=255)                       gCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCa                                                                                                                          >  2:673682/1‑151 (MQ=255)                       gCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCa                                                                                                                          >  1:75100/1‑151 (MQ=255)                              aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAcc                                                                                                                    <  2:733893/150‑1 (MQ=255)                              aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAcc                                                                                                                    <  2:1087358/150‑1 (MQ=255)                              aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAcc                                                                                                                    <  2:1041545/150‑1 (MQ=255)                              aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAcc                                                                                                                    <  1:205804/150‑1 (MQ=255)                                  catcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACg                                                                                                                >  2:982666/1‑150 (MQ=255)                                    tcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCACCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGcc                                                                                                              >  2:658125/1‑150 (MQ=255)                                         gTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                                                                                                         <  2:524793/150‑1 (MQ=255)                                         gTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                                                                                                         <  2:823821/150‑1 (MQ=255)                                               cgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAg                                                                                                  >  1:280248/1‑151 (MQ=255)                                                 tGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGaa                                                                                                >  1:60243/1‑151 (MQ=255)                                                   gCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAg                                                                                               <  2:69259/150‑1 (MQ=255)                                                    cGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGc                                                                                              <  1:1089422/150‑1 (MQ=255)                                                      aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                                           >  2:144763/1‑151 (MQ=255)                                                          gggTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                                                        >  2:87153/1‑150 (MQ=255)                                                          gggTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                                                        >  1:1062539/1‑150 (MQ=255)                                                            gtgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                                                                     >  2:582227/1‑151 (MQ=255)                                                              gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                                                   >  2:518558/1‑151 (MQ=255)                                                               tAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                                                   <  2:66086/150‑1 (MQ=255)                                                                    aTCTCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                                                                             >  2:432561/1‑151 (MQ=255)                                                                     tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                                                                            <  1:592615/151‑1 (MQ=255)                                                                           tAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAACCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                                                       <  2:72345/150‑1 (MQ=255)                                                                           tAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                                                       <  1:780305/150‑1 (MQ=255)                                                                           tAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                                                       <  1:728883/150‑1 (MQ=255)                                                                           tAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCt                                                                      <  1:189576/151‑1 (MQ=255)                                                                                 aagTTATCTTTGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCAATCCTCGACgg                                                                >  2:876918/1‑151 (MQ=255)                                                                                 aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                                                  <  1:125049/149‑1 (MQ=255)                                                                                 aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                                                  <  1:1083745/149‑1 (MQ=255)                                                                                 aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                                                  <  2:1130635/149‑1 (MQ=255)                                                                                   gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTa                                                              >  1:390203/1‑151 (MQ=255)                                                                                   gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCACCCCCGACGGTa                                                              >  2:767673/1‑151 (MQ=255)                                                                                        cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                                                         <  1:24875/151‑1 (MQ=255)                                                                                          tCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGAAAAGTCGTTATCACGGCGTCCCAGATGCTGGATTCCATGATAAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                                                       <  2:163210/151‑1 (MQ=255)                                                                                            gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                                                     >  2:286744/1‑151 (MQ=255)                                                                                                 gaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                >  1:1097461/1‑151 (MQ=255)                                                                                                 gaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                >  1:541629/1‑151 (MQ=255)                                                                                                    gatgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                >  1:1188064/1‑148 (MQ=255)                                                                                                    gatgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                >  1:589376/1‑148 (MQ=255)                                                                                                     atgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                                                >  2:1053777/1‑147 (MQ=255)                                                                                                     atgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTg                                             <  1:205499/150‑1 (MQ=255)                                                                                                         tgatCGACATTTGTATCCGTTCACGTAAAGTCGCTATCACTGCGTCCCAGCTGCTGGATTCCATTATCAAAAACCCCCGCCCCACTCGCGCAGAAACCGGCGTCGTTGCAAACGCCATCGTCGACAGAACTGACGCAGTGATGCTGTCTgg                                        >  2:842597/1‑151 (MQ=255)                                                                                                          catCGAAAAATGTATCCGTGCACGTAAAGACGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGGGATGCTGTCTGGt                                       <  2:1189479/150‑1 (MQ=255)                                                                                                           atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCCTCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg                                      >  2:528472/1‑151 (MQ=255)                                                                                                           atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg                                      <  1:1115444/151‑1 (MQ=255)                                                                                                           atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg                                      <  2:67566/151‑1 (MQ=255)                                                                                                            tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGCCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGa                                     <  1:994161/151‑1 (MQ=255)                                                                                                                 aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACg                                                        <  1:310271/127‑1 (MQ=255)                                                                                                                 aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCg                                >  1:504373/1‑151 (MQ=255)                                                                                                                 aaaTGTACCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAACCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACg                                                        >  2:310177/1‑127 (MQ=255)                                                                                                                   aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa                              >  2:660839/1‑151 (MQ=255)                                                                                                                   aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa                              >  1:761889/1‑151 (MQ=255)                                                                                                                   aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa                              >  1:545785/1‑151 (MQ=255)                                                                                                                      tATCCGTGCACGTAAAGTCGTTATCACTGCGTCCGATATGCTGGATCCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGATGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaa                            <  2:572758/150‑1 (MQ=255)                                                                                                                              cACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATa                   <  2:190012/151‑1 (MQ=255)                                                                                                                               aCGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                  <  1:1014796/151‑1 (MQ=255)                                                                                                                               aCGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                  <  1:824277/151‑1 (MQ=255)                                                                                                                                cGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAcc                 <  1:415073/151‑1 (MQ=255)                                                                                                                                 gTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCCTCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAcc                 >  2:918395/1‑150 (MQ=255)                                                                                                                                      gTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTgg           >  2:718304/1‑151 (MQ=255)                                                                                                                                        cGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGaa         >  1:15308/1‑151 (MQ=255)                                                                                                                                          ttATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAg        >  1:673153/1‑150 (MQ=255)                                                                                                                                          ttATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGc       <  1:582516/151‑1 (MQ=255)                                                                                                                                              cACTGCGTCCCAGGTGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCTGGTGACGGTGCAAACGCCATCCTCTACGGTACCGACCCAGAGATGCTGTCTGGTGAATGCCCAAGAGGTCAATACCCGCTGGAAGCGGtt   >  2:904231/1‑151 (MQ=255)                                                                                                                                               aCTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGttt  >  2:618521/1‑151 (MQ=255)                                                                                                                                                cTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGttt  >  1:513810/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTT  >  NC_000913/1756381‑1756673

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