Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 2,122,777 (T)5→4 coding (4/1413 nt) pykF ← pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,122,7730T.98.2% 1712.1 / 29.5 387coding (8/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base . (182/194);  total (185/202)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 403 total aligned reads displayed.

ACGTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCAGAGGA  >  DH1/2122623‑2122917
                                                                                                                                                      |                                                                                                                                                
ccGTGGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATACTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCt                                                                                                                                                  <  6:315618/150‑1 (MQ=255)
  gTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCttt                                                                                                                                                <  5:228842/151‑1 (MQ=255)
  gTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCttt                                                                                                                                                >  6:497009/1‑151 (MQ=255)
      cgcAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTTAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCttttcat                                                                                                                                            <  3:65369/151‑4 (MQ=255)
         aGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGAc                                                                                                                                        >  3:1482207/1‑151 (MQ=255)
         aGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGAc                                                                                                                                        <  4:1030891/151‑1 (MQ=255)
         aGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATATTAGCTAACATCTCTTCAGATTCTGTTATCGGTCCGATGGTGCAAACAATTTTGCTC‑ATTTCATGa                                                                                                                                         >  8:1685170/1‑150 (MQ=255)
         aGATTCTGAATGCGCTGACCGTGTTCTGAATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGAc                                                                                                                                        >  8:1607205/1‑151 (MQ=255)
           aTTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACag                                                                                                                                      <  1:1816560/151‑1 (MQ=255)
             tCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtc                                                                                                                                    >  8:1536508/1‑151 (MQ=255)
                   tGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtcttagtc                                                                                                                              >  1:788932/1‑151 (MQ=255)
                    gcgcTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtcttagtc                                                                                                                              <  1:948875/150‑1 (MQ=255)
                     cgcTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtcttagtctt                                                                                                                            <  5:907378/151‑1 (MQ=255)
                      gcTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCttt                                                                                                                           >  7:450243/1‑151 (MQ=255)
                       cTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCttt                                                                                                                           >  1:10279/1‑150 (MQ=255)
                         gACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTa                                                                                                                          <  5:1705805/149‑1 (MQ=255)
                         gACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTa                                                                                                                          <  4:981116/149‑1 (MQ=255)
                         gACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGGTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTaa                                                                                                                         <  1:479980/150‑1 (MQ=255)
                           ccGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGtt                                                                                                                      <  3:422667/151‑1 (MQ=255)
                            cGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATACCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTg                                                                                                                     <  3:1112201/151‑1 (MQ=255)
                               gtTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTATCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTgaga                                                                                                                  <  5:138767/151‑1 (MQ=255)
                               gtTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTgaga                                                                                                                  <  8:834962/151‑1 (MQ=255)
                               gtTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTg                                                                                                                     >  3:88232/1‑148 (MQ=255)
                                 tCTGCATAGTCACCATGAGAGAAGTGCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGaa                                                                                                                 <  6:785154/150‑1 (MQ=255)
                                    gCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAgg                                                                                                               <  1:335205/149‑1 (MQ=255)
                                    gCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGa                                                                                                              <  7:1874386/150‑1 (MQ=255)
                                      aTAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATgg                                                                                                           <  4:1243740/151‑1 (MQ=255)
                                        aGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATggg                                                                                                          <  2:488780/150‑1 (MQ=255)
                                        aGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATggg                                                                                                          <  3:1714013/150‑1 (MQ=255)
                                        aGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATggg                                                                                                          >  5:83025/1‑150 (MQ=255)
                                        aGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGga                                                                                                         <  8:550116/151‑1 (MQ=255)
                                         gTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTATCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGgag                                                                                                        <  1:1119898/151‑1 (MQ=255)
                                          tCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATggg                                                                                                          >  2:1362471/1‑148 (MQ=255)
                                           cACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGgaga                                                                                                       >  5:1868892/1‑150 (MQ=255)
                                           cACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGaa                                                                                                      <  7:715173/151‑1 (MQ=255)
                                               aTGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACtt                                                                                                  >  8:364999/1‑151 (MQ=255)
                                                   gagaAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCtt                                                                                              >  1:65669/1‑151 (MQ=255)
                                                    agaAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCttt                                                                                             <  4:444655/151‑1 (MQ=255)
                                                     gaAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTTAGTTGAGAAGGATGGGAGAAACTTGCTTTc                                                                                            <  5:758195/151‑1 (MQ=255)
                                                      aatttCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGACTTACTATTTAAGTTTCGAAGGATGGTAGAAACTTCCTTTCt                                                                                           >  6:1575658/4‑151 (MQ=255)
                                                      aaGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCt                                                                                           <  7:281049/151‑1 (MQ=255)
                                                       aGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTg                                                                                          <  4:1120192/151‑1 (MQ=255)
                                                        gTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTgg                                                                                         >  4:358183/1‑151 (MQ=255)
                                                          tCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGgc                                                                                       >  1:1599802/1‑151 (MQ=255)
                                                          tCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGgc                                                                                       >  7:637601/1‑151 (MQ=255)
                                                           cAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGgc                                                                                       <  3:805597/150‑1 (MQ=255)
                                                           cAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGgc                                                                                       <  8:1725861/150‑1 (MQ=255)
                                                            aGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGgcgc                                                                                     >  5:1187673/1‑151 (MQ=255)
                                                               cGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACt                                                                                  <  6:250530/151‑1 (MQ=255)
                                                                gCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTg                                                                                 >  4:1455559/1‑151 (MQ=255)
                                                                       gTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGaa                                                                           >  4:584745/1‑150 (MQ=255)
                                                                          cATGCCTGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTTGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCg                                                                       <  6:1773923/151‑1 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGTTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGCAGAAACTTGCTTTCTCGGCACACTGCCGGAAAGcc                                                                       >  2:1763985/1‑150 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGc                                                                        <  6:719453/150‑1 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGc                                                                        <  1:218486/150‑1 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCg                                                                       >  2:1410227/1‑151 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCg                                                                       >  5:151157/1‑151 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCg                                                                       <  3:1907842/151‑1 (MQ=255)
                                                                             gCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAAt                                                                    <  5:45119/151‑1 (MQ=255)
                                                                              ccAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATc                                                                   <  5:680508/151‑1 (MQ=255)
                                                                               cAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATcc                                                                  >  7:1592743/1‑151 (MQ=255)
                                                                                aGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCg                                                                 <  7:87592/151‑1 (MQ=255)
                                                                                  cGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGt                                                               >  6:1435700/1‑151 (MQ=255)
                                                                                   gaccAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGAGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTa                                                              <  1:1166380/149‑1 (MQ=255)
                                                                                     ccATCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTAt                                                             <  2:1906330/150‑1 (MQ=255)
                                                                                      cAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAATGATGGGAGAAACTTTCTTTCTGGGCGAACTGCCGGAAATCGAAGCCGGTAtt                                                            >  2:1829502/1‑150 (MQ=255)
                                                                                      cAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTAtt                                                            >  6:1174066/1‑150 (MQ=255)
                                                                                      cAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTAtt                                                            >  8:1769356/1‑150 (MQ=255)
                                                                                       aGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTAc                                                          <  5:368554/151‑1 (MQ=255)
                                                                                       aGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTAc                                                          >  2:938074/1‑151 (MQ=255)
                                                                                        gCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACg                                                         <  4:683611/151‑1 (MQ=255)
                                                                                            tttAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGACAGCGAATCCGGTATTACGaaa                                                      <  1:1735146/150‑1 (MQ=255)
                                                                                             ttAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgt                                                    >  5:293952/1‑151 (MQ=255)
                                                                                              tAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtg                                                   >  3:551051/1‑151 (MQ=255)
                                                                                                   aaCATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtggtgg                                               <  4:1379460/150‑1 (MQ=255)
                                                                                                       tctcTTCAGATTCGGTTTTCGGTCAGATGGTGCAAACAATTTTGGAC‑TATTCATGACAGACTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGATTTCTGGGCGCACTGCCGGAAAGCGAAACCGGTATTACGAAAGTGGTGGTGCCt                                          >  8:730777/1‑151 (MQ=255)
                                                                                                            tCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGctc                                     >  1:567378/1‑151 (MQ=255)
                                                                                                              aGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGctct                                    >  2:651981/1‑150 (MQ=255)
                                                                                                                aTTCGGTTTTCGGTCCGATGGTGCAAAGAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGCATGCGGTATTAGGAAAGTGGTGGTGGCTCGCGCTatg                                 >  4:1703235/1‑149 (MQ=255)
                                                                                                                   cGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAgg                                <  7:1269252/150‑1 (MQ=255)
                                                                                                                   cGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGAt                              <  1:1302176/151‑1 (MQ=255)
                                                                                                                    ggTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGAt                              >  1:1408246/1‑150 (MQ=255)
                                                                                                                    ggTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATa                             <  2:1531151/151‑1 (MQ=255)
                                                                                                                    ggTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATa                             <  2:792683/151‑1 (MQ=255)
                                                                                                                     gTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGAt                              >  4:255090/1‑149 (MQ=255)
                                                                                                                        ttCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGt                          <  8:1895210/150‑1 (MQ=255)
                                                                                                                          cGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGAc                       >  8:786270/1‑151 (MQ=255)
                                                                                                                           ggTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGAc                       >  8:1293718/1‑150 (MQ=255)
                                                                                                                               cGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCg                   >  7:1489645/1‑150 (MQ=255)
                                                                                                                               cGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACAGCg                   >  8:1258267/1‑150 (MQ=255)
                                                                                                                                gATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaa                 <  1:504209/151‑1 (MQ=255)
                                                                                                                                 aTGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaaa                <  3:355970/151‑1 (MQ=255)
                                                                                                                                 aTGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaa                 >  7:953572/1‑150 (MQ=255)
                                                                                                                                  tGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaaa                >  4:78181/1‑150 (MQ=255)
                                                                                                                                   ggTGCAAACAATATTGGTC‑TTTTCATGACAGGCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAAGTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaaaa               <  6:922991/150‑1 (MQ=255)
                                                                                                                                      gCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGc            <  7:116996/150‑1 (MQ=255)
                                                                                                                                        aaaCAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGt          >  5:1306304/1‑151 (MQ=255)
                                                                                                                                          aCAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGc        >  5:1928799/1‑150 (MQ=255)
                                                                                                                                           cAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGc        >  2:315963/1‑149 (MQ=255)
                                                                                                                                            aaTTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGc        >  7:575942/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
ACGTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCAGAGGA  >  DH1/2122623‑2122917

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