Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,622	G→A	V309M (GTG→ATG)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,622	0	G	A	98.5%	273.9 / ‑3.5	68	V309M (GTG→ATG)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/34);  minor base T (1/0);  total (34/34)
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.53e-01

AGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGG  >  NC_000913/1756476‑1756769                                                                                                                                                   |                                                                                                                                                    agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGAt                                                                                                                                                 <  1:22366/151‑1 (MQ=255)         atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTg                                                                                                                                         >  1:417134/1‑151 (MQ=255)          tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTgg                                                                                                                                        >  1:198411/1‑151 (MQ=255)           gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTgg                                                                                                                                        <  1:507455/150‑1 (MQ=255)            atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTg                                                                                                                                      <  2:393179/151‑1 (MQ=255)              cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGaa                                                                                                                                    <  2:982072/151‑1 (MQ=255)              cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGaa                                                                                                                                    <  1:1054959/151‑1 (MQ=255)              cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGa                                                                                                                                     >  1:1113692/1‑150 (MQ=255)              cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGa                                                                                                                                     <  1:538485/150‑1 (MQ=255)               gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAAt                                                                                                                                   >  2:91370/1‑151 (MQ=255)                aaaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATc                                                                                                                                  >  1:35577/1‑151 (MQ=255)                   aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGATGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTATGGCGAATCCg                                                                                                                                >  2:902028/1‑150 (MQ=255)                    aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGc                                                                                                                               >  2:880998/1‑150 (MQ=255)                      gTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCaaa                                                                                                                            <  1:348422/151‑1 (MQ=255)                      gTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCaaa                                                                                                                            <  1:254220/151‑1 (MQ=255)                           cGTGCACGTAAAGTCGTTATCACTACGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAgg                                                                                                                         >  2:1100338/1‑149 (MQ=255)                            gTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTaa                                                                                                                      <  1:415820/151‑1 (MQ=255)                            gTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTa                                                                                                                       >  2:313873/1‑150 (MQ=255)                              gCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAAt                                                                                                                    <  2:1021835/151‑1 (MQ=255)                               cACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAAt                                                                                                                    <  2:202523/150‑1 (MQ=255)                               cACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAAt                                                                                                                    >  1:305744/1‑150 (MQ=255)                                aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                                                                                                  >  2:565948/1‑151 (MQ=255)                                aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATAc                                                                                                                  >  2:200511/1‑151 (MQ=255)                                    aaaGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGc                                                                                                              >  1:951420/1‑151 (MQ=255)                                      aGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTg                                                                                                            <  1:1022240/151‑1 (MQ=255)                                       gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTgg                                                                                                           <  2:121115/151‑1 (MQ=255)                                         cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGa                                                                                                          >  1:45385/1‑150 (MQ=255)                                          gTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGa                                                                                                          >  1:770866/1‑149 (MQ=255)                                            tATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCg                                                                                                      >  1:439332/1‑151 (MQ=255)                                            tATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGAAAATACCCGCTGGAAGCg                                                                                                      >  1:265779/1‑151 (MQ=255)                                                  tGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCt                                                                                                >  2:442901/1‑151 (MQ=255)                                                    cGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTAt                                                                                              <  1:604281/151‑1 (MQ=255)                                                                 ggATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCa                                                                                   <  1:380976/149‑1 (MQ=255)                                                                  gATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCt                                                                                <  2:429809/151‑1 (MQ=255)                                                                    ttCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCATTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGc                                                                              >  1:31912/1‑151 (MQ=255)                                                                      ccATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCg                                                                             >  2:1155342/1‑150 (MQ=255)                                                                        aTGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAAc                                                                          <  1:330117/151‑1 (MQ=255)                                                                                aaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccgac                                                                  <  2:333655/151‑1 (MQ=255)                                                                                aaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTaccgac                                                                  <  1:91370/151‑1 (MQ=255)                                                                                         cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGtgatg                                                         <  1:575101/151‑1 (MQ=255)                                                                                         cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGtgatg                                                         <  2:886665/151‑1 (MQ=255)                                                                                            gACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGaa                                                       >  2:512913/1‑150 (MQ=255)                                                                                            gACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAAc                                                      <  1:628818/151‑1 (MQ=255)                                                                                                cgtgcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGcc                                                  <  2:300813/148‑1 (MQ=255)                                                                                                cgcgCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGcc                                                  >  2:967693/1‑151 (MQ=255)                                                                                                 gcgcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCg                                                 >  1:980976/1‑151 (MQ=255)                                                                                                  cgcAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGt                                                <  1:588177/151‑1 (MQ=255)                                                                                                    cAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGtct                                              <  1:1047033/151‑1 (MQ=255)                                                                                                       aaGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCg                                            >  2:1161884/1‑150 (MQ=255)                                                                                                        aGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCg                                            >  1:950701/1‑149 (MQ=255)                                                                                                            ggTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTca                                      <  2:160118/151‑1 (MQ=255)                                                                                                             gTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTca                                      >  2:1125641/1‑150 (MQ=255)                                                                                                                     gCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACa                             <  1:838408/151‑1 (MQ=255)                                                                                                                     gCAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACa                             <  1:712725/151‑1 (MQ=255)                                                                                                                      cAAACGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACa                             >  2:1112587/1‑150 (MQ=255)                                                                                                                         aCGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCg                         <  1:189167/151‑1 (MQ=255)                                                                                                                          cGCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGt                        <  2:462639/151‑1 (MQ=255)                                                                                                                           gCCATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTa                       <  1:37501/151‑1 (MQ=255)                                                                                                                            ccATCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTa                       >  2:1136716/1‑150 (MQ=255)                                                                                                                              aTCCTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAAc                    >  2:70921/1‑151 (MQ=255)                                                                                                                                 cTCGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTgc                 >  1:898709/1‑151 (MQ=255)                                                                                                                                   cGACGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTgcgc               >  2:307340/1‑151 (MQ=255)                                                                                                                                      cGGTACTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCAtt            <  1:704300/151‑1 (MQ=255)                                                                                                                                           cTGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCg        >  1:596644/1‑150 (MQ=255)                                                                                                                                            tGACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGa       >  2:26245/1‑150 (MQ=255)                                                                                                                                             gACGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGaa      <  2:152874/150‑1 (MQ=255)                                                                                                                                              aCGCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGc    >  1:1044485/1‑151 (MQ=255)                                                                                                                                                gCAATGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCgg  <  2:708763/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAGGTAAATACCCGCTGGAAGCGGTTTCTATCATGGCGACCATCTGCGAACGTACCGACCGCGTGATGAACAGCCGTCTCGAGTTCAACAATGACAACCGTAAACTGCGCATTACCGAAGCGG  >  NC_000913/1756476‑1756769

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