Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,547,816 A→G Y80H (TAC→CAC)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,547,8160AG98.6% 1581.2 / 5.3 426Y80H (TAC→CAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (211/205);  total (215/211)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 426 total aligned reads displayed.

TGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTT  >  DH1/4547669‑4547963
                                                                                                                                                   |                                                                                                                                                   
tgCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGAt                                                                                                                                                  >  1:1704266/1‑151 (MQ=255)
 gCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATg                                                                                                                                                 <  4:819771/151‑1 (MQ=255)
    ggCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGtca                                                                                                                                              <  7:1441022/151‑1 (MQ=255)
       gTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGtcatc                                                                                                                                            >  1:1204029/1‑150 (MQ=255)
        ttttCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGtcatc                                                                                                                                            >  7:1720479/1‑149 (MQ=255)
         tttCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAgg                                                                                                                                         >  4:1330222/1‑151 (MQ=255)
          ttCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAgg                                                                                                                                         <  7:1816528/150‑1 (MQ=255)
           gcAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGGTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGa                                                                                                                                        <  5:17877/149‑1 (MQ=255)
             aGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCg                                                                                                                                     >  1:1044012/1‑151 (MQ=255)
             aGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCg                                                                                                                                     <  3:1318203/151‑1 (MQ=255)
              ggTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACATTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCAAGTGGAAACGTGGATGTCATCAGGACAgg                                                                                                                                    >  2:340679/1‑151 (MQ=255)
               gTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGc                                                                                                                                   <  5:1364362/151‑1 (MQ=255)
                 ttGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCg                                                                                                                                  <  4:26744/150‑1 (MQ=255)
                  tGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCg                                                                                                                                  >  3:152635/1‑149 (MQ=255)
                   gTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAgg                                                                                                                               >  1:47220/1‑151 (MQ=255)
                      gAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGt                                                                                                                              >  4:1649622/1‑149 (MQ=255)
                      gAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTACGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAAGt                                                                                                                              >  1:216406/1‑149 (MQ=255)
                        aGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGa                                                                                                                          <  5:1142270/151‑1 (MQ=255)
                         gTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAg                                                                                                                         <  7:1707846/151‑1 (MQ=255)
                           tAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGc                                                                                                                        <  2:600697/150‑1 (MQ=255)
                             aCTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTGATCAGGACCGGCGAGGTAGGAGCTGt                                                                                                                     >  7:748306/1‑151 (MQ=255)
                             aCTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGt                                                                                                                     >  1:1096130/1‑151 (MQ=255)
                             aCTCCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGt                                                                                                                     <  4:68815/151‑1 (MQ=255)
                                tCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCt                                                                                                                   <  5:987275/150‑1 (MQ=255)
                                 cGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGc                                                                                                                 <  4:1503558/151‑1 (MQ=255)
                                  gTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCg                                                                                                                <  7:1387808/151‑1 (MQ=255)
                                  gTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCg                                                                                                                >  7:1199215/1‑151 (MQ=255)
                                   ttAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCgg                                                                                                               >  4:560927/1‑151 (MQ=255)
                                   ttAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCgg                                                                                                               >  7:499931/1‑151 (MQ=255)
                                    tAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGAAGGAGCTGTCTGCgg                                                                                                               >  2:1323965/1‑150 (MQ=255)
                                    tAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATAACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGCTGTCATCAGGACCGGCGAGGTAGGACCTGTGTGCggc                                                                                                              >  2:1180822/1‑150 (MQ=255)
                                     aaCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGa                                                                                                              >  2:475973/1‑150 (MQ=255)
                                     aaCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGa                                                                                                              >  5:1607513/1‑150 (MQ=255)
                                     aaCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGa                                                                                                              <  6:1175983/150‑1 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                            <  5:342748/151‑1 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                            >  8:1227838/1‑151 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                            <  2:1092067/151‑1 (MQ=255)
                                         ttCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACgg                                                                                                          <  7:851313/150‑1 (MQ=255)
                                         ttCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACg                                                                                                           <  8:1371281/149‑1 (MQ=255)
                                         ttCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACg                                                                                                           <  2:866471/149‑1 (MQ=255)
                                           cagcagCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggag                                                                                                        <  7:954617/150‑1 (MQ=255)
                                            agcagcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCTGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggagga                                                                                                      <  6:1728426/151‑1 (MQ=255)
                                            agcagcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggagga                                                                                                      <  8:1853122/151‑1 (MQ=255)
                                            agcagcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGACGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggagga                                                                                                      <  3:446820/151‑1 (MQ=255)
                                               agcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAc                                                                                                   >  8:60374/1‑151 (MQ=255)
                                                gcgcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAACCGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAcc                                                                                                  >  8:860298/1‑151 (MQ=255)
                                                gcgcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGCACGGAGGAAAcc                                                                                                  >  8:241924/1‑151 (MQ=255)
                                                 cgcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAcc                                                                                                  >  6:466481/1‑150 (MQ=255)
                                                  gcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCaa                                                                                                <  6:1866140/151‑1 (MQ=255)
                                                  gcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCa                                                                                                 >  4:1213329/1‑150 (MQ=255)
                                                      aaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAACCCAAatcc                                                                                            >  4:1464904/1‑151 (MQ=255)
                                                      aaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatcc                                                                                            >  3:679410/1‑151 (MQ=255)
                                                      aaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatcc                                                                                            >  1:1850623/1‑151 (MQ=255)
                                                      aaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTAGGCTAGCGGAAACGTGGATGTCATCAGGACAGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatcc                                                                                            >  7:1555144/1‑151 (MQ=255)
                                                       aTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatcca                                                                                           <  3:963444/151‑1 (MQ=255)
                                                           cGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGTAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCt                                                                                       <  6:771332/151‑1 (MQ=255)
                                                           cGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTCGATGTCATCAGGACCGGCCAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatccatcat                                                                                       >  8:163443/1‑149 (MQ=255)
                                                            gTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTg                                                                                      <  1:307784/151‑1 (MQ=255)
                                                            gTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTg                                                                                      >  5:893540/1‑151 (MQ=255)
                                                            gTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTg                                                                                      >  2:1858180/1‑151 (MQ=255)
                                                               cACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAATCGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCaa                                                                                   <  5:445692/151‑1 (MQ=255)
                                                               cACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCaa                                                                                   <  3:1780322/151‑1 (MQ=255)
                                                                aCCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAat                                                                                  >  5:1509751/1‑151 (MQ=255)
                                                                  ctacttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAatat                                                                                >  5:229370/4‑151 (MQ=255)
                                                                    tcttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAatat                                                                                <  1:1528022/149‑1 (MQ=255)
                                                                     cttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATAtctc                                                                             >  5:748879/1‑151 (MQ=255)
                                                                            cgggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACg                                                                      >  7:316109/1‑151 (MQ=255)
                                                                                cgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCa                                                                   >  6:1436731/1‑150 (MQ=255)
                                                                                 gAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAGACCAAATCCATCCTGCAATATCTCCAGTCCGGCATc                                                                 >  6:912668/1‑151 (MQ=255)
                                                                                 gAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCCTc                                                                 >  5:1102698/1‑151 (MQ=255)
                                                                                 gAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATc                                                                 <  4:435329/151‑1 (MQ=255)
                                                                                  aaTCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATc                                                                 >  4:903141/1‑150 (MQ=255)
                                                                                  aaTCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCa                                                                <  7:1022670/151‑1 (MQ=255)
                                                                                       tACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCAcc                                                              <  1:533122/148‑1 (MQ=255)
                                                                                         ccAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGa                                                         <  3:250637/151‑1 (MQ=255)
                                                                                               aTGGTTTCACCAGTGCGGATTTTGAAACGGCGGATTTTGCTATGGGATACTTTGATCTCATCAGGACCGGCGAGGTAGTAGCTGTCTGCGGAACGGAGGAAACCAACTCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCtt                                                   <  1:829432/151‑1 (MQ=255)
                                                                                               aTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCtt                                                   <  8:1457626/151‑1 (MQ=255)
                                                                                                  gTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGc                                                <  3:1202715/151‑1 (MQ=255)
                                                                                                   tATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGcc                                               >  2:160858/1‑151 (MQ=255)
                                                                                                     tCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCAc                                             >  5:554766/1‑151 (MQ=255)
                                                                                                      cACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCAct                                            <  7:580356/151‑1 (MQ=255)
                                                                                                         cAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCtt                                         <  4:1130575/151‑1 (MQ=255)
                                                                                                         cAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCtt                                         <  6:1570673/151‑1 (MQ=255)
                                                                                                                   gTTGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGtgctgc                               >  2:1711723/1‑151 (MQ=255)
                                                                                                                     tGAAACGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTgctgct                              <  1:950655/150‑1 (MQ=255)
                                                                                                                        aaCGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCGATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGGTTCAg                          <  1:1252403/151‑1 (MQ=255)
                                                                                                                        aaCGGCGGATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGAGTGCTGCTTc                            >  3:1012178/1‑149 (MQ=255)
                                                                                                                           ggcggATTTGGCTAGGGGAAAGGTGGATGTCATCAGGACGGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAATCCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCGACTCTTTGGGTGCTGCTTCAGGAt                       <  8:347241/151‑1 (MQ=255)
                                                                                                                           ggcggATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGAt                       <  7:712859/151‑1 (MQ=255)
                                                                                                                            gcggATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAATTCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATg                      <  8:714320/151‑1 (MQ=255)
                                                                                                                             cggATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATgg                     >  8:1160065/1‑151 (MQ=255)
                                                                                                                             cggATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATgg                     >  2:1795425/1‑151 (MQ=255)
                                                                                                                              ggATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGc                    <  4:1878969/151‑1 (MQ=255)
                                                                                                                               gATTTGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGc                    >  3:1371261/1‑150 (MQ=255)
                                                                                                                                   tGGCTAGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCaaaa                <  2:1404425/150‑1 (MQ=255)
                                                                                                                                        aGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataa            <  4:1814766/149‑1 (MQ=255)
                                                                                                                                        aGGGGAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATg          >  3:1098351/1‑151 (MQ=255)
                                                                                                                                         ggggAAACGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGt         <  1:1381807/151‑1 (MQ=255)
                                                                                                                                         ggggAAACGTAGATGTCATCAGGAGCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCGATCCTGCAATATCTCCAGTACGCCATCACGAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGt         <  6:569589/151‑1 (MQ=255)
                                                                                                                                               aCGTGGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCt   >  3:725814/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTT  >  DH1/4547669‑4547963

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