Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 1,284,869 G→A G290S (GGT→AGT)  dapA → 4‑hydroxy‑tetrahydrodipicolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,284,8690GA98.8% 1554.9 / 3.0 406G290S (GGT→AGT) dapA4‑hydroxy‑tetrahydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (186/213);  total (190/216)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 406 total aligned reads displayed.

CAACAAACTATTTGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCGGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAGTGGTGATGAAGC  >  DH1/1284721‑1285015
                                                                                                                                                    |                                                                                                                                                  
  acaaACTATTTGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGttt                                                                                                                                                <  7:840607/151‑1 (MQ=255)
   caaACTATTTGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTg                                                                                                                                               >  1:1272454/1‑151 (MQ=255)
     aaCTATTTGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCt                                                                                                                                             <  8:752774/151‑1 (MQ=255)
     aaCTATTTGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCt                                                                                                                                             >  7:593492/1‑151 (MQ=255)
       cTATTTGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGATTGCTGt                                                                                                                                           >  8:626166/1‑151 (MQ=255)
           ttGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTaa                                                                                                                                         >  5:948328/1‑149 (MQ=255)
               cGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGttt                                                                                                                                    >  1:1504723/1‑150 (MQ=255)
               cGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTa                                                                                                                                   <  5:1172704/151‑1 (MQ=255)
                gAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGttt                                                                                                                                    >  1:187178/1‑149 (MQ=255)
                gAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGttt                                                                                                                                    >  7:1189568/1‑149 (MQ=255)
                gAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAg                                                                                                                                  <  3:378561/151‑1 (MQ=255)
                  aCCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAggg                                                                                                                                >  4:659159/1‑151 (MQ=255)
                  aCCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAggg                                                                                                                                <  6:32951/151‑1 (MQ=255)
                  aCCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAggg                                                                                                                                <  8:1288772/151‑1 (MQ=255)
                   cccAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGga                                                                                                                               <  3:484072/151‑1 (MQ=255)
                    ccaatccaatCCCGGTGAAATGGGCATGTAAGGAACTGGGGCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGgag                                                                                                                              <  6:1511360/151‑1 (MQ=255)
                     caatccaatcCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGgaga                                                                                                                             <  7:1496359/151‑1 (MQ=255)
                       atccaatccCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGAtt                                                                                                                           <  1:712623/151‑1 (MQ=255)
                         gcattccCGGTGATAAGGGCATGTAAGTTACTGGGTCTTGTGGCGTCCGATACGATGCGCCTGCCATTGTCACCAATCTCCGAGAGTGGTCGTGAGACGGTAAGAGCGTCGCTTAAGAATGCCAGTTTGCTGTAAAGTTTAGGGAGtttt                                                                                                                          <  5:747346/146‑1 (MQ=255)
                         ccaatccCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTg                                                                                                                         >  3:1004644/1‑151 (MQ=255)
                         ccaatccCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCGGCCAATGACACCAATCACCGCCAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTg                                                                                                                         >  6:739603/1‑151 (MQ=255)
                          caatccCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTg                                                                                                                         >  5:1010949/1‑150 (MQ=255)
                          caatccCGGTGAAATGGGCATGTAAGGAACTGGGGCTTGTGGCGACCGAGACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGa                                                                                                                        <  6:996554/151‑1 (MQ=255)
                            atccCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATg                                                                                                                      <  5:210524/151‑1 (MQ=255)
                             tccCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATgg                                                                                                                     >  8:1043729/1‑151 (MQ=255)
                              cccGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATgg                                                                                                                     <  1:410675/150‑1 (MQ=255)
                              cccGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGc                                                                                                                    >  1:1000553/1‑151 (MQ=255)
                                cGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCtt                                                                                                                  <  6:54940/151‑1 (MQ=255)
                                 ggTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGC‑‑TAAAGTTTAGGGAGATTTGATGGCTTAct                                                                                                               <  3:1097328/151‑1 (MQ=255)
                                    gAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTAct                                                                                                               <  4:161078/150‑1 (MQ=255)
                                    gAAATGGGCATGTAAGGAACTGGGGCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTAct                                                                                                               <  1:1109933/150‑1 (MQ=255)
                                     aaaTGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTActct                                                                                                             <  3:592078/151‑1 (MQ=255)
                                       aTGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGt                                                                                                           <  4:359002/151‑1 (MQ=255)
                                         gggCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGtt                                                                                                          <  8:217504/150‑1 (MQ=255)
                                         gggCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGtt                                                                                                          <  7:697461/150‑1 (MQ=255)
                                              tGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCaaa                                                                                                      >  6:81321/1‑149 (MQ=255)
                                                 aaGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGCCGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTcg                                                                                                 >  3:442904/1‑151 (MQ=255)
                                                 aaGGAACTGGGTCTTGTGGCGACCGATACGCTGCCCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTcg                                                                                                 >  2:152581/1‑151 (MQ=255)
                                                      ccTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACTGTGGTCGGGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCt                                                                                            <  2:1442834/150‑1 (MQ=255)
                                                      aCTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCt                                                                                            <  3:1425610/151‑1 (MQ=255)
                                                       cTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTg                                                                                           >  1:293606/1‑151 (MQ=255)
                                                         gggTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTg                                                                                           <  8:174212/149‑1 (MQ=255)
                                                         gggTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGc                                                                                         <  8:803236/151‑1 (MQ=255)
                                                             cTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCaa                                                                                       >  6:927422/1‑149 (MQ=255)
                                                               tgtgGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAg                                                                                     <  7:530225/149‑1 (MQ=255)
                                                                gtgGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGtt                                                                                  <  8:1549980/151‑1 (MQ=255)
                                                                 tgGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTg                                                                                 >  2:827866/1‑151 (MQ=255)
                                                                 gggCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTg                                                                                 <  6:1126601/150‑1 (MQ=255)
                                                                   gCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCg                                                                               <  2:416426/151‑1 (MQ=255)
                                                                    cGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCgg                                                                              <  3:1654922/151‑1 (MQ=255)
                                                                     gACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCggg                                                                             >  5:467948/1‑151 (MQ=255)
                                                                      aCCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGgt                                                                            <  2:579475/151‑1 (MQ=255)
                                                                        cGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGgtgt                                                                          >  7:309280/1‑151 (MQ=255)
                                                                         gATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGtt                                                                         <  8:890236/151‑1 (MQ=255)
                                                                          aTACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGttt                                                                        <  4:1239295/151‑1 (MQ=255)
                                                                          aTACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGttt                                                                        <  1:755668/151‑1 (MQ=255)
                                                                              gCTGCCCATGCCAATGACACCAATCACCGACAGTGGACCTCAGACGGTCAGAGCGGCGCTAAAGAATGCAAGATTCCTGTAAAGATTACGGAGAATAGATGGCTTACTCTGTTCAAAAGTCGAGACTGCCAAAGGTAGCGGCTGttaccct                                                                    >  4:1383359/1‑146 (MQ=255)
                                                                                tGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTg                                                                  >  5:1476552/1‑151 (MQ=255)
                                                                                 gcgcCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCtt                                                                   <  3:717687/149‑1 (MQ=255)
                                                                                 gcgcCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTg                                                                  <  1:143168/150‑1 (MQ=255)
                                                                                  cgcCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGtt                                                                >  5:1526309/1‑151 (MQ=255)
                                                                                  cgcCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGtt                                                                <  4:870509/151‑1 (MQ=255)
                                                                                   gcCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGttt                                                               <  2:1075067/151‑1 (MQ=255)
                                                                                    ccTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCGTGGCAAAGGTTGCGGGTGTTTCGCTTGtttt                                                              >  7:1231351/1‑151 (MQ=255)
                                                                                     cTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGtttt                                                              >  4:503299/1‑150 (MQ=255)
                                                                                      tGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGGTttat                                                            >  6:713195/1‑151 (MQ=255)
                                                                                       gCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTttatt                                                           <  3:244154/151‑1 (MQ=255)
                                                                                       gCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTtta                                                             <  5:526876/149‑1 (MQ=255)
                                                                                         cAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTGGCGGGGGTTTCGCTGGTTTTATTAc                                                         >  4:195108/1‑151 (MQ=255)
                                                                                          aaTGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACt                                                        <  5:1660312/151‑1 (MQ=255)
                                                                                             gACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGc                                                     >  7:621955/1‑151 (MQ=255)
                                                                                              acacCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCGGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCt                                                    >  5:1286043/1‑151 (MQ=255)
                                                                                               cacCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCTAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTg                                                   >  3:300817/1‑151 (MQ=255)
                                                                                               cacCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTg                                                   <  2:796704/151‑1 (MQ=255)
                                                                                                acCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGc                                                  <  1:1177241/151‑1 (MQ=255)
                                                                                                  cAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCt                                                >  2:296446/1‑151 (MQ=255)
                                                                                                      cACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTa                                             <  3:1296817/150‑1 (MQ=255)
                                                                                                      cACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTa                                             <  7:1356487/150‑1 (MQ=255)
                                                                                                      cACCGACAGGGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTa                                             <  3:891121/150‑1 (MQ=255)
                                                                                                         cGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTc                                         <  5:630941/151‑1 (MQ=255)
                                                                                                          gACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTCAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCt                                        <  3:78699/151‑1 (MQ=255)
                                                                                                           aCAGTGGTCGTGAGACGGTCAGAGCTGCTCTTAATCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTCCCTGTAGTTCTg                                       <  6:443471/151‑1 (MQ=255)
                                                                                                           aCAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTg                                       <  3:1539042/151‑1 (MQ=255)
                                                                                                                   cGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGc                               >  5:1145896/1‑151 (MQ=255)
                                                                                                                      gagaCGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCtat                            <  2:187536/151‑1 (MQ=255)
                                                                                                                        gaCGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTCGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCtata                           >  6:1619108/1‑150 (MQ=255)
                                                                                                                         aCGGTCAGAGCGGCTCTTAAGCATGCGAGTTTGCTGTAAAGTTTAGGGAGATTTGATCGCTTACTCTGTTAAAAAGTCGCTCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTAATCGCTGCCTGTAGTTCTGACTCACGCTATaa                          <  6:1294805/150‑1 (MQ=255)
                                                                                                                         aCGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCtata                           >  3:664910/1‑149 (MQ=255)
                                                                                                                          cGGTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTTGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTTGCAAAGGTTGCTGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGc                        <  4:725388/151‑1 (MQ=255)
                                                                                                                           ggTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCg                       >  5:783094/1‑151 (MQ=255)
                                                                                                                            gTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCgt                      >  2:1360278/1‑151 (MQ=255)
                                                                                                                            gTCAGAGCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATCGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCgt                      <  4:1105095/151‑1 (MQ=255)
                                                                                                                               agagCGGCGCTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAATTCTGACTCACGCTATAAGCgtc                     >  8:396364/1‑149 (MQ=255)
                                                                                                                                      cgcTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAgtg            >  1:986285/1‑151 (MQ=255)
                                                                                                                                       gcTTAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAgtgg           >  5:1391364/1‑151 (MQ=255)
                                                                                                                                          tAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAgtggtg         >  2:726348/1‑150 (MQ=255)
                                                                                                                                          tAAGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAgtggtg         <  1:634442/150‑1 (MQ=255)
                                                                                                                                           aaGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAGTGGtga        <  6:118245/150‑1 (MQ=255)
                                                                                                                                           aaGCATGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAGTGGtga        <  4:567179/150‑1 (MQ=255)
                                                                                                                                               aTGCCAGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAGTGGTGATGAAg   <  3:1253371/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAACAAACTATTTGTCGAACCCAATCCAATCCCGGTGAAATGGGCATGTAAGGAACTGGGTCTTGTGGCGACCGATACGCTGCGCCTGCCAATGACACCAATCACCGACAGTGGTCGTGAGACGGTCAGAGCGGCGCTTAAGCATGCCGGTTTGCTGTAAAGTTTAGGGAGATTTGATGGCTTACTCTGTTCAAAAGTCGCGCCTGGCAAAGGTTGCGGGTGTTTCGCTTGTTTTATTACTCGCTGCCTGTAGTTCTGACTCACGCTATAAGCGTCAGGTCAGTGGTGATGAAGC  >  DH1/1284721‑1285015

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