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,8150TC95.7% 1631.7 / 36.1 461Y80C (TAC→TGC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base C (220/220);  total (231/230)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 461 total aligned reads displayed.

TGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCT  >  DH1/4547669‑4547962
                                                                                                                                                  |                                                                                                                                                   
tgCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGa                                                                                                                                                  <  3:273669/150‑1 (MQ=255)
    ggCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAATATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTCGCTAGGGGAAACGTAGATGtca                                                                                                                                             <  6:846123/151‑1 (MQ=255)
    ggCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTGCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtca                                                                                                                                             <  2:45696/151‑1 (MQ=255)
      cGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtcatc                                                                                                                                           >  5:704915/1‑151 (MQ=255)
       gTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGtcatc                                                                                                                                           >  3:398924/1‑150 (MQ=255)
         tttCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAg                                                                                                                                         <  7:1450477/150‑1 (MQ=255)
         tttCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAg                                                                                                                                         <  4:1588481/150‑1 (MQ=255)
            catgTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGAcc                                                                                                                                     <  6:521236/148‑1 (MQ=255)
            cAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGAcc                                                                                                                                     >  5:1769917/1‑151 (MQ=255)
               gTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGATGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGc                                                                                                                                  <  8:1040811/151‑1 (MQ=255)
                 ttGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAGACGGCGGATTTGGCTAGGGGAAACGCAGCTGTCATCAGGGCCGGCgg                                                                                                                                >  6:1798462/1‑150 (MQ=255)
                   gTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAgg                                                                                                                              <  5:614462/151‑1 (MQ=255)
                    tCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGCCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGt                                                                                                                             >  8:1740569/1‑151 (MQ=255)
                       aaGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAgg                                                                                                                          >  5:300068/1‑151 (MQ=255)
                       aaGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAgg                                                                                                                          <  6:1061718/151‑1 (MQ=255)
                       aaGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGATAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAgg                                                                                                                          >  4:201150/1‑151 (MQ=255)
                        aGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAgg                                                                                                                          >  1:582271/1‑150 (MQ=255)
                        aGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAgg                                                                                                                          >  7:1668972/1‑150 (MQ=255)
                         gTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAg                                                                                                                        <  7:689374/151‑1 (MQ=255)
                         gTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAg                                                                                                                        <  7:636866/151‑1 (MQ=255)
                          ttAACTTCGTTAACTTTCAGAAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGTGGAAACGCAGATGTCATCAGGAACGGCGAGGTAGGAGc                                                                                                                       >  3:1307537/1‑151 (MQ=255)
                            aaCTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTg                                                                                                                     >  5:880232/1‑151 (MQ=255)
                            aaCTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTg                                                                                                                     >  6:1316616/1‑151 (MQ=255)
                             aCTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGt                                                                                                                    >  8:1230567/1‑151 (MQ=255)
                                 cGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCTGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTg                                                                                                                 <  8:347164/150‑1 (MQ=255)
                                 cGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGc                                                                                                                >  3:693246/1‑151 (MQ=255)
                                 cGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGc                                                                                                                <  3:616246/151‑1 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                           <  6:437893/151‑1 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                           <  2:1551892/151‑1 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                           <  7:1094597/151‑1 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                           <  2:1201461/151‑1 (MQ=255)
                                      aCTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAAc                                                                                                           <  4:1216488/151‑1 (MQ=255)
                                       cTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACg                                                                                                          >  5:1549986/1‑151 (MQ=255)
                                        tttCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACgg                                                                                                         >  4:593906/1‑151 (MQ=255)
                                        tttCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACgg                                                                                                         >  4:1294263/1‑151 (MQ=255)
                                           cagcagCGCAAAATAGCGTTCAGCTTGTTTCGGCGGGCTAATATTACCAGAGATGGTATCACCAGTGCGGAGGTTCAACCGGCGTATTTGGCTAGGGTAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggag                                                                                                       <  1:912462/150‑1 (MQ=255)
                                            agcagcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggagga                                                                                                     <  4:700069/151‑1 (MQ=255)
                                            agcagcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACggagg                                                                                                      >  3:1117267/1‑150 (MQ=255)
                                              cagcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGaaa                                                                                                   <  1:652678/151‑1 (MQ=255)
                                               agcGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAAc                                                                                                  >  7:231571/1‑151 (MQ=255)
                                                 cgcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCa                                                                                                <  1:1003629/151‑1 (MQ=255)
                                                  gcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCaa                                                                                               <  3:1656718/151‑1 (MQ=255)
                                                  gcAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCa                                                                                                >  5:1241842/1‑150 (MQ=255)
                                                    aaaaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAGACCAAat                                                                                             >  8:696819/1‑151 (MQ=255)
                                                      aaTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatcc                                                                                           >  1:387963/1‑151 (MQ=255)
                                                         aGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatccatc                                                                                        >  7:1164598/1‑151 (MQ=255)
                                                          gCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatccatcc                                                                                       <  4:901031/151‑1 (MQ=255)
                                                           cGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAatccatccn                                                                                      <  3:561841/151‑2 (MQ=255)
                                                               cACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCa                                                                                   <  5:1514452/150‑1 (MQ=255)
                                                                aCCTTCTTTCGGCTGTCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGTTGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGATCGGATGAAACCAAATCCATCCTGCAat                                                                                 <  2:1704376/151‑1 (MQ=255)
                                                                aCCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAACCCAAATTCATCCTGCAat                                                                                 >  3:1811246/1‑151 (MQ=255)
                                                                aCCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAat                                                                                 <  4:959155/151‑1 (MQ=255)
                                                                 ccTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAat                                                                                 >  2:1333462/1‑150 (MQ=255)
                                                                    tcttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATAtc                                                                              <  8:182134/150‑1 (MQ=255)
                                                                      tttCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTcc                                                                           >  5:1657418/1‑151 (MQ=255)
                                                                        tCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCa                                                                          <  2:865795/150‑1 (MQ=255)
                                                                           gcgggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGtcc                                                                      >  2:775770/1‑149 (MQ=255)
                                                                            cgggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGCAGCTGTCTGCGGAACGGAGGAAAACAAATCCATCCTGCAATATCTCCAGTAcc                                                                     >  2:352138/1‑150 (MQ=255)
                                                                              ggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGcc                                                                   >  7:1007139/1‑151 (MQ=255)
                                                                              ggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGcc                                                                   >  1:136305/1‑151 (MQ=255)
                                                                                   aTCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGCAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGCCCGGCGCGGTAGGAGCTGTCTCCGGAACGGAGGCACCCAAATCCATCCTGCAATATCTCCAGTACGCCATCAc                                                              >  1:1163941/1‑151 (MQ=255)
                                                                                     cTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCAcc                                                             >  7:535797/1‑150 (MQ=255)
                                                                                       tACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCaa                                                           <  5:533049/150‑1 (MQ=255)
                                                                                       tACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCa                                                            <  3:1015212/149‑1 (MQ=255)
                                                                                        aCCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCaa                                                           >  8:1399354/1‑149 (MQ=255)
                                                                                        aCCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAg                                                         <  8:75224/151‑1 (MQ=255)
                                                                                        aCCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAg                                                         <  7:767356/151‑1 (MQ=255)
                                                                                        aCCAGAGATGGTATCACCAGTGCAGAGGTTGAAACGGCGGATTTGGCTAGGGGAGACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAg                                                         <  6:1510803/151‑1 (MQ=255)
                                                                                         ccAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGa                                                        <  6:1771135/151‑1 (MQ=255)
                                                                                          cAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGat                                                       <  3:1519763/151‑1 (MQ=255)
                                                                                             agaTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATc                                                    >  6:97353/1‑151 (MQ=255)
                                                                                             agaTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATc                                                    <  2:958302/151‑1 (MQ=255)
                                                                                             agaTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACACCATCACCAAAGATATc                                                    >  7:1562215/1‑151 (MQ=255)
                                                                                               aTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCtt                                                  >  8:722901/1‑151 (MQ=255)
                                                                                                  gTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGc                                               >  1:231140/1‑151 (MQ=255)
                                                                                                      cACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGACCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGACATCACCAAAGATATCTTCGACAct                                           >  6:281359/1‑151 (MQ=255)
                                                                                                         cAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAACCCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCtt                                        >  1:1812260/1‑151 (MQ=255)
                                                                                                         cAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCtt                                        >  3:910184/1‑151 (MQ=255)
                                                                                                         cAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCtt                                        <  8:422266/151‑1 (MQ=255)
                                                                                                         cAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCtt                                        >  1:782689/1‑151 (MQ=255)
                                                                                                                gAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGtgc                                 <  2:1046182/151‑1 (MQ=255)
                                                                                                                      gAAACGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTc                           >  4:1137862/1‑151 (MQ=255)
                                                                                                                       aaaCGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCa                          <  4:476991/151‑1 (MQ=255)
                                                                                                                        aaCGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTc                           >  6:712482/1‑149 (MQ=255)
                                                                                                                        aaCGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAg                         <  1:741412/151‑1 (MQ=255)
                                                                                                                        aaCGGCGGATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAg                         <  8:156609/151‑1 (MQ=255)
                                                                                                                          cggcggATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGa                       <  3:983785/151‑1 (MQ=255)
                                                                                                                           ggcggATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGAt                      <  6:751235/151‑1 (MQ=255)
                                                                                                                           ggcggATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGGCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGAt                      <  8:536948/151‑1 (MQ=255)
                                                                                                                              ggATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGc                   >  7:285026/1‑151 (MQ=255)
                                                                                                                              ggATTTGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCGTCAGGATGGc                   >  6:1161419/1‑151 (MQ=255)
                                                                                                                                   tGGCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCaaaa               <  3:120350/150‑1 (MQ=255)
                                                                                                                                     gCTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaata            <  1:1655644/151‑1 (MQ=255)
                                                                                                                                      cTAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaata            >  6:1461135/1‑150 (MQ=255)
                                                                                                                                       tAGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataat          >  3:1419527/1‑151 (MQ=255)
                                                                                                                                        aGGGGAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGGCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataa           <  1:445614/149‑1 (MQ=255)
                                                                                                                                            gAAACGCATATGTCATCAGGACCGGCGATGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCt     <  7:1739380/151‑1 (MQ=255)
                                                                                                                                            gAAACGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCt     <  2:18499/151‑1 (MQ=255)
                                                                                                                                             aaaCGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTg    >  1:1491689/1‑151 (MQ=255)
                                                                                                                                             aaaCGCAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTg    <  4:737047/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGCGGGCGTTTTCAGGTTTGTCGAAGTTAACTTCGTTAACTTTCAGCAGCGCAAAATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCT  >  DH1/4547669‑4547962

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: