Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,173,457 T→G W169G (TGG→GGG)  melR → melibiose operon regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,173,4570TG98.7% 1790.0 / 4.7 472W169G (TGG→GGG) melRmelibiose operon regulatory protein
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (219/246);  total (222/250)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 495 total aligned reads displayed.

GCATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCTGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCACGTCA  >  DH1/4173308‑4173604
                                                                                                                                                     |                                                                                                                                                   
gCATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTctggctg                                                                                                                                                    >  5:910241/1‑151 (MQ=255)
 cATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTctggcggg                                                                                                                                                   >  8:1617635/1‑151 (MQ=255)
  aTGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCgggg                                                                                                                                                  <  6:905021/151‑1 (MQ=255)
        aTCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCg                                                                                                                                             >  1:1370769/1‑150 (MQ=255)
         tCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGAt                                                                                                                                           >  3:1788718/1‑151 (MQ=255)
          cAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGAt                                                                                                                                           <  3:1728105/150‑1 (MQ=255)
            actCACTGGCCACACAGCAACTTAGCCCGTATGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCt                                                                                                                                        <  1:202344/149‑1 (MQ=255)
            aaTCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCt                                                                                                                                        >  1:1281026/1‑151 (MQ=255)
            aaTCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCt                                                                                                                                        >  7:1374570/1‑151 (MQ=255)
             aTCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTg                                                                                                                                       <  8:1570767/151‑1 (MQ=255)
               cACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGATTTAAACAGTCCGAACGAGCTAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGt                                                                                                                                     >  8:372332/1‑151 (MQ=255)
                aCTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTc                                                                                                                                    >  8:1714666/1‑151 (MQ=255)
                   ggCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCa                                                                                                                                   <  4:881857/149‑1 (MQ=255)
                   ggCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCa                                                                                                                                   <  2:571588/149‑1 (MQ=255)
                     ccACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATa                                                                                                                                >  4:220499/1‑150 (MQ=255)
                      cacacaGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATaaa                                                                                                                              <  8:867930/151‑1 (MQ=255)
                      cacacaGCAACTTACCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATaaa                                                                                                                              <  2:191870/151‑1 (MQ=255)
                          caGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCt                                                                                                                          <  6:793739/151‑1 (MQ=255)
                           aGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTc                                                                                                                         >  5:441809/1‑151 (MQ=255)
                             cAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTcgc                                                                                                                       <  5:850305/151‑1 (MQ=255)
                              aaCTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGCTTCTGGTCAATAAAACCTcgcg                                                                                                                      >  4:1278875/1‑151 (MQ=255)
                              aaCTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTcgcg                                                                                                                      >  5:1344942/1‑151 (MQ=255)
                                  tAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACTAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGcac                                                                                                                   <  7:433925/150‑1 (MQ=255)
                                  tAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGcacc                                                                                                                  >  6:16949/1‑150 (MQ=255)
                                      ccGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGcacacaca                                                                                                              <  7:610170/151‑1 (MQ=255)
                                      ccGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGcacacac                                                                                                               >  3:361714/1‑150 (MQ=255)
                                           tGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAAc                                                                                                         >  1:469716/1‑151 (MQ=255)
                                             aaGTGCGCCGCTGGCAGCATGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAg                                                                                                       <  3:250270/151‑1 (MQ=255)
                                              aGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGc                                                                                                      <  4:1105084/151‑1 (MQ=255)
                                                tGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGt                                                                                                    >  4:1330518/1‑151 (MQ=255)
                                                  cgccgcTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGtct                                                                                                  <  7:339151/151‑1 (MQ=255)
                                                        tGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGcc                                                                                            >  2:408351/1‑151 (MQ=255)
                                                            agcagTAATTAAACAGGCCGAACGAGAAAATTCGCCAGCTCGCAATTGATTAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgc                                                                                        <  5:677979/151‑1 (MQ=255)
                                                            agcagGAATTAAACAGTCCGATCGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgc                                                                                        <  7:760014/151‑1 (MQ=255)
                                                            agcagGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgc                                                                                        >  3:732247/1‑151 (MQ=255)
                                                            agcagGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgc                                                                                        >  7:1071964/1‑151 (MQ=255)
                                                             gcagGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgcg                                                                                       <  7:1636605/151‑1 (MQ=255)
                                                             gcagGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATgcg                                                                                       <  1:281736/151‑1 (MQ=255)
                                                               agGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCa                                                                                     <  7:1237571/151‑1 (MQ=255)
                                                                ggAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCaa                                                                                    <  1:1614108/151‑1 (MQ=255)
                                                                   aTTAAACAGTGCGAACGCGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAAttt                                                                                 <  4:1031405/151‑1 (MQ=255)
                                                                     tAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTAAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAAtttt                                                                                <  3:103380/150‑1 (MQ=255)
                                                                        aCAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGt                                                                            >  2:1364294/1‑151 (MQ=255)
                                                                          aGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCTGGGAACCGATTCTGGGCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTa                                                                          <  2:1742460/151‑1 (MQ=255)
                                                                          aGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTa                                                                          <  7:1764665/151‑1 (MQ=255)
                                                                            gccGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGc                                                                        <  1:718973/150‑1 (MQ=255)
                                                                             ccGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGc                                                                        <  6:1139308/150‑1 (MQ=255)
                                                                              cGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGcc                                                                       >  2:543052/1‑150 (MQ=255)
                                                                              cGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCa                                                                      <  7:900693/151‑1 (MQ=255)
                                                                                 aCGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGAt                                                                   <  8:1192537/151‑1 (MQ=255)
                                                                                  cGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATg                                                                  <  5:1503017/151‑1 (MQ=255)
                                                                                   gAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGc                                                                 <  5:476897/151‑1 (MQ=255)
                                                                                     gCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTg                                                               >  6:221889/1‑151 (MQ=255)
                                                                                     ccAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTg                                                               <  1:95441/150‑1 (MQ=255)
                                                                                      cAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTgg                                                              <  4:321111/151‑1 (MQ=255)
                                                                                       aaaTTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTggg                                                             >  1:1183416/1‑151 (MQ=255)
                                                                                       aaaTTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTggg                                                             >  4:1373638/1‑151 (MQ=255)
                                                                                       aaaTTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTggg                                                             >  8:1388913/1‑151 (MQ=255)
                                                                                         aTTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCt                                                           <  7:1537054/151‑1 (MQ=255)
                                                                                             gCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTAt                                                       >  3:1565664/1‑151 (MQ=255)
                                                                                               cAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAAGCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCTGATGCTGGGCTTTATTg                                                     <  8:1250324/151‑1 (MQ=255)
                                                                                                aGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGc                                                    <  3:1278723/151‑1 (MQ=255)
                                                                                                aGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGAGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGc                                                    <  3:908256/151‑1 (MQ=255)
                                                                                                 gCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCACACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGcc                                                   <  8:1072960/151‑1 (MQ=255)
                                                                                                   tCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCg                                                  >  2:756133/1‑150 (MQ=255)
                                                                                                   tCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCg                                                  >  2:1192339/1‑150 (MQ=255)
                                                                                                    cGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGa                                                 >  1:1449141/1‑150 (MQ=255)
                                                                                                    cGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGa                                                 >  3:1387392/1‑150 (MQ=255)
                                                                                                    cGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGa                                                 >  8:6527/1‑150 (MQ=255)
                                                                                                    cGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGa                                                 <  2:1447647/150‑1 (MQ=255)
                                                                                                    cGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGa                                                 <  7:21734/150‑1 (MQ=255)
                                                                                                        aTTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACt                                            <  6:1232704/151‑1 (MQ=255)
                                                                                                             tgaAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCGTTATTGCCGAATACTATGAt                                       >  6:1466096/1‑151 (MQ=255)
                                                                                                               aaaTTGGCCTGATGATCAAGCGATTTAGCCTCTCTGGCGGTGAACCGATTCTGTTCAATTAAACATCGCGCACACACAAAAACAGCGTCTCGCTCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCa                                     <  2:909960/151‑1 (MQ=255)
                                                                                                                aaTTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAg                                    <  5:316216/151‑1 (MQ=255)
                                                                                                                 aTTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAgg                                   >  1:1015131/1‑151 (MQ=255)
                                                                                                                 aTTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCAGGTCAATAAAACCTCGCGCACACAGAAAAACAGCGTCTCGCCCCAAGGGCAATTTTATGTAAGCCAGATGCTGGGCTTTATTGCCGAAAACAAGGATCAgc                                   >  4:777042/1‑150 (MQ=255)
                                                                                                                  ttGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGgc                                  >  7:971734/1‑151 (MQ=255)
                                                                                                                   tGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGgcg                                 >  4:1764318/1‑151 (MQ=255)
                                                                                                                   tGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGgc                                  <  4:1447625/150‑1 (MQ=255)
                                                                                                                   tGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGgc                                  <  5:1723115/150‑1 (MQ=255)
                                                                                                                      ccTGATGCTCAAGCGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCTCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTg                              <  5:1105744/151‑1 (MQ=255)
                                                                                                                                   cGATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGAtg                 >  4:573095/1‑151 (MQ=255)
                                                                                                                                    gATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAACCCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGAtgt                >  6:626725/1‑151 (MQ=255)
                                                                                                                                    gATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGAtgt                >  5:1027067/1‑151 (MQ=255)
                                                                                                                                    gATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGAtgt                <  3:1167709/151‑1 (MQ=255)
                                                                                                                                    gATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGAtgt                >  3:518798/1‑151 (MQ=255)
                                                                                                                                    gATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGAtgt                <  4:4622/151‑1 (MQ=255)
                                                                                                                                    gATTTAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGAtgt                >  8:778691/1‑151 (MQ=255)
                                                                                                                                       ttAGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGATAACTATGATCAGGCGCTGACCATAAACGATGTGGc             >  2:82119/1‑151 (MQ=255)
                                                                                                                                         aGCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCt            <  5:167496/150‑1 (MQ=255)
                                                                                                                                          gCCTCTCTGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGa          <  8:207952/151‑1 (MQ=255)
                                                                                                                                           ccTCTCTGGCGGGGAACGGATTCTGGTCAATAAAACCTCGCGCACACACAAAAAGAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGGGCTGACCATCAACGATGTGGCTGAg         <  6:472933/151‑1 (MQ=255)
                                                                                                                                             tctctGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGc        >  5:1246568/1‑150 (MQ=255)
                                                                                                                                              ctctGGCGGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGc        <  4:52638/149‑1 (MQ=255)
                                                                                                                                                ctggcgggGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCa       <  2:1796218/148‑1 (MQ=255)
                                                                                                                                                ctggcgggGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCAc      <  2:1546198/149‑1 (MQ=255)
                                                                                                                                                  ggcgggGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACCGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAATCTATGATCACGCGCTGACCATCAACGATGTGGCTGAGCACGTc   <  6:1605118/150‑1 (MQ=255)
                                                                                                                                                  ggcgggGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATNCGCNNNNNTNNNTNNGANAGANGCTGGGANTTAATGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCAAGTc   >  2:1004632/1‑150 (MQ=255)
                                                                                                                                                  ggcgggGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGGCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCACGTCa  <  6:323434/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCATGGTGATCAAATCACTGGCCACACAGCAACTTAGCCCGTTTGAAGTGCGCCGCTGGCAGCAGGAATTAAACAGTCCGAACGAGCAAATTCGCCAGCTCGCCATTGATGAAATTGGCCTGATGCTCAAGCGATTTAGCCTCTCTGGCTGGGAACCGATTCTGGTCAATAAAACCTCGCGCACACACAAAAACAGCGTCTCGCGCCATGCGCAATTTTATGTTAGCCAGATGCTGGGCTTTATTGCCGAAAACTATGATCAGGCGCTGACCATCAACGATGTGGCTGAGCACGTCA  >  DH1/4173308‑4173604

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: