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.97.2% 1732.4 / 44.6 393coding (8/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base . (217/164);  total (222/171)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 417 total aligned reads displayed.

ACGTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCAGAGGAAT  >  DH1/2122623‑2122919
                                                                                                                                                      |                                                                                                                                                  
aCGTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCt                                                                                                                                                    >  1:1781620/1‑151 (MQ=255)
 cGTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCtt                                                                                                                                                   <  7:619754/151‑1 (MQ=255)
 cGTTGCGCAGATTCTGAATGCGCTGACCGTCTTCTGCATAGTCACCATGATAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAATCAATTTTGGTCtt                                                                                                                                                   <  8:1459024/151‑1 (MQ=255)
       gcAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCttttcatg                                                                                                                                             <  8:1169501/151‑5 (MQ=255)
         aGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGAc                                                                                                                                          >  5:168647/1‑151 (MQ=255)
          gATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGCCGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGAc                                                                                                                                          >  6:515615/1‑150 (MQ=255)
             tCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtc                                                                                                                                      >  7:363992/1‑151 (MQ=255)
             tCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtc                                                                                                                                      >  7:1185874/1‑151 (MQ=255)
               tGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtctt                                                                                                                                    >  5:211708/1‑151 (MQ=255)
                gAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtctt                                                                                                                                    <  3:88287/150‑1 (MQ=255)
                gAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtctt                                                                                                                                    <  4:1116712/150‑1 (MQ=255)
                gAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtctt                                                                                                                                    <  1:1505184/150‑1 (MQ=255)
                  aTGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACagtcttagt                                                                                                                                 <  8:1726133/151‑1 (MQ=255)
                       cTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCttt                                                                                                                             >  4:544948/1‑150 (MQ=255)
                         gACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCttt                                                                                                                             <  8:461681/149‑1 (MQ=255)
                         gACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTaat                                                                                                                          >  6:1392289/1‑150 (MQ=255)
                         gACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAA‑TTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGt                                                                                                                         <  8:188492/151‑1 (MQ=255)
                           ccGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGtt                                                                                                                        >  1:349966/1‑151 (MQ=255)
                            cGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTg                                                                                                                       >  6:773123/1‑151 (MQ=255)
                            cGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTg                                                                                                                       >  7:1577314/1‑151 (MQ=255)
                             gtgtTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTg                                                                                                                       >  5:477952/1‑150 (MQ=255)
                               gtTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTgaga                                                                                                                    <  4:1473606/151‑1 (MQ=255)
                               gtTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTg                                                                                                                       >  2:452039/1‑148 (MQ=255)
                                ttCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGATgaga                                                                                                                    >  5:1373893/1‑150 (MQ=255)
                                    gCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGAt                                                                                                               <  2:127366/151‑1 (MQ=255)
                                    gCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGAt                                                                                                               >  1:911713/1‑151 (MQ=255)
                                     cATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATg                                                                                                              >  5:387505/1‑151 (MQ=255)
                                      aTAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATgg                                                                                                             >  6:1278952/1‑151 (MQ=255)
                                      aTAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATgg                                                                                                             >  8:1264917/1‑151 (MQ=255)
                                      aTAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATgg                                                                                                             <  6:257943/151‑1 (MQ=255)
                                          tCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATggg                                                                                                            >  6:585834/1‑148 (MQ=255)
                                                 gagagaAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTg                                                                                                   <  6:204442/150‑1 (MQ=255)
                                                   gagaAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCtt                                                                                                >  6:1649595/1‑151 (MQ=255)
                                                        gTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTgg                                                                                           >  3:1562001/1‑151 (MQ=255)
                                                          tCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTggg                                                                                          <  3:1044102/150‑1 (MQ=255)
                                                          tCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGgc                                                                                         >  6:115560/1‑151 (MQ=255)
                                                           cAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGgcg                                                                                        <  2:236667/151‑1 (MQ=255)
                                                                gCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑GTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTg                                                                                   >  3:1612976/1‑151 (MQ=255)
                                                                  aTAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGcc                                                                                 >  7:26516/1‑151 (MQ=255)
                                                                  aTAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGcc                                                                                 >  2:539145/1‑151 (MQ=255)
                                                                    aaCGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCgg                                                                               >  4:248953/1‑151 (MQ=255)
                                                                    aaCGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCg                                                                                <  5:579138/150‑1 (MQ=255)
                                                                     aCGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCgg                                                                               >  2:1391306/1‑151 (MQ=255)
                                                                      cgatCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGa                                                                              >  7:1448147/4‑150 (MQ=255)
                                                                      cGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGa                                                                              >  5:723606/1‑150 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCg                                                                         <  4:1263207/151‑1 (MQ=255)
                                                                          cATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCg                                                                         <  7:456924/151‑1 (MQ=255)
                                                                            tGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGa                                                                        >  4:386149/1‑150 (MQ=255)
                                                                             gCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAAt                                                                      >  4:697722/1‑151 (MQ=255)
                                                                             gCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAAt                                                                      >  2:942324/1‑151 (MQ=255)
                                                                                   gTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTa                                                                <  7:679122/151‑1 (MQ=255)
                                                                                   gTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTa                                                                <  7:846967/151‑1 (MQ=255)
                                                                                   gTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTa                                                                <  5:926567/151‑1 (MQ=255)
                                                                                       aGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTAc                                                            >  1:1075856/1‑151 (MQ=255)
                                                                                         cATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACg                                                           >  3:1703289/1‑150 (MQ=255)
                                                                                          aTTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGa                                                          >  4:1621174/1‑150 (MQ=255)
                                                                                          aTTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGa                                                          >  8:1111606/1‑150 (MQ=255)
                                                                                            tttAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGaa                                                         >  5:333162/1‑150 (MQ=255)
                                                                                                  tAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtgctgg                                                 >  7:924273/1‑151 (MQ=255)
                                                                                                   aaCATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtggtggt                                                >  8:354337/1‑151 (MQ=255)
                                                                                                   aaCATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtggtggt                                                >  1:1245771/1‑151 (MQ=255)
                                                                                                    aCATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtggtggtg                                               >  5:69532/1‑151 (MQ=255)
                                                                                                    aCATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtggtggtg                                               >  7:1070549/1‑151 (MQ=255)
                                                                                                    aCATCACTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAgtggtggtg                                               <  1:750700/151‑1 (MQ=255)
                                                                                                        ctctTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTc                                           >  1:1665584/1‑151 (MQ=255)
                                                                                                          ctTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTc                                           <  5:756406/149‑1 (MQ=255)
                                                                                                            tCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGTTC‑TTTTCATGACAGTCTTAGTATTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTATGGGCGCACTGCCGGAAAGCGAATCCGCTATTACAAAAGAGCTGGTGACTCGct                                        >  8:599971/1‑150 (MQ=255)
                                                                                                            tCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGctc                                       >  7:1485576/1‑151 (MQ=255)
                                                                                                            tCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGCTGGTGCCTCGct                                        >  4:1762805/1‑150 (MQ=255)
                                                                                                             cAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGctct                                      >  6:1703812/1‑151 (MQ=255)
                                                                                                             cAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGctct                                      >  3:1378550/1‑151 (MQ=255)
                                                                                                             cAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGctct                                      <  3:1480348/151‑1 (MQ=255)
                                                                                                              aGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGgtct                                      >  5:11554/1‑150 (MQ=255)
                                                                                                               gATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTa                                     >  2:407746/1‑150 (MQ=255)
                                                                                                                    ggTTTTCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGAt                                >  5:1807842/1‑150 (MQ=255)
                                                                                                                        ttCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTg                           >  3:234347/1‑151 (MQ=255)
                                                                                                                        ttCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTCGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTg                           >  6:354952/1‑151 (MQ=255)
                                                                                                                         tCGGTCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTg                           <  2:1703114/150‑1 (MQ=255)
                                                                                                                             tCCGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATc                      >  3:727265/1‑151 (MQ=255)
                                                                                                                              ccGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATc                      <  7:1272669/150‑1 (MQ=255)
                                                                                                                              ccGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGCATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAACGCGAAACCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGANAGGTGACATCg                     >  3:350416/1‑151 (MQ=255)
                                                                                                                               cGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCg                     >  2:724784/1‑150 (MQ=255)
                                                                                                                               cGATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCg                     >  4:446416/1‑150 (MQ=255)
                                                                                                                                gATGGTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaa                   <  1:506246/151‑1 (MQ=255)
                                                                                                                                  tGGTGCAAACAATTTTGGTC‑TTTTCATAACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaaa                  >  5:1632168/1‑150 (MQ=255)
                                                                                                                                   ggTGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGaaaa                 <  4:1362345/150‑1 (MQ=255)
                                                                                                                                     tGCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGc              >  6:1615894/1‑151 (MQ=255)
                                                                                                                                      gCAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCg             <  3:111098/151‑1 (MQ=255)
                                                                                                                                       cAAACAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAAGATGGGAGAAACTTGCTATCTGGGAGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGACGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGt            >  2:1096340/1‑151 (MQ=255)
                                                                                                                                          aCAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCa         <  5:854589/151‑1 (MQ=255)
                                                                                                                                           cAATTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGc          >  3:1284755/1‑149 (MQ=255)
                                                                                                                                            aaTTTTGGTC‑TTTTCATGACAGTGTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCGCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCaga       <  4:119504/151‑1 (MQ=255)
                                                                                                                                            aaTTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGc          >  7:1703994/1‑148 (MQ=255)
                                                                                                                                            aaTTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGc          >  2:1609044/1‑148 (MQ=255)
                                                                                                                                            aaTTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCaga       <  4:1015695/151‑1 (MQ=255)
                                                                                                                                             aTTTTGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCagag      <  1:1379513/151‑1 (MQ=255)
                                                                                                                                              ttttGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAACGATAGGTGACATCGAAAAAGCGTGCAGAgg     >  5:1046338/1‑151 (MQ=255)
                                                                                                                                               tttGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCAGAgg     <  2:828094/150‑1 (MQ=255)
                                                                                                                                                ttGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCAGAGGa    <  3:605762/150‑1 (MQ=255)
                                                                                                                                                 tGGTC‑TTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCAGAGGAAt  >  5:1517552/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ACGTTGCGCAGATTCTGAATGCGCTGACCGTGTTCTGCATAGTCACCATGAGAGAAGTTCAGACGCATAACGTTCATGCCAGCGTCCAGCATTTTAGCTAACATCTCTTCAGATTCGGTTTTCGGTCCGATGGTGCAAACAATTTTGGTCTTTTTCATGACAGTCTTAGTCTTTAAGTTGAGAAGGATGGGAGAAACTTGCTTTCTGGGCGCACTGCCGGAAAGCGAATCCGGTATTACGAAAGTGGTGGTGCCTCGCTCTAAGGATAGGTGACATCGAAAAAGCGTGCAGAGGAAT  >  DH1/2122623‑2122919

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: ‑