Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,547,561 C→G A165P (GCA→CCA)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,547,5610CG99.2% 2340.7 / 4.6 613A165P (GCA→CCA) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base G (268/336);  total (270/343)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 613 total aligned reads displayed.

CCGGACGTTCGTCGATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGCCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAAC  >  DH1/4547414‑4547708
                                                                                                                                                   |                                                                                                                                                     
ccGGACGTTCGTCGATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCa                                                                                                                                                     <  5:251073/150‑1 (MQ=255)
    aCGTTCGTCGATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATcc                                                                                                                                                >  5:1814171/1‑151 (MQ=255)
     cGTTCGTCGATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCa                                                                                                                                               <  1:1344516/151‑1 (MQ=255)
      gTTCGTCGATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATcc                                                                                                                                                >  6:398216/1‑149 (MQ=255)
         cgtcgATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTa                                                                                                                                            <  6:123051/150‑1 (MQ=255)
         cgtcgATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTAc                                                                                                                                           >  2:865224/1‑151 (MQ=255)
         cgtcgATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTAc                                                                                                                                           >  5:680611/1‑151 (MQ=255)
          gtcgATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTAcg                                                                                                                                          <  2:529251/151‑1 (MQ=255)
            cgATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTAcgcg                                                                                                                                        <  2:1166916/151‑1 (MQ=255)
              aTCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTAcgcg                                                                                                                                        >  7:441434/1‑149 (MQ=255)
               tCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGc                                                                                                                                     >  1:1035221/1‑151 (MQ=255)
                cagcagAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGc                                                                                                                                     <  5:561836/150‑1 (MQ=255)
                cagcagAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCa                                                                                                                                    <  7:135636/151‑1 (MQ=255)
                   cagAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGCTCCAGTACGCGAGCAGtt                                                                                                                                 >  6:448698/1‑151 (MQ=255)
                       aCCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATTGCAGATCCAGTACGCGAGCAGTTAAAt                                                                                                                             <  1:1871067/151‑1 (MQ=255)
                       aCCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAAt                                                                                                                             <  1:556270/151‑1 (MQ=255)
                       aCCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTATACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACAACGACCGATAGGTGATCGCAGATCCAGTACGCGAGCAGTTaaa                                                                                                                              >  6:1477769/1‑150 (MQ=255)
                        ccATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATc                                                                                                                            >  2:1431694/1‑151 (MQ=255)
                          aTCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCtt                                                                                                                          >  4:1705721/1‑151 (MQ=255)
                            cAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTc                                                                                                                         <  6:598680/150‑1 (MQ=255)
                            cAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTGCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCa                                                                                                                        <  3:1739684/151‑1 (MQ=255)
                                  acaATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGaa                                                                                                                  <  5:1261268/151‑1 (MQ=255)
                                  acaATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGa                                                                                                                   >  4:1229468/1‑150 (MQ=255)
                                     aTCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCg                                                                                                               <  1:1145194/151‑1 (MQ=255)
                                     aTCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCg                                                                                                               >  3:1010565/1‑151 (MQ=255)
                                       ccGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGt                                                                                                              <  7:1746926/150‑1 (MQ=255)
                                         gggTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTAc                                                                                                           <  3:957240/151‑1 (MQ=255)
                                             ggtTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACg                                                                                                       <  4:18953/151‑1 (MQ=255)
                                                tGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTc                                                                                                    >  4:215525/1‑151 (MQ=255)
                                                 gTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGAGATGGCAGATCCAGTACGCGAGCAGTTAAATCATCAGTAGATACGTAACCACGTTc                                                                                                    >  8:404667/1‑150 (MQ=255)
                                                   aaGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCa                                                                                                  <  2:76411/150‑1 (MQ=255)
                                                   aaGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCAt                                                                                                 <  8:1092705/151‑1 (MQ=255)
                                                    aGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCAt                                                                                                 >  2:208503/1‑150 (MQ=255)
                                                     ccAATGGGCTGCGCAATGTTCTGCAGCAGCATGGGTTTACCGGCTTTCGGCGGTGCCACAATCACACCACGCTTACCACGACCCATAGGTGATGGCACATGCAGTACGCGCGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATAc                                                                                               <  1:399702/150‑1 (MQ=255)
                                                       aaTGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGc                                                                                             >  4:1815885/1‑151 (MQ=255)
                                                        aTGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCa                                                                                            <  5:1090767/151‑1 (MQ=255)
                                                          gCTCTGAGCAATGTTGTGCATCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGTCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGa                                                                                          <  3:798587/151‑1 (MQ=255)
                                                             ctGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGAAGACg                                                                                        >  3:1311903/1‑150 (MQ=255)
                                                                  cAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGtt                                                                                 >  6:175319/1‑151 (MQ=255)
                                                                   aaTGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGttt                                                                                >  7:372179/1‑151 (MQ=255)
                                                                    aTGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTg                                                                               >  5:855989/1‑151 (MQ=255)
                                                                     tGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGc                                                                              >  4:956866/1‑151 (MQ=255)
                                                                           gcagcagcaTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCa                                                                         <  8:1247694/150‑1 (MQ=255)
                                                                           gcagcagcaTGGTTTTACCGGCTTTCGGCGGTGC‑‑CAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGc                                                                       <  7:47183/150‑1 (MQ=255)
                                                                                agcaTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCgggg                                                                   >  7:1550348/1‑151 (MQ=255)
                                                                                 gcaTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGt                                                                  <  4:1100324/151‑1 (MQ=255)
                                                                                 gcaTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGt                                                                  <  6:828516/151‑1 (MQ=255)
                                                                                 gcaTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGt                                                                  >  5:458773/1‑151 (MQ=255)
                                                                                  caTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGtt                                                                 <  5:349234/151‑1 (MQ=255)
                                                                                  caTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGtt                                                                 <  3:45566/151‑1 (MQ=255)
                                                                                   aTGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGtt                                                                 >  8:743744/1‑150 (MQ=255)
                                                                                     ggTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAg                                                              <  7:832175/151‑1 (MQ=255)
                                                                                       ttttACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGtt                                                            >  2:1720251/1‑151 (MQ=255)
                                                                                        tttACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGtt                                                            <  5:1530179/150‑1 (MQ=255)
                                                                                         ttACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTtc                                                           <  7:948205/150‑1 (MQ=255)
                                                                                         ttACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGCCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTtc                                                           <  7:205625/150‑1 (MQ=255)
                                                                                          tACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTtctc                                                         >  8:1475445/1‑151 (MQ=255)
                                                                                          tACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTtct                                                          <  2:10506/150‑1 (MQ=255)
                                                                                          tACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTtct                                                          <  3:1113326/150‑1 (MQ=255)
                                                                                           aCCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCa                                                        <  2:1323354/151‑1 (MQ=255)
                                                                                            ccGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCaa                                                       <  3:1465065/151‑1 (MQ=255)
                                                                                               gCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCaa                                                       >  3:624880/1‑148 (MQ=255)
                                                                                                     ggcggTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCATTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGAttt                                             <  3:741039/151‑1 (MQ=255)
                                                                                                     ggcggTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGAttt                                             <  4:431782/151‑1 (MQ=255)
                                                                                                       cggTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTttg                                           >  7:1118560/1‑151 (MQ=255)
                                                                                                        ggTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTttgt                                          <  3:1607524/151‑1 (MQ=255)
                                                                                                        ggTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTttgt                                          >  3:24445/1‑151 (MQ=255)
                                                                                                        ggTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTttgt                                          <  8:605846/151‑1 (MQ=255)
                                                                                                              acaATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTCCACCTTGCATACGCACGCGAGACGTTTGGTGCCACCGGGTTAACTTGTGAAACACGGTTTTTGTGGCggc                                    >  8:261575/1‑150 (MQ=255)
                                                                                                               caATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGc                                   <  4:324541/151‑1 (MQ=255)
                                                                                                               caATCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGc                                   <  5:1649939/151‑1 (MQ=255)
                                                                                                                aaTCAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCg                                  <  8:1630717/151‑1 (MQ=255)
                                                                                                                   cAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGttt                               >  4:920995/1‑151 (MQ=255)
                                                                                                                   cAGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGttt                               <  7:614074/151‑1 (MQ=255)
                                                                                                                    cgACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGtttt                              <  6:1044973/150‑1 (MQ=255)
                                                                                                                    aGACCACGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGtttt                              <  8:1762369/151‑1 (MQ=255)
                                                                                                                         aCGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGt                         >  2:376316/1‑151 (MQ=255)
                                                                                                                          cGCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGtt                        <  1:1648965/151‑1 (MQ=255)
                                                                                                                           gCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGGAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGttt                       <  3:1896925/151‑1 (MQ=255)
                                                                                                                           gCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGttt                       <  1:231447/151‑1 (MQ=255)
                                                                                                                           gCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGttt                       <  6:210888/151‑1 (MQ=255)
                                                                                                                           gCTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGGTTTCAGGttt                       >  4:811956/1‑151 (MQ=255)
                                                                                                                            cTGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGtt                        <  1:1472839/149‑1 (MQ=255)
                                                                                                                             tGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGt                     >  2:306409/1‑151 (MQ=255)
                                                                                                                             tGACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAATAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGt                     >  6:1391488/1‑151 (MQ=255)
                                                                                                                              gACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTg                      <  8:1922963/149‑1 (MQ=255)
                                                                                                                              gACCACGACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTg                      <  3:1521223/149‑1 (MQ=255)
                                                                                                                                  aCGACCGCTATGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAg                <  2:1617389/151‑1 (MQ=255)
                                                                                                                                    gACCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGtt              <  7:1502390/151‑1 (MQ=255)
                                                                                                                                     agctaTAGGTGATGGCAGATCCAGTACGCGAGCAGTTAATTCTTCATTAGAAGCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGtt              <  6:973636/146‑1 (MQ=255)
                                                                                                                                     aCCGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTa             <  7:1891752/151‑1 (MQ=255)
                                                                                                                                      ccGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTaa            <  6:1630055/151‑1 (MQ=255)
                                                                                                                                       cGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTaa            >  4:1438216/1‑150 (MQ=255)
                                                                                                                                       cGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAAc           <  4:1550176/151‑1 (MQ=255)
                                                                                                                                       cGATAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAAc           <  1:850434/151‑1 (MQ=255)
                                                                                                                                         aTAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACtt         >  5:1896765/1‑151 (MQ=255)
                                                                                                                                         aTAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACtt         >  2:1820160/1‑151 (MQ=255)
                                                                                                                                         aTAGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACtt         >  4:665901/1‑151 (MQ=255)
                                                                                                                                           aGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCg       <  5:797008/151‑1 (MQ=255)
                                                                                                                                           aGGTGATGGCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCg       >  5:143666/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGGACGTTCGTCGATCAGCAGAACCATCAGCACACAATCCGGGTGGTTGTAAGCAATGCTCTGAGCAATGTTCTGCAGCAGCATGGTTTTACCGGCTTTCGGCGGTGCCACAATCAGACCACGCTGACCACGACCGATAGGTGATGCCAGATCCAGTACGCGAGCAGTTAAATCTTCAGTAGAACCGTTACCACGTTCCATACGCAGACGAGAGTTTGCGTGCAGCGGGGTTAAGTTCTCAAAGAGGATTTTGTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAAC  >  DH1/4547414‑4547708

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: