Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,547,872 C→T G61E (GGA→GAA)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,547,8720CT99.2% 1447.5 / ‑1.1 382G61E (GGA→GAA) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (182/196);  total (182/200)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 398 total aligned reads displayed.

ATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATCAGCTC  >  DH1/4547724‑4548017
                                                                                                                                                    |                                                                                                                                                 
aTAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAattca                                                                                                                                                 <  7:1075312/151‑1 (MQ=255)
 tAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAattcat                                                                                                                                                >  5:1343010/1‑151 (MQ=255)
  aGCGTTCACCTTCTTTCGGCTGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAattcatc                                                                                                                                               <  4:1241640/151‑1 (MQ=255)
  aGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAattcatc                                                                                                                                               <  4:459327/151‑1 (MQ=255)
  aGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAattcatc                                                                                                                                               >  4:647191/1‑151 (MQ=255)
    cGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCt                                                                                                                                             <  8:667200/151‑1 (MQ=255)
     gTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTg                                                                                                                                            <  4:302908/151‑1 (MQ=255)
       tCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGc                                                                                                                                           >  7:452857/1‑150 (MQ=255)
        cACCTTCTTTCGGCGGGCGAATCTTACCAGAGATTGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCaa                                                                                                                                         <  7:767460/151‑1 (MQ=255)
         aCCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAat                                                                                                                                        >  3:1181196/1‑151 (MQ=255)
         aCCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAat                                                                                                                                        <  5:1674657/151‑1 (MQ=255)
          ccTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAata                                                                                                                                       <  1:968503/151‑1 (MQ=255)
            ttcttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCCTCCTGCAATAtc                                                                                                                                     >  3:272370/1‑151 (MQ=255)
            ttcttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATAtc                                                                                                                                     >  6:1096875/1‑151 (MQ=255)
              cttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATAtctc                                                                                                                                   >  3:470530/1‑151 (MQ=255)
              cttTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCAGCCTGCAATAtctc                                                                                                                                   >  5:1543182/1‑151 (MQ=255)
                ttCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATAtct                                                                                                                                    <  7:1380231/148‑1 (MQ=255)
                 tCGGCGGGCGAATCTTACCAGAGATGGTATCATCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTcc                                                                                                                                  >  6:616972/1‑149 (MQ=255)
                 tCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCa                                                                                                                                 <  1:766051/150‑1 (MQ=255)
                 tCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCa                                                                                                                                 <  7:70279/150‑1 (MQ=255)
                  cggcggGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGt                                                                                                                               >  1:683691/1‑151 (MQ=255)
                   ggcgggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAAGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTa                                                                                                                              <  6:96896/151‑1 (MQ=255)
                    gcgggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTAc                                                                                                                             <  8:298138/151‑1 (MQ=255)
                     cgggccAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTAcc                                                                                                                            >  8:529728/1‑150 (MQ=255)
                      gggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGCCCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGCAACCAAATTCATCCTGCAATATCTCCAGTACGc                                                                                                                           >  7:1303789/1‑151 (MQ=255)
                      gggcgAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGAAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGc                                                                                                                           >  8:1586538/1‑151 (MQ=255)
                                   cAGAGATTGTATCAACAGTGCGGAGGTTGAAACGGCGGATTTGGGTAGGGGAAACGTTGATGTCATCAGGAGGGGCGTGGTAGGAGCTGTCTGCTGAACTGATGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCTAAGat                                                                                                              <  4:1392133/151‑1 (MQ=255)
                                   cAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCaa                                                                                                                  >  6:694292/1‑147 (MQ=255)
                                    agagATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGata                                                                                                             <  1:1628146/151‑1 (MQ=255)
                                        aTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCtt                                                                                                         <  2:1473045/151‑1 (MQ=255)
                                        aTGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCtt                                                                                                         <  6:30977/151‑1 (MQ=255)
                                          ggTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTc                                                                                                        <  3:1226565/150‑1 (MQ=255)
                                              tCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGcccc                                                                                                   >  8:1456190/1‑149 (MQ=255)
                                                 ccAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCActct                                                                                                >  2:1395035/1‑151 (MQ=255)
                                                  cAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCtt                                                                                               >  3:1322953/1‑151 (MQ=255)
                                                     tGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGc                                                                                            >  6:749177/1‑151 (MQ=255)
                                                       cGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGt                                                                                          >  7:1417960/1‑151 (MQ=255)
                                                        ggAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGtg                                                                                         >  2:1035165/1‑151 (MQ=255)
                                                         gAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGtgc                                                                                        <  5:279297/151‑1 (MQ=255)
                                                          aGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGtgct                                                                                       <  6:568272/151‑1 (MQ=255)
                                                          aGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGtgct                                                                                       <  5:755412/151‑1 (MQ=255)
                                                                     gcggATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATg                                                                            >  5:1188340/1‑151 (MQ=255)
                                                                     gcggATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATg                                                                            >  1:478363/1‑151 (MQ=255)
                                                                           ttGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCaaa                                                                       <  6:1363389/150‑1 (MQ=255)
                                                                             ggCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCaaaaa                                                                     <  1:1322420/150‑1 (MQ=255)
                                                                               cTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaata                                                                   >  8:423947/1‑150 (MQ=255)
                                                                                tAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataat                                                                 >  3:155987/1‑151 (MQ=255)
                                                                                 aGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataa                                                                  <  7:1246659/149‑1 (MQ=255)
                                                                                 aGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataa                                                                  <  4:276366/149‑1 (MQ=255)
                                                                                 aGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataa                                                                  <  2:169319/149‑1 (MQ=255)
                                                                                 aGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGGCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAaataat                                                                 <  1:53577/150‑1 (MQ=255)
                                                                                  ggggACACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGt                                                               <  2:48066/151‑1 (MQ=255)
                                                                                   gggAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTc                                                              >  8:772993/1‑151 (MQ=255)
                                                                                     gAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCt                                                            <  2:1619686/151‑1 (MQ=255)
                                                                                      aaaCGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTg                                                           >  3:1103536/1‑151 (MQ=255)
                                                                                       aaCGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCTGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGc                                                          <  1:854367/151‑1 (MQ=255)
                                                                                       aaCGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGc                                                          <  7:974599/151‑1 (MQ=255)
                                                                                          gTAGATGTCATCAGTACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGTAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTa                                                       <  4:1172064/151‑1 (MQ=255)
                                                                                          gTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCtt                                                        >  3:1023959/1‑150 (MQ=255)
                                                                                          gTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTa                                                       <  4:586616/151‑1 (MQ=255)
                                                                                            aGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTAc                                                      <  2:1271494/150‑1 (MQ=255)
                                                                                             gATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGc                                                    >  5:437464/1‑151 (MQ=255)
                                                                                                gTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATa                                                 <  4:156571/151‑1 (MQ=255)
                                                                                                gTCATCAGGACCGGCGAGCTAGGAGCTGTCTGCGGACCGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGGGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATa                                                 <  2:370337/151‑1 (MQ=255)
                                                                                                 tcatcaGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATAc                                                >  6:464583/1‑151 (MQ=255)
                                                                                                   atcaGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGa                                              <  3:1002476/151‑1 (MQ=255)
                                                                                                    tcaGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAATAATAATGTCCTGCTTACGCATACg                                               >  8:1646134/1‑149 (MQ=255)
                                                                                                        gACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAg                                         <  1:1532130/151‑1 (MQ=255)
                                                                                                           cGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGtt                                      >  7:1661095/1‑151 (MQ=255)
                                                                                                           cGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGtt                                      <  2:936906/151‑1 (MQ=255)
                                                                                                            ggCGAGGTAGGAGCTGTCTGCGGAACGGATGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGttt                                     <  2:12078/151‑1 (MQ=255)
                                                                                                            ggCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAGTTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGttt                                     >  8:1290231/1‑151 (MQ=255)
                                                                                                            ggCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGttt                                     <  2:666197/151‑1 (MQ=255)
                                                                                                               gAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTcc                                  <  6:1246419/151‑1 (MQ=255)
                                                                                                               gAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTcc                                  >  2:818333/1‑151 (MQ=255)
                                                                                                               gAGGTAGGAGCTGTCTGCCGAACGGAGCAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCCCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTcc                                  <  2:300967/151‑1 (MQ=255)
                                                                                                                aGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCa                                 <  3:729821/151‑1 (MQ=255)
                                                                                                                 ggTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTcc                                  >  7:607782/1‑149 (MQ=255)
                                                                                                                 ggTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTcc                                  >  6:902372/1‑149 (MQ=255)
                                                                                                                 ggTAGGAGCTGTCTGCGGAACCGAGGAAACCAAATTCATCCTGCAATATCTCCCGTACGCCATCACCAAAGATATCTTCGCCCCTCTTTGCGTGCTGCTTCAGGATGCCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGtgttgc                                  >  8:1247079/1‑147 (MQ=255)
                                                                                                                      gAGCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGcccc                            >  4:1117658/1‑150 (MQ=255)
                                                                                                                       aGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATCGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCat                          <  2:233269/151‑1 (MQ=255)
                                                                                                                        gCTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCata                         <  6:1537429/151‑1 (MQ=255)
                                                                                                                         cTGTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCatat                        >  4:1071944/1‑151 (MQ=255)
                                                                                                                           gTCTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATAttt                      >  5:1516362/1‑151 (MQ=255)
                                                                                                                             cTGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTc                    >  8:1112816/1‑151 (MQ=255)
                                                                                                                              tGCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCg                   >  7:1588352/1‑151 (MQ=255)
                                                                                                                               gCGGAACGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTc                    <  1:1392516/149‑1 (MQ=255)
                                                                                                                                     cGGAGGAAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGagag            <  1:218096/151‑1 (MQ=255)
                                                                                                                                      tgatgaAAACATATTCATCCTTCAATATCTCCAGTACTCCATCACTAAAGATATCTTCTGCTCTCTTTGAGTGCTGCTTCAGGATGGTAAAAATATTGTCCTGCTTACGCATACGAGCCATTTTTTCGAGCCCCATATTTTCGCGGagag            <  8:1044427/146‑1 (MQ=255)
                                                                                                                                      ggaggaAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGt           >  1:1139780/1‑151 (MQ=255)
                                                                                                                                        aggaAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTg          >  1:1187556/1‑150 (MQ=255)
                                                                                                                                        acgaAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGa         <  6:1651007/149‑1 (MQ=255)
                                                                                                                                          gaAACCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATc       >  4:815883/1‑151 (MQ=255)
                                                                                                                                           aaaCCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATCa      <  2:1117280/151‑1 (MQ=255)
                                                                                                                                            aaCCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATc       >  2:880940/1‑149 (MQ=255)
                                                                                                                                            aaCCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATCAg     <  1:1002867/151‑1 (MQ=255)
                                                                                                                                             aCCAAATTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATCAGc    >  3:953534/1‑151 (MQ=255)
                                                                                                                                                aaaTTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATCAGCTc  >  8:1199865/1‑150 (MQ=255)
                                                                                                                                                aaaTTCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGAGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATCAGCTc  >  5:816336/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATAGCGTTCACCTTCTTTCGGCGGGCGAATCTTACCAGAGATGGTATCACCAGTGCGGAGGTTGAAACGGCGGATTTGGCTAGGGGAAACGTAGATGTCATCAGGACCGGCGAGGTAGGAGCTGTCTGCGGAACGGAGGAAACCAAATCCATCCTGCAATATCTCCAGTACGCCATCACCAAAGATATCTTCGCCACTCTTTGCGTGCTGCTTCAGGATGGCAAAAATAATGTCCTGCTTACGCATACGAGCCAGGTTTTCCAGCCCCATATTTTCGCCGAGAGTGATCAGCTC  >  DH1/4547724‑4548017

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: ‑