Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,756,497	G→T	C267F (TGT→TTT)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,756,497	0	G	T	95.6%	314.7 / 2.8	91	C267F (TGT→TTT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (39/47);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCC  >  NC_000913_3_bme_pgi/1756353‑1756642                                                                                                                                                 |                                                                                                                                                  ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATcc                                                                                                                                              <  3:588663/150‑1 (MQ=255)  cAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTCGAAGCCGTTCTCTTCGCCCAGAGGGTGATGATCGAAGAATTTATCCg                                                                                                                                             >  1:39391/1‑150 (MQ=255)     aaTCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGc                                                                                                                                          >  1:452889/1‑150 (MQ=255)                aGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCg                                                                                                                               <  1:396653/150‑1 (MQ=255)                 ggAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGt                                                                                                                              <  2:361313/150‑1 (MQ=255)                      gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCa                                                                                                                         <  1:753263/150‑1 (MQ=255)                        cTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGt                                                                                                                              >  2:552189/1‑143 (MQ=255)                        cTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGt                                                                                                                              <  1:552189/143‑1 (MQ=255)                         tCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTg                                                                                                                      <  1:197138/150‑1 (MQ=255)                           aacaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCa                                                                                                                         >  1:721690/1‑145 (MQ=255)                           aacaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCa                                                                                                                         <  2:721690/145‑1 (MQ=255)                            acaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGTGTTTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGa                                                                                                                   <  1:18855/150‑1 (MQ=255)                            acaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGa                                                                                                                   >  4:521557/1‑150 (MQ=255)                            acaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGa                                                                                                                   <  1:97645/150‑1 (MQ=255)                             caacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGAc                                                                                                                  >  2:927762/1‑150 (MQ=255)                                cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCa                                                                                                               <  1:918942/150‑1 (MQ=255)                                cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCa                                                                                                               >  2:613609/1‑150 (MQ=255)                                cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCa                                                                                                               <  2:324333/150‑1 (MQ=255)                                cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCa                                                                                                               >  2:994576/1‑150 (MQ=255)                                  tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGc                                                                                                                                          <  1:729903/121‑1 (MQ=255)                                  tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGc                                                                                                                                          >  2:729903/1‑121 (MQ=255)                                  tCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGa                                                                                                             <  1:814654/150‑1 (MQ=255)                                      cgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCt                                                                                                         <  3:314885/150‑1 (MQ=255)                                       gaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTg                                                                                                        >  4:475594/1‑150 (MQ=255)                                        aaaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTgg                                                                                                       <  3:521557/150‑1 (MQ=255)                                         aaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGa                                                                                                      >  2:782076/1‑150 (MQ=255)                                           tCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGAtt                                                                                                    <  2:452889/150‑1 (MQ=255)                                              tCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTAt                                                                                                                           <  3:368494/124‑1 (MQ=255)                                              tCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTAt                                                                                                                           >  4:368494/1‑124 (MQ=255)                                              tCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCa                                                                                                 <  1:927762/150‑1 (MQ=255)                                              tCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCa                                                                                                 >  4:384137/1‑150 (MQ=255)                                              tCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGACCGAAAAATATATCACTGAAAGTAACATAGTTTTCACTGAGACCACGAAACTGGATTCCa                                                                                                 >  1:753808/1‑150 (MQ=255)                                                     ctctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaa                                                                                          >  2:733013/1‑150 (MQ=255)                                                      tctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGt                                                                                                                              <  1:1039445/113‑1 (MQ=255)                                                      tctGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGt                                                                                                                              >  2:1039445/1‑113 (MQ=255)                                                       ctGACGGCATCATGGTAGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTATAGTGGTTATCACGCCGATCAATATGCTGGATACCTTGATCaaaa                                                                                        >  3:519840/1‑150 (MQ=255)                                                        tGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCaaaaa                                                                                       >  3:413429/1‑150 (MQ=255)                                                          aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                                                                     >  3:596090/1‑150 (MQ=255)                                                           ccgcttcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAccc                                                                                    <  4:519840/145‑1 (MQ=255)                                                               atcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGc                                                                                >  3:491158/1‑150 (MQ=255)                                                                  atGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                                             >  3:521154/1‑150 (MQ=255)                                                                  atGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGAACGTCAAGTAGTTATCACCGCGACCCAGATGCTGGATTGCATGATAAAAAACACACGCCCg                                                                             >  3:298467/1‑150 (MQ=255)                                                                      ttGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTc                                                                         >  4:220345/1‑150 (MQ=255)                                                                        gcgcGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgc                                                                       <  1:491654/150‑1 (MQ=255)                                                                          gcgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTa                                                                                                                                     <  3:255414/86‑1 (MQ=255)                                                                          gcgTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTa                                                                                                                                     >  4:255414/1‑86 (MQ=255)                                                                                cGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGc                                                               >  2:950085/1‑150 (MQ=255)                                                                                  aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                             <  4:625205/150‑1 (MQ=255)                                                                                  aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                             <  3:220345/150‑1 (MQ=255)                                                                                  aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                             <  1:950085/150‑1 (MQ=255)                                                                                  aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                             <  3:272596/150‑1 (MQ=255)                                                                                     tGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                          <  4:491158/150‑1 (MQ=255)                                                                                         tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                                      <  1:194482/150‑1 (MQ=255)                                                                                          gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGtt                                                     <  4:580086/150‑1 (MQ=255)                                                                                           tAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                    >  2:40768/1‑150 (MQ=255)                                                                                            aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                   <  4:139854/150‑1 (MQ=255)                                                                                            aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                   >  2:454957/1‑150 (MQ=255)                                                                                             gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGa                                                                            >  4:394961/1‑124 (MQ=255)                                                                                             gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGa                                                                            <  3:394961/124‑1 (MQ=255)                                                                                             gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCa                                                  >  2:14862/1‑150 (MQ=255)                                                                                              aaaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCaa                                                 >  1:567791/1‑150 (MQ=255)                                                                                                 tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                                              <  1:994576/150‑1 (MQ=255)                                                                                                 tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTCCACGTCAAGTCGTTCTCACAGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                                              <  1:1067831/150‑1 (MQ=255)                                                                                                    cGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                                        <  2:897417/137‑1 (MQ=255)                                                                                                    cGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                                        >  1:897417/1‑137 (MQ=255)                                                                                                          aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                     >  4:614757/1‑150 (MQ=255)                                                                                                          aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                     >  3:572765/1‑150 (MQ=255)                                                                                                             aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACg                                  >  1:307721/1‑150 (MQ=255)                                                                                                               gTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTg                                                                                                        >  3:246410/1‑78 (MQ=255)                                                                                                               gTTATCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTg                                                                                                        <  4:246410/78‑1 (MQ=255)                                                                                                                  aTCTTCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                             >  2:1016279/1‑150 (MQ=255)                                                                                                                     ttCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                          <  1:613609/150‑1 (MQ=255)                                                                                                                      tCGCCCAGAAGATGATGATCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACg                         <  3:475594/150‑1 (MQ=255)                                                                                                                         cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg                            <  4:68195/144‑1 (MQ=255)                                                                                                                         cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg                            >  3:68195/1‑144 (MQ=255)                                                                                                                                   gatgatCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc            <  1:14862/150‑1 (MQ=255)                                                                                                                                   gatgatCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc            <  2:901371/150‑1 (MQ=255)                                                                                                                                    atgatCGAAAAATTTATCCGTGCACGTAAAGTCGGTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt           <  1:782076/150‑1 (MQ=255)                                                                                                                                     tgatCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCa                                                                    >  1:499450/1‑92 (MQ=255)                                                                                                                                     tgatCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCa                                                                    <  2:499450/92‑1 (MQ=255)                                                                                                                                       atCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt        <  3:186120/150‑1 (MQ=255)                                                                                                                                        tCGAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg       <  1:1016279/150‑1 (MQ=255)                                                                                                                                          gAAAAATTTATTCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa     >  2:268583/1‑150 (MQ=255)                                                                                                                                          gAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                     <  3:522716/134‑1 (MQ=255)                                                                                                                                          gAAAAATTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                     >  4:522716/1‑134 (MQ=255)                                                                                                                                           aaaaaTTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                   <  2:786180/103‑1 (MQ=255)                                                                                                                                           aaaaaTTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                   >  1:786180/1‑103 (MQ=255)                                                                                                                                           aaaaaTTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAAt    >  2:172912/1‑150 (MQ=255)                                                                                                                                           aaaaaTTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAAt    <  1:301238/150‑1 (MQ=255)                                                                                                                                           aaaaaTTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAAt    <  1:197509/150‑1 (MQ=255)                                                                                                                                             aaaTTTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATcc  <  2:567791/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCC  >  NC_000913_3_bme_pgi/1756353‑1756642

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