Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 3,913,413 T→G intergenic (+79/‑175) yjjA → / → opgB hypothetical protein/phosphoglycerol transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,913,4130TG98.6% 2251.0 / 9.3 580intergenic (+79/‑175)yjjA/opgBhypothetical protein/phosphoglycerol transferase I
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base G (264/308);  total (268/312)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

Only 100 of 582 total aligned reads displayed.

GCTGGCAGAAAAGGTGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTTACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTGGCTA  >  DH1/3913265‑3913559
                                                                                                                                                    |                                                                                                                                                  
     cAGAAAAGGTGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTg                                                                                                                                             <  7:171540/151‑1 (MQ=255)
      aGAAAAGGTGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGt                                                                                                                                            <  4:939986/151‑1 (MQ=255)
       gAAAAGGTGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGtt                                                                                                                                           <  4:1778091/151‑1 (MQ=255)
        aaaaGGTGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGtt                                                                                                                                           >  1:1765595/1‑150 (MQ=255)
           aGGTGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTacac                                                                                                                                       <  3:425183/151‑1 (MQ=255)
              tGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACggg                                                                                                                                    >  3:622779/1‑151 (MQ=255)
                 aaaCCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGCCGGTTTGACCCCTGATACACGGGa                                                                                                                                   >  2:1580084/1‑149 (MQ=255)
                   aCCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTg                                                                                                                               >  6:493726/1‑151 (MQ=255)
                    acAAAGCCTGCGATCTGGTGTTAAAACAGGGCCTGCACTTCATTTCCAGAACCGCGAAACGCCCCGTGATAGCTGTTACCCGGCTGACGCGCCTCTCTACTGCCCTATTACCTGCTCCAGGAACGTTTGCCCACTCTTACACGGCAAGtt                                                                                                                               >  3:1333922/2‑149 (MQ=255)
                      aaaGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGAcc                                                                                                                            <  4:1575073/151‑1 (MQ=255)
                        aGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGt                                                                                                                          >  1:279226/1‑151 (MQ=255)
                         gCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGt                                                                                                                          >  1:1513435/1‑150 (MQ=255)
                         gCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGt                                                                                                                          >  1:137898/1‑150 (MQ=255)
                          ccTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTaa                                                                                                                        <  8:1069346/151‑1 (MQ=255)
                            tGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATc                                                                                                                      >  7:547607/1‑151 (MQ=255)
                             gCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAAt                                                                                                                       <  1:502078/149‑1 (MQ=255)
                              cGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCa                                                                                                                     >  5:1187068/1‑150 (MQ=255)
                                  cTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCaaaaa                                                                                                                 >  3:1616814/1‑150 (MQ=255)
                                    ggTGTTAAAACAGGGGCTGAAGTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGATCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGt                                                                                                              <  8:557555/151‑1 (MQ=255)
                                    ggTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCaaa                                                                                                                   <  4:78550/146‑1 (MQ=255)
                                    ggTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCaaa                                                                                                                   <  7:1358673/146‑1 (MQ=255)
                                    ggTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCaaa                                                                                                                   <  7:384735/146‑1 (MQ=255)
                                    ggTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTCGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGt                                                                                                              <  5:518418/151‑1 (MQ=255)
                                    ggTGTTAAAACAGGGGCTGAACTTCATCTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGt                                                                                                              <  6:1240289/151‑1 (MQ=255)
                                       gtTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGTGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCt                                                                                                            <  5:1900545/150‑1 (MQ=255)
                                       gtTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCt                                                                                                            <  4:580570/150‑1 (MQ=255)
                                       gtTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTg                                                                                                           <  5:1498158/151‑1 (MQ=255)
                                       gtTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTg                                                                                                           <  1:1947730/151‑1 (MQ=255)
                                         tAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGtt                                                                                                         >  8:460401/1‑151 (MQ=255)
                                          aaaaCAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGttt                                                                                                        >  8:780242/1‑151 (MQ=255)
                                           aaaCAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGttt                                                                                                        >  3:1476892/1‑150 (MQ=255)
                                           aaaCAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTta                                                                                                       <  3:1256122/151‑1 (MQ=255)
                                            aaCAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTtat                                                                                                      <  3:1132176/151‑1 (MQ=255)
                                            aaCAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTtat                                                                                                      <  7:1164445/151‑1 (MQ=255)
                                                ggggCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAAc                                                                                                  <  4:1169936/151‑1 (MQ=255)
                                                    cTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAAt                                                                                              >  8:1843458/1‑151 (MQ=255)
                                                      gAACTTCATTTCCTGAAGCTCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGATTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTc                                                                                            <  2:311391/151‑1 (MQ=255)
                                                      gAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAAt                                                                                              <  6:783971/149‑1 (MQ=255)
                                                           tCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTaac                                                                                       >  8:1298201/1‑151 (MQ=255)
                                                           tCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAGTCTGTTTATAAACTAATTCCTaac                                                                                       <  6:1632413/150‑1 (MQ=255)
                                                            cATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTaaca                                                                                      <  5:316732/151‑1 (MQ=255)
                                                            cATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTaac                                                                                       >  8:1159830/1‑150 (MQ=255)
                                                             aTTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTaacaa                                                                                     <  5:1318717/151‑1 (MQ=255)
                                                             aTTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTaaca                                                                                      >  1:63580/1‑150 (MQ=255)
                                                               ttCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTACCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTaacaac                                                                                    <  1:1583532/150‑1 (MQ=255)
                                                                 ccTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACAt                                                                                  <  8:1366690/150‑1 (MQ=255)
                                                                   tGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGc                                                                               >  7:1626925/1‑151 (MQ=255)
                                                                     aaCCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCAt                                                                             >  7:1494765/1‑151 (MQ=255)
                                                                     aaCCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCGTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCAt                                                                             >  4:1060090/1‑151 (MQ=255)
                                                                       ccGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCAttt                                                                           <  7:72872/151‑1 (MQ=255)
                                                                        gtcgAAACGCCGCGTTTTAGCTGTTACGCTGCTGACGCGTCTCTCTTCTGACCTTTTCCCTGCTACTGGACGTTTTGACTCATGTTACACGTGAAGTGACCGTAATCAAAATAGTCTTTTTATAAACTAATTCCTAACAACATCGCAtttt                                                                          <  6:1929785/149‑1 (MQ=255)
                                                                         gcgAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGGCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCAttttt                                                                         <  6:911528/151‑1 (MQ=255)
                                                                            aaaCGGCGCGGTGTAGCTGTTACCCGGCTGACAGGGCTCGCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGAGCCCTGTTACACGGGAAGTGACCGTAATCAAAAACGTCTGTTTATAAACTAATTCCTAACAAAATCGCATTTTTTTg                                                                      <  1:1082100/151‑1 (MQ=255)
                                                                             aaCGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACACCATCGCATTTTTTTg                                                                      >  8:1663222/1‑150 (MQ=255)
                                                                                    cgTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTTGa                                                              <  7:443598/151‑1 (MQ=255)
                                                                                     gTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGa                                                              >  1:1002077/1‑150 (MQ=255)
                                                                                       tttAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAAtt                                                           >  3:277647/1‑151 (MQ=255)
                                                                                        ttAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAAtt                                                           <  6:1105191/150‑1 (MQ=255)
                                                                                         tAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTgc                                                         >  1:722943/1‑151 (MQ=255)
                                                                                             tGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTgcagc                                                      >  3:1796992/1‑150 (MQ=255)
                                                                                             tGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTgcagc                                                      >  3:833602/1‑150 (MQ=255)
                                                                                                tACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAat                                                   >  1:839297/1‑150 (MQ=255)
                                                                                                 aCCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAata                                                  <  2:984604/150‑1 (MQ=255)
                                                                                                  cccGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATg                                                >  5:603957/1‑151 (MQ=255)
                                                                                                  cccGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATg                                                >  2:1386635/1‑151 (MQ=255)
                                                                                                    cGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGc                                              <  4:123480/151‑1 (MQ=255)
                                                                                                      gCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGa                                            <  7:1262610/151‑1 (MQ=255)
                                                                                                      gCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGAAGCAATATGGCg                                             >  2:1143097/1‑150 (MQ=255)
                                                                                                        tGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACg                                          >  2:1863921/1‑151 (MQ=255)
                                                                                                         gACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGc                                         >  5:1050084/1‑151 (MQ=255)
                                                                                                         gACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGc                                         >  5:733252/1‑151 (MQ=255)
                                                                                                           cgcgGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCa                                       <  6:1519485/151‑1 (MQ=255)
                                                                                                             cgGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCAtt                                     >  8:1013121/1‑151 (MQ=255)
                                                                                                              ggCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTa                                    <  5:891672/151‑1 (MQ=255)
                                                                                                               gCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTAc                                   >  5:1806989/1‑151 (MQ=255)
                                                                                                               gCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTAc                                   <  6:1376903/151‑1 (MQ=255)
                                                                                                                ctctctCCTGCCCTTTTCCCTGCTCCTGGACGGTTTGCCCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCGGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTAc                                   >  2:101217/1‑150 (MQ=255)
                                                                                                                  ctctcCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACa                                  <  6:170249/149‑1 (MQ=255)
                                                                                                                  ctctcCTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACAtt                                >  1:455617/1‑151 (MQ=255)
                                                                                                                      ccTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATtgt                              <  2:258025/149‑1 (MQ=255)
                                                                                                                      ccTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATtgt                              <  1:396058/149‑1 (MQ=255)
                                                                                                                       cTGCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCt                           >  1:1243009/1‑151 (MQ=255)
                                                                                                                         gCCCTTTTCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCt                           <  7:1858982/149‑1 (MQ=255)
                                                                                                                               ttCCCTGCTCCTGGACTGTTTGACCCCTGTTACACGTGAAGTTACCGTAATCAAAAAAGTCTTTTTAAAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCaaa                    <  5:1634729/150‑1 (MQ=255)
                                                                                                                               ttCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCaaa                    <  5:839237/150‑1 (MQ=255)
                                                                                                                               ttCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCaaa                    >  7:890740/1‑150 (MQ=255)
                                                                                                                                tCCCTTCTCCTGTACTGTTTGACCCCTGTTTTATGTGAATTTATCGTAATCAAAAAATTCTGTTTATAAACTAATTCCTAACAACATCTCATTTTTTTGACACTAGAATTGCAGCAATATGGCTACGCCATTACATTGTGCTTCTCaaaa                   <  8:363248/150‑1 (MQ=255)
                                                                                                                                tCCCTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACTCCATTACATTGTGCTTCTCaaaa                   >  4:1662563/1‑150 (MQ=255)
                                                                                                                                 cccTGGTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAg                 <  6:626243/151‑1 (MQ=255)
                                                                                                                                 cccTGCTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCaaaaa                  <  6:255834/150‑1 (MQ=255)
                                                                                                                                      cTCCTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAAGCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGACATTACATTGTGGTTCTCAAAAAGCATa             >  7:1020538/1‑150 (MQ=255)
                                                                                                                                         cTGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACg         >  1:1167847/1‑151 (MQ=255)
                                                                                                                                          tGGACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGt        >  8:687010/1‑151 (MQ=255)
                                                                                                                                           ggACGGTTTGACTCCTGTTACAAGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTATCAACATCGCATTTTTTTGACACTAGAATTGCAGCATTATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAAc          <  4:304233/148‑1 (MQ=255)
                                                                                                                                            gACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTg      <  2:803181/151‑1 (MQ=255)
                                                                                                                                            gACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTg      <  2:482724/151‑1 (MQ=255)
                                                                                                                                            gACGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTg      >  4:475310/1‑151 (MQ=255)
                                                                                                                                              cGGTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTGGc    >  3:234367/1‑151 (MQ=255)
                                                                                                                                               ggTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTGGCt   >  2:1704836/1‑151 (MQ=255)
                                                                                                                                                gTTTGACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTGGCt   >  5:1523051/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTGGCAGAAAAGGTGAAAACCAAAGCCTGCGATCTGGTGTTAAAACAGGGGCTGAACTTCATTTCCTGAACCGCGAAACGCCGCGTTTTAGCTGTTACCCGGCTGACGCGGCTCTCTCCTGCCCTTTTCCCTGCTCCTGGACGGTTTTACCCCTGTTACACGGGAAGTGACCGTAATCAAAAAAGTCTGTTTATAAACTAATTCCTAACAACATCGCATTTTTTTGACACTAGAATTGCAGCAATATGGCGACGCCATTACATTGTGCTTCTCAAAAAGCTTAAACGTTGGCTA  >  DH1/3913265‑3913559

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: