Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,755,903	G→A	G69D (GGT→GAT)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,903	0	G	A	98.8%	307.6 / ‑4.2	80	G69D (GGT→GAT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (40/39);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTC  >  NC_000913/1755756‑1756050                                                                                                                                                       |                                                                                                                                                    tAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATcc                                                                                                                                                  >  2:1038086/1‑151 (MQ=255)  agaaaaaaTGCGGGCCGCTGGAATGACCTTTATGCTTCTAACCTTCTCTCA‑GGTGAACTATGAAAAACCCGGCCAGCGCATTAAAATTCTGCAAAAGGTAAA‑GAGAAAAACGGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATcc                                                                                                                                                  <  1:968154/146‑1 (MQ=255)   gcaaaaaGGCTGGACGGGGGAATGAACTTTTTGCGTCTGATCTTCCCTCGTGGTG‑ATTAGGCAGACCGCGGCCAGCAAATTCAAATTCTGCGCAACGTGAT‑GACAAAAATTGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCg                                                                                                                                                 <  2:236313/147‑1 (MQ=255)   gaaaaaaGGCGGGCCGCTGGCATAACCGTTATGCGTCTGACCTTCTCTCAGGGTG‑ACTATGCAGAACACGGGCAGCCCATTCAGATTCGGCGCAACGTGA‑‑GAAAAAAACCGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGATc                                                                                                                                                   <  1:707382/145‑1 (MQ=255)      aaaaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGaa                                                                                                                                              >  2:66840/1‑150 (MQ=255)         aTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATcc                                                                                                                                          >  1:87357/1‑151 (MQ=255)          tGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCg                                                                                                                                         >  1:1480016/1‑151 (MQ=255)           gcgggACGCTGGCATGATCTTTATGCGTCTGACCTTCTCTCGTGGTG‑ACTGTGCAAACCCCGGACAGCGCATCAAAAAATTGCGGAACG‑GATGGAGAAAATGTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATcc                                                                                                                                          <  2:185049/146‑1 (MQ=255)              ggACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCAcc                                                                                                                                     >  2:575926/1‑151 (MQ=255)                  gCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATg                                                                                                                                  >  2:617035/1‑150 (MQ=255)                    tGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGaa                                                                                                                                >  2:1311265/1‑150 (MQ=255)                    tGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGaa                                                                                                                                >  1:263014/1‑150 (MQ=255)                        aTGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTgg                                                                                                                           >  2:1477775/1‑151 (MQ=255)                           aaCTTTATGCGTCTAACCTTCCCTCATGTGA‑ACTATGAAGACCAGGGCCGGCGCATTCAGGTTCTGCGCACGTTAAT‑GAAAAAAACTGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGaa                                                                                                                         <  1:704374/150‑1 (MQ=255)                             ctttACGCGTTTTATCTTCTCTATGGGTA‑ACTATAAAGAACACGGTCAGCCCATCAAAAACCTGCGCATGGTGAG‑GAAAAAAGTTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAgg                                                                                                                       <  1:766784/148‑1 (MQ=255)                              tttATGCGTCTAACCTCCTCTCAGGGTG‑ACTATGCAGACCACGGCCAGCGCATTAAGAATCTGCGCAACGTGAG‑GAGAAAAACTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                                                      <  2:1586259/149‑1 (MQ=255)                               ttATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                                                     >  1:1628621/1‑150 (MQ=255)                                      tCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                                                             >  2:181896/1‑151 (MQ=255)                                        tGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGt                                                                                                           >  2:793803/1‑151 (MQ=255)                                          aaCTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                                                         >  1:303285/1‑151 (MQ=255)                                              tctctcATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                                                     >  1:1342991/1‑151 (MQ=255)                                                   cATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaa                                                                                                  >  1:180879/1‑149 (MQ=255)                                                    aTGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                                                               >  1:758091/1‑151 (MQ=255)                                                     tGGTG‑ATTATGAAAAACGGGGCCAGCGCATACAAAATCGGCGCAAGGTGAG‑AAAAAAAACTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAGATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                                                                 <  1:677434/148‑1 (MQ=255)                                                     tGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCt                                                                                              >  1:1578736/1‑151 (MQ=255)                                                      gggg‑ACAATGCAGACCCCGGACAGCCCATTCAGAGTCTGCGAAGAGTGAA‑GAACAAACCATA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                                                               <  1:89862/146‑1 (MQ=255)                                                      ggTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                                                             >  2:918090/1‑151 (MQ=255)                                                             tatgaaaaccccGGTCAGCGCATTCAAAATCGGCGAAACGTGAG‑GAGAAAAAGTGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                                                           <  2:1577076/136‑1 (MQ=255)                                                              atcaaGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAAAAAAACTGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGa                                                                                       <  2:1205937/146‑1 (MQ=255)                                                              aTGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAc                                                                                      >  1:1435520/1‑151 (MQ=255)                                                                  aGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                                                                  >  2:625885/1‑151 (MQ=255)                                                                    aaCACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAACCTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGCCCTTTAc                                                                                >  2:240845/1‑151 (MQ=255)                                                                    aaCACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                                                                >  1:1415611/1‑151 (MQ=255)                                                                    aaCACGGTCAGCGCATAAAAAATCGGCGGAGCGGGAT‑GAAAAAAATTAG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTa                                                                                 <  1:1403524/150‑1 (MQ=255)                                                                        cGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                                                            >  2:259462/1‑151 (MQ=255)                                                                         ggacAGCGCATCCAGAGTCTCCGAAACGTGAG‑GAGAAAAACTGT‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAct                                                                                     <  1:765968/138‑2 (MQ=255)                                                                         ggTCAGCGCATTAAAATTCTGCGCAAGGGAAA‑GAAAAAACCGGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                                                            <  2:1697917/150‑1 (MQ=255)                                                                           tCAGCGCATTAAGAATATACCAAAGGAAAC‑GAAAAAAACTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACt                                                                               <  1:1541793/145‑1 (MQ=255)                                                                            ccgcgcATCCAGAATCGGCGAAACGTAAG‑AAAAAAAACTGA‑AAAAACCGCCGCTATCTTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcacc                                                                         <  2:1497068/148‑1 (MQ=255)                                                                             aGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                                                       >  2:1702511/1‑151 (MQ=255)                                                                               cgtaTTCAGAATCTGCACAACGTGAG‑GAGAAAAAGTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                                                       <  2:1280484/146‑1 (MQ=255)                                                                                 cATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                                                   >  2:723477/1‑151 (MQ=255)                                                                                  aTTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                                                   >  2:825386/1‑150 (MQ=255)                                                                                           cggccaaaCGTGAA‑GAAAAAAATTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                                                          <  1:1167845/144‑1 (MQ=255)                                                                                               gcAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCgg                                                     >  2:275522/1‑151 (MQ=255)                                                                                                   agggaaagagaAAAATTGT‑AAAAGCCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCTCTGATAAATCTGTTATCg                                                      <  1:1349262/140‑1 (MQ=255)                                                                                                      gat‑gaGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGc                                              >  2:1560184/1‑151 (MQ=255)                                                                                                       at‑gaGAAAAAATGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACa                                                <  2:109689/148‑1 (MQ=255)                                                                                                           agaaaaaCTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAg                                               <  1:1158752/143‑1 (MQ=255)                                                                                                           aacAAAACTTG‑AAACACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGa                                            <  1:271921/147‑1 (MQ=255)                                                                                                            gCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAAt                                         >  2:1058274/1‑151 (MQ=255)                                                                                                              aaaaCTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg                                       >  2:1198718/1‑151 (MQ=255)                                                                                                               aaaCTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGt                                      <  1:470436/151‑1 (MQ=255)                                                                                                                aaCTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt                                     >  2:799506/1‑151 (MQ=255)                                                                                                                aaCTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt                                     >  1:566087/1‑151 (MQ=255)                                                                                                                aaCTGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt                                     <  2:21557/151‑1 (MQ=255)                                                                                                                      aaaaagCCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCgg                                 <  2:887397/148‑1 (MQ=255)                                                                                                                       aaaaaCCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTa                               <  1:535243/150‑1 (MQ=255)                                                                                                                        aaaaCCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCgg                                 <  1:93906/148‑1 (MQ=255)                                                                                                                        aaaaCCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCg                                  <  1:274576/147‑1 (MQ=255)                                                                                                                         aaaCCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTaa                              >  2:437145/1‑150 (MQ=255)                                                                                                                              gccgcTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATg                        >  2:1210690/1‑151 (MQ=255)                                                                                                                                cgcTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                      <  2:199868/151‑1 (MQ=255)                                                                                                                                cgcTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                      <  1:66837/151‑1 (MQ=255)                                                                                                                                  cTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAg                     <  1:55607/150‑1 (MQ=255)                                                                                                                                    aTCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAgg                    <  2:1584088/149‑1 (MQ=255)                                                                                                                                        tGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCac              >  2:1692251/1‑151 (MQ=255)                                                                                                                                         gCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCact             >  1:812481/1‑151 (MQ=255)                                                                                                                                         gCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCact             >  2:50915/1‑151 (MQ=255)                                                                                                                                         gCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCact             <  1:379042/151‑1 (MQ=255)                                                                                                                                           ttGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCactac           >  1:625033/1‑151 (MQ=255)                                                                                                                                            tGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCactac           <  1:1038008/150‑1 (MQ=255)                                                                                                                                             gATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTg         <  1:791816/151‑1 (MQ=255)                                                                                                                                             gATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTg         >  2:1104345/1‑151 (MQ=255)                                                                                                                                              aTACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTg         >  2:526619/1‑150 (MQ=255)                                                                                                                                                aCCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGAc       <  1:1477677/150‑1 (MQ=255)                                                                                                                                                 ccAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACct     <  1:1487069/151‑1 (MQ=255)                                                                                                                                                 ccAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACct     >  2:822691/1‑151 (MQ=255)                                                                                                                                                   aaaGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACctgt   <  1:155936/151‑1 (MQ=255)                                                                                                                                                    aaGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACctgtc  <  2:1505133/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGAT‑GAGCAAAACTGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTC  >  NC_000913/1755756‑1756050

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