Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,755,8710TA42.1% 109.1 / 52.3 76G58G (GGT→GGA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (44/0);  new base A (0/32);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 3.71e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTT  >  NC_000913/1755721‑1756015
                                                                                                                                                           |                                                                                                                                                
cACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑t                                                                                                                                                  <  2:594643/151‑1 (MQ=255)
   cATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑Taa                                                                                                                                                >  1:1648240/1‑150 (MQ=255)
     tCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAAc                                                                                                                                             >  2:54875/1‑151 (MQ=255)
       ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAAccg                                                                                                                                           >  2:766031/1‑151 (MQ=255)
             aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCt                                                                                                                                      >  2:1541889/1‑150 (MQ=255)
             aaaaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCt                                                                                                                                      >  1:749735/1‑150 (MQ=255)
               aaCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATc                                                                                                                                   >  1:1673365/1‑151 (MQ=255)
                aCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATcc                                                                                                                                  >  1:200145/1‑151 (MQ=255)
                 ccGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCt                                                                                                                                 >  1:21557/1‑151 (MQ=255)
                   gAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGc                                                                                                                               >  1:887329/1‑151 (MQ=255)
                          aaGAGTTGTTAGAAAAAATGCGGGACGCGGGCATGACCTTTTGGCGTTTAAACTTCTCTCATGGTG‑ACTATGCAAAACACGGTCAGCGAATTAAAAATCTGCGAAA‑GGTGAT‑GAGAAAAAT‑TGA‑AAAAACCGCCGCTATCCTGCTTGATAc                                                                                                                        <  2:525494/151‑1 (MQ=255)
                          aaGAGATTTTAGAAAAATTGCTGGCCGCTGGCATGAACTTTATGCGTCTAAACTTCTCTCTTGGTG‑ACTATGAAAAACCGGGTCAGCGCATTAAGAATCTGCGAAA‑CGTGAT‑GAAAAAATC‑TGG‑AAAAACCGCCGCTATCCTGCTTGAt                                                                                                                          <  2:1095949/149‑1 (MQ=255)
                          aaGAGATGTTAGCAAAAATGCGGGCCGCGGGCATGAACGTTATGCTTTTGACCTTCTCTCAGGGTG‑ACTATGAAAAACCGGGTCAGCGCATTCAGAATCTGCGCAA‑CGTAAT‑GAGAAAAAC‑TGG‑AAAAACCGCCGCTATCCTGCTTGATa                                                                                                                         <  2:1705971/150‑1 (MQ=255)
                           agagATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATAcc                                                                                                                       >  2:1167932/1‑151 (MQ=255)
                             agaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                      >  2:93912/1‑150 (MQ=255)
                              gaTGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCaa                                                                                                                     >  2:922136/1‑150 (MQ=255)
                                tgttaaaaaaaaTGCGGGCCGCGGGAATGACCTTTATGCGTTTAACCTTCTCTATTGGTG‑ACTATGAAAAACGGGGCCAGCGCATTAAAAACCGGCGAAA‑GGTGAT‑GAAAAAAAT‑GGG‑AAAAACCGCCGCTATCCTGCTTGATACCaa                                                                                                                     <  1:1616764/140‑1 (MQ=255)
                                  ttAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGACCTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTAAGATCCTGCGCAA‑CGTGAT‑GAAAAAAAC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCa                                                                                                                      <  1:833564/145‑1 (MQ=255)
                                  ttAGCTAAAATGCGGGACGCGGGCATGAACTTTTTGCGTCTGACCTTCTCTCATGGTG‑ATTATGCAAAACACGGTCAGCACTTTCAGATTCGGCAAAA‑GGTGAT‑GAAAAAAGC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATAcc                                                                                                                       <  1:745849/144‑1 (MQ=255)
                                  tgaccaaaaaTGTTGGGCGCTGCCGGGAACGTTCTGCGTCTAACTTTCTCTAAGGGTG‑ACTAGCAAGAACACGGCCCCCTCATACAGAATC‑GCGAAATCGTG‑T‑GAAAAAAACTTGA‑AAAAACCGCCGCTATCCTGCTTGATAcc                                                                                                                       <  2:76858/138‑1 (MQ=255)
                                   tAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATcc                                                                                                               >  2:1038086/1‑151 (MQ=255)
                                    agtaaaagTGCTGGACGCTGATAAGAACGTTAGACTTCTTACCTTCTCTTTGGGTG‑ATAATGAAAAACGCGGTCAGCCAATTCAAAATCGACGCAT‑GGTGAA‑GA‑CAAAACATGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAga                                                                                                                  <  2:1272515/143‑2 (MQ=255)
                                    agaaaaaaTGCGGGCCGCTGGAATGACCTTTATGCTTCTAACCTTCTCTCA‑GGTGAACTATGAAAAACCCGGCCAGCGCATTAAAATTCTGCAAAA‑GGTAAA‑GAGAAAAAC‑GGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATcc                                                                                                               <  1:968154/146‑1 (MQ=255)
                                     gcaaaaaGGCTGGACGGGGGAATGAACTTTTTGCGTCTGATCTTCCCTCGTGGTG‑ATTAGGCAGACCGCGGCCAGCAAATTCAAATTCTGCGCAA‑CGTGAT‑GACAAAAAT‑TGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCg                                                                                                              <  2:236313/147‑1 (MQ=255)
                                     gaaaaaaGGCGGGCCGCTGGCATAACCGTTATGCGTCTGACCTTCTCTCAGGGTG‑ACTATGCAGAACACGGGCAGCCCATTCAGATTCGGCGCAA‑CGTGA‑‑GAAAAAAAC‑CGGAAAAAACCGCCGCTATCCTGCTTGATACCAAAGATc                                                                                                                <  1:707382/145‑1 (MQ=255)
                                        aaaaTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGaa                                                                                                           >  2:66840/1‑150 (MQ=255)
                                           aTGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATcc                                                                                                       >  1:87357/1‑151 (MQ=255)
                                            tGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCg                                                                                                      >  1:1480016/1‑151 (MQ=255)
                                             gcgggACGCTGGCATGATCTTTATGCGTCTGACCTTCTCTCGTGGTG‑ACTGTGCAAACCCCGGACAGCGCATCAAAAAATTGCGGAA‑CG‑GATGGAGAAAATG‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATcc                                                                                                       <  2:185049/146‑1 (MQ=255)
                                                ggACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCAcc                                                                                                  >  2:575926/1‑151 (MQ=255)
                                                    gCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATg                                                                                               >  2:617035/1‑150 (MQ=255)
                                                      tGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGaa                                                                                             >  2:1311265/1‑150 (MQ=255)
                                                      tGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGaa                                                                                             >  1:263014/1‑150 (MQ=255)
                                                          aTGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTgg                                                                                        >  2:1477775/1‑151 (MQ=255)
                                                             aaCTTTATGCGTCTAACCTTCCCTCATGTGA‑ACTATGAAGACCAGGGCCGGCGCATTCAGGTTCTGCGCAC‑GTTAAT‑GAAAAAAAC‑TGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGaa                                                                                      <  1:704374/150‑1 (MQ=255)
                                                               ctttACGCGTTTTATCTTCTCTATGGGTA‑ACTATAAAGAACACGGTCAGCCCATCAAAAACCTGCGCAT‑GGTGAG‑GAAAAAAGT‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAgg                                                                                    <  1:766784/148‑1 (MQ=255)
                                                                tttATGCGTCTAACCTCCTCTCAGGGTG‑ACTATGCAGACCACGGCCAGCGCATTAAGAATCTGCGCAA‑CGTGAG‑GAGAAAAAC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGc                                                                                   <  2:1586259/149‑1 (MQ=255)
                                                                 ttATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCg                                                                                  >  1:1628621/1‑150 (MQ=255)
                                                                        tCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAacgac                                                                          >  2:181896/1‑151 (MQ=255)
                                                                          tGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGt                                                                        >  2:793803/1‑151 (MQ=255)
                                                                            aaCTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGttt                                                                      >  1:303285/1‑151 (MQ=255)
                                                                                tctctcATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTctct                                                                  >  1:1342991/1‑151 (MQ=255)
                                                                                     cATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaa                                                               >  1:180879/1‑149 (MQ=255)
                                                                                      aTGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                            >  1:758091/1‑151 (MQ=255)
                                                                                       tGGTG‑ATTATGAAAAACGGGGCCAGCGCATACAAAATCGGCGCAA‑GGTGAG‑AAAAAAAAC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAGATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGaaa                                                              <  1:677434/148‑1 (MQ=255)
                                                                                       tGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCt                                                           >  1:1578736/1‑151 (MQ=255)
                                                                                        gggg‑ACAATGCAGACCCCGGACAGCCCATTCAGAGTCTGCGAAG‑AGTGAA‑GAACAAACC‑ATA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGc                                                            <  1:89862/146‑1 (MQ=255)
                                                                                        ggTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTg                                                          >  2:918090/1‑151 (MQ=255)
                                                                                               tatgaaaaccccGGTCAGCGCATTCAAAATCGGCGAAA‑CGTGAG‑GAGAAAAAG‑TGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGt                                                        <  2:1577076/136‑1 (MQ=255)
                                                                                                atcaaGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAAAAAAAC‑TGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGa                                                    <  2:1205937/146‑1 (MQ=255)
                                                                                                aTGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAc                                                   >  1:1435520/1‑151 (MQ=255)
                                                                                                    aGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCttt                                               >  2:625885/1‑151 (MQ=255)
                                                                                                      aaCACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAACC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGCCCTTTAc                                             >  2:240845/1‑151 (MQ=255)
                                                                                                      aaCACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTAc                                             >  1:1415611/1‑151 (MQ=255)
                                                                                                      aaCACGGTCAGCGCATAAAAAATCGGCGGAG‑CGGGAT‑GAAAAAAAT‑TAG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTa                                              <  1:1403524/150‑1 (MQ=255)
                                                                                                          cGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                         >  2:259462/1‑151 (MQ=255)
                                                                                                           ggacAGCGCATCCAGAGTCTCCGAAA‑CGTGAG‑GAGAAAAAC‑TGT‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGAct                                                  <  1:765968/138‑2 (MQ=255)
                                                                                                           ggTCAGCGCATTAAAATTCTGCGCAA‑GGGAAA‑GAAAAAACC‑GGG‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTc                                         <  2:1697917/150‑1 (MQ=255)
                                                                                                             tCAGCGCATTAAGAATATACCAAA‑GGAAAC‑GAAAAAAAC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACt                                            <  1:1541793/145‑1 (MQ=255)
                                                                                                              ccgcgcATCCAGAATCGGCGAAA‑CGTAAG‑AAAAAAAAC‑TGA‑AAAAACCGCCGCTATCTTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcacc                                      <  2:1497068/148‑1 (MQ=255)
                                                                                                               aGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                    >  2:1702511/1‑151 (MQ=255)
                                                                                                                 cgtaTTCAGAATCTGCACAA‑CGTGAG‑GAGAAAAAG‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTcaccac                                    <  2:1280484/146‑1 (MQ=255)
                                                                                                                   cATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                >  2:723477/1‑151 (MQ=255)
                                                                                                                    aTTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGAt                                >  2:825386/1‑150 (MQ=255)
                                                                                                                             cggccaaa‑CGTGAA‑GAAAAAAAT‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGtt                       <  1:1167845/144‑1 (MQ=255)
                                                                                                                                 gcAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCgg                  >  2:275522/1‑151 (MQ=255)
                                                                                                                                      agggaaagagaAAAAT‑TGT‑AAAAGCCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCTCTGATAAATCTGTTATCg                   <  1:1349262/140‑1 (MQ=255)
                                                                                                                                         gat‑gaGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGc           >  2:1560184/1‑151 (MQ=255)
                                                                                                                                          at‑gaGAAAAAA‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACa             <  2:109689/148‑1 (MQ=255)
                                                                                                                                              agaaaaaC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAg            <  1:1158752/143‑1 (MQ=255)
                                                                                                                                              aacAAAAC‑TTG‑AAACACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGa         <  1:271921/147‑1 (MQ=255)
                                                                                                                                               gCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAAt      >  2:1058274/1‑151 (MQ=255)
                                                                                                                                                 aaaaC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATgg    >  2:1198718/1‑151 (MQ=255)
                                                                                                                                                  aaaC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGt   <  1:470436/151‑1 (MQ=255)
                                                                                                                                                   aaC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt  >  1:566087/1‑151 (MQ=255)
                                                                                                                                                   aaC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt  >  2:799506/1‑151 (MQ=255)
                                                                                                                                                   aaC‑TGA‑AAAAACCGCCGCTATCCTGCTTGATACCAAAGATCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGtt  <  2:21557/151‑1 (MQ=255)
                                                                                                                                                           |                                                                                                                                                
CACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTG‑ACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAA‑CGTGAT‑GAGCAAAAC‑TGG‑TAAAACCGCCGCTATCCTGCTTGATACCAAAGGTCCGGAAATCCGCACCATGAAACTGGAAGGCGGTAACGACGTTTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTT  >  NC_000913/1755721‑1756015

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