Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,547,815 T→C Y80C (TAC→TGC)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,547,8150TC96.1% 1176.7 / 21.1 332Y80C (TAC→TGC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (7/6);  new base C (166/152);  total (174/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 332 total aligned reads displayed.

GTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCT  >  DH1/4547667‑4547962
                                                                                                                                                    |                                                                                                                                                   
   gCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGAt                                                                                                                                                 <  5:651915/150‑1 (MQ=255)
    cGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGt                                                                                                                                               <  6:1076964/151‑1 (MQ=255)
     gggCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtc                                                                                                                                              >  2:559061/1‑151 (MQ=255)
      gtcGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtca                                                                                                                                             <  5:123229/149‑1 (MQ=255)
      ggCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtca                                                                                                                                             <  5:70395/151‑1 (MQ=255)
       gCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtcat                                                                                                                                            >  8:558764/1‑151 (MQ=255)
        cGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGCTGtcatc                                                                                                                                           >  7:1247871/1‑151 (MQ=255)
        cGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtcatc                                                                                                                                           >  6:520563/1‑151 (MQ=255)
          ttttCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtcatc                                                                                                                                           >  7:1099591/1‑149 (MQ=255)
            ttCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAggc                                                                                                                                       >  5:489716/1‑150 (MQ=255)
             tCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGATACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGa                                                                                                                                       <  4:197741/150‑1 (MQ=255)
              cAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGAcc                                                                                                                                     <  8:1087295/151‑1 (MQ=255)
                ggTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCgg                                                                                                                                   <  6:877563/151‑1 (MQ=255)
                   ttGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCg                                                                                                                                 >  8:971584/1‑150 (MQ=255)
                    tGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCg                                                                                                                                 >  2:1059326/1‑149 (MQ=255)
                    tGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGa                                                                                                                                <  5:913558/150‑1 (MQ=255)
                     gTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAgg                                                                                                                              >  2:216235/1‑151 (MQ=255)
                      tCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAACCGCAGATGTCATCAGGACCGGCGAGGt                                                                                                                             >  6:694726/1‑151 (MQ=255)
                       cGAAGTTAACTTCGTTAACTTTCAGCAGCTCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGTGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGt                                                                                                                             <  4:93892/150‑1 (MQ=255)
                       cGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTa                                                                                                                            <  5:989691/151‑1 (MQ=255)
                         aaGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAgg                                                                                                                          >  4:1193014/1‑151 (MQ=255)
                         aaGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAgg                                                                                                                          >  3:769276/1‑151 (MQ=255)
                             tAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGTTAGGAGCt                                                                                                                      >  1:516699/1‑151 (MQ=255)
                              aaCTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTg                                                                                                                     >  7:888853/1‑151 (MQ=255)
                               aCTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGt                                                                                                                    <  7:628344/151‑1 (MQ=255)
                                 ttCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTcg                                                                                                                  >  5:1015358/1‑150 (MQ=255)
                                  tCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCt                                                                                                                  <  2:883506/150‑1 (MQ=255)
                                  tCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCt                                                                                                                  <  3:672680/150‑1 (MQ=255)
                                   cGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTg                                                                                                                 <  1:1117680/150‑1 (MQ=255)
                                   cGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGc                                                                                                                >  1:3617/1‑151 (MQ=255)
                                   cGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGc                                                                                                                >  4:381979/1‑151 (MQ=255)
                                    gTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCg                                                                                                               >  6:1200888/1‑151 (MQ=255)
                                        aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                           >  4:573576/1‑151 (MQ=255)
                                          tttCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACgg                                                                                                         >  3:355663/1‑151 (MQ=255)
                                           ttCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACg                                                                                                          <  5:859561/149‑1 (MQ=255)
                                            tCAGCAGCGCAAAATAGCGGTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACgga                                                                                                        <  7:1169092/150‑1 (MQ=255)
                                             cagcagCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggagg                                                                                                      >  3:916792/1‑151 (MQ=255)
                                             cagcagCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggag                                                                                                       <  1:1180528/150‑1 (MQ=255)
                                              agcagcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggagga                                                                                                     <  4:321916/151‑1 (MQ=255)
                                                 agcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAc                                                                                                  <  2:399735/151‑1 (MQ=255)
                                                 agcGCAAAATAGCGTTCACCTTCTTTCGGCGGCCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGAGTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAc                                                                                                  <  6:296123/151‑1 (MQ=255)
                                                  gcgcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAcc                                                                                                 <  3:879991/151‑1 (MQ=255)
                                                   cgcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGCAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAcc                                                                                                 >  3:1232183/1‑150 (MQ=255)
                                                    gcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCaa                                                                                               <  1:641096/151‑1 (MQ=255)
                                                    gcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCaa                                                                                               <  7:313130/151‑1 (MQ=255)
                                                    gcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCCGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGAAGCAGCTGTCTGCGGAAGCGACGAAAcaa                                                                                                >  2:277243/1‑148 (MQ=255)
                                                     cAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCaa                                                                                               >  8:880489/1‑150 (MQ=255)
                                                      aaaaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAATCCAAat                                                                                             <  5:215528/151‑1 (MQ=255)
                                                       aaaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatc                                                                                            <  5:840818/151‑1 (MQ=255)
                                                           aGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGTAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatccat                                                                                         <  5:25584/150‑1 (MQ=255)
                                                           aGCGTGCCCCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGTATTTGGCTAGGGGAAACGCAGATGTCATCAGTACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatccatc                                                                                        <  4:275212/151‑1 (MQ=255)
                                                             cGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGTCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCt                                                                                      >  4:124280/1‑151 (MQ=255)
                                                              gTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTg                                                                                     >  3:563439/1‑151 (MQ=255)
                                                               ttCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGc                                                                                    >  1:997667/1‑151 (MQ=255)
                                                                 cACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGACCGGAGGAAACCAAATCCATCCTGCa                                                                                   >  2:1185016/1‑150 (MQ=255)
                                                                 cACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCa                                                                                   >  3:1055318/1‑150 (MQ=255)
                                                                      tcttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATAtct                                                                             >  4:1102170/1‑151 (MQ=255)
                                                                         ttCGGCGGGCGAATCTTACCAGAGATGGTAACACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTcc                                                                           >  7:1030101/1‑150 (MQ=255)
                                                                          tCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTcc                                                                           >  1:841532/1‑149 (MQ=255)
                                                                           cggcggGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGt                                                                        >  1:185516/1‑151 (MQ=255)
                                                                            ggcgggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGt                                                                        >  8:202901/1‑150 (MQ=255)
                                                                             gcgggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAACCCAAATCCATCCTGCAATATCTCCAGTAc                                                                      >  3:156328/1‑151 (MQ=255)
                                                                                         tACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCAc                                                              <  7:15580/147‑1 (MQ=255)
                                                                                           ccAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGGCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGa                                                        <  2:631031/151‑1 (MQ=255)
                                                                                            cAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGat                                                       <  3:606397/151‑1 (MQ=255)
                                                                                             agagATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGat                                                       >  3:510599/1‑150 (MQ=255)
                                                                                              gagaTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGatat                                                     <  7:964243/151‑1 (MQ=255)
                                                                                              gagaTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGatat                                                     >  1:7439/1‑151 (MQ=255)
                                                                                                  tGGTATCAGCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCtt                                                  <  8:1153937/150‑1 (MQ=255)
                                                                                                  tGGTATCACCAGTGCGGAGGTTCAAACGGCGGATTTGGCTAGGGGAACCGCAGATCTCATCAGGACCCGCGAGGTAGGAGGTGTCTGCGGAACCGACGAAACCAAATACATCCTGCAATATCTCCAGAACGCCATCAACAAAGATATCTTc                                                 >  6:956424/1‑151 (MQ=255)
                                                                                                   ggTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCg                                                <  8:28872/151‑1 (MQ=255)
                                                                                                     tATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGcc                                              >  3:423565/1‑151 (MQ=255)
                                                                                                      aTCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCa                                             <  7:288826/151‑1 (MQ=255)
                                                                                                       tCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCAc                                            >  8:425240/1‑151 (MQ=255)
                                                                                                       tCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCAc                                            >  8:473122/1‑151 (MQ=255)
                                                                                                        cACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCAct                                           >  2:1130526/1‑151 (MQ=255)
                                                                                                         aCCAGTGCGGAGGTTGAAACGGCGGATTTGGGTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCActc                                          <  8:834409/151‑1 (MQ=255)
                                                                                                          ccAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCCCCAAAGATATCTTCGCCActct                                         <  5:1049744/151‑1 (MQ=255)
                                                                                                            aGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTTCAATATCTCCAGTACGCCATCACCAAAGATATCTTCTCCACTCttt                                       >  8:1024954/1‑151 (MQ=255)
                                                                                                            aGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCttt                                       <  3:628793/151‑1 (MQ=255)
                                                                                                                 ggAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGtg                                  <  4:410844/151‑1 (MQ=255)
                                                                                                                       tGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCtt                            >  2:235485/1‑151 (MQ=255)
                                                                                                                        gAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTc                           >  3:252030/1‑151 (MQ=255)
                                                                                                                         aaaCGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCa                          <  8:188371/151‑1 (MQ=255)
                                                                                                                          aaCGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAg                         <  3:526051/151‑1 (MQ=255)
                                                                                                                           aCGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAgg                        <  3:1189848/151‑1 (MQ=255)
                                                                                                                           aCGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAgg                        <  6:143653/151‑1 (MQ=255)
                                                                                                                           aCGGCGGAGTTGGCGAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAgg                        <  6:469880/151‑1 (MQ=255)
                                                                                                                              gcggATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATg                     >  7:583227/1‑151 (MQ=255)
                                                                                                                                 gATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCa                  <  1:410000/151‑1 (MQ=255)
                                                                                                                                     tGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCaaaa               <  6:234695/150‑1 (MQ=255)
                                                                                                                                     tGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCaaaa               >  4:843196/1‑150 (MQ=255)
                                                                                                                                       gCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaat             >  7:473705/1‑150 (MQ=255)
                                                                                                                                       gCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaat             >  4:1005212/1‑150 (MQ=255)
                                                                                                                                         tAGGGGAAACGCAGATGTCATCAGCACCGGCGACGAAGGAGCTGACTGCGCAACGGAGGAAACCAAATCCATGCAGCAATATCTCCAGTACGCCAACACCAAAGATATCTTCGCCACTCATTGCGTGCTGCTTCAGCATGCCAAAaataat          >  6:772384/1‑151 (MQ=255)
                                                                                                                                          aGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATg         <  1:1060973/151‑1 (MQ=255)
                                                                                                                                           ggggAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGt        <  4:511352/151‑1 (MQ=255)
                                                                                                                                             ggAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTcc      >  2:131077/1‑151 (MQ=255)
                                                                                                                                               aaaCGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTg    <  7:753361/151‑1 (MQ=255)
                                                                                                                                               aaaCGCAGATGTCATCAGGACCGGCGAGGGAGGAGCTGGCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTg    <  4:1157038/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCT  >  DH1/4547667‑4547962

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