Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,546,838 T→G T406P (ACC→CCC)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,546,8380TG96.2% 1764.7 / 22.4 472T406P (ACC→CCC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base G (223/230);  total (232/240)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 488 total aligned reads displayed.

AAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAG  >  DH1/4546692‑4546984
                                                                                                                                                  |                                                                                                                                                  
aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                <  3:908783/151‑1 (MQ=255)
aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                <  5:100007/151‑1 (MQ=255)
 aaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTg                                                                                                                                               >  8:666342/1‑151 (MQ=255)
 aaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTg                                                                                                                                               >  3:106961/1‑151 (MQ=255)
   ttCCAGGAATTTTCCTAACATCGCCAGCGCTGCAATAAGATTACAGATATAAAAGCAAAACGACACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAACTAGTCAGTGGTCTTGGGCATTGcc                                                                                                                                             >  2:1514846/1‑151 (MQ=255)
    tCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGcc                                                                                                                                             <  8:356031/150‑1 (MQ=255)
     ccAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGcc                                                                                                                                             >  3:588592/1‑149 (MQ=255)
      cAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGt                                                                                                                                          <  2:1162001/151‑1 (MQ=255)
       aGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGtt                                                                                                                                         >  2:907896/1‑151 (MQ=255)
       aGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGtt                                                                                                                                         >  2:1296187/1‑151 (MQ=255)
              ttCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTaa                                                                                                                                   <  7:1586559/150‑1 (MQ=255)
              ttCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAAt                                                                                                                                  >  3:955544/1‑151 (MQ=255)
               tCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAAt                                                                                                                                  <  3:403800/150‑1 (MQ=255)
                ccTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATg                                                                                                                                 <  8:841309/150‑1 (MQ=255)
                ccTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATg                                                                                                                                 >  5:921567/1‑150 (MQ=255)
                ccTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATg                                                                                                                                 <  6:560855/150‑1 (MQ=255)
                ccTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGa                                                                                                                                <  5:1313377/151‑1 (MQ=255)
                 cTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATg                                                                                                                                 >  2:988570/1‑149 (MQ=255)
                        cGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTcc                                                                                                                        <  7:620319/151‑1 (MQ=255)
                        cGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTcc                                                                                                                        <  3:255600/151‑1 (MQ=255)
                        cGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTc                                                                                                                         <  4:625274/150‑1 (MQ=255)
                         gCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTcc                                                                                                                        >  4:283976/1‑150 (MQ=255)
                         gCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACCTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTcc                                                                                                                        >  5:1199771/1‑150 (MQ=255)
                              cgcgGCATTAAGATTAAAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGc                                                                                                                   >  6:1419426/1‑150 (MQ=255)
                               gcgGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCAt                                                                                                                 <  5:292357/151‑1 (MQ=255)
                                     ttAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAttt                                                                                                           >  4:798916/1‑151 (MQ=255)
                                       aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCg                                                                                                         <  7:1069398/151‑1 (MQ=255)
                                       aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCg                                                                                                         >  6:382711/1‑151 (MQ=255)
                                        aGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGc                                                                                                        >  5:202053/1‑151 (MQ=255)
                                           ttACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGccc                                                                                                      >  8:1005288/1‑150 (MQ=255)
                                            tACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCa                                                                                                     <  2:31592/150‑1 (MQ=255)
                                            tACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCAt                                                                                                    >  7:1396852/1‑151 (MQ=255)
                                                gATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCggg                                                                                                >  1:360448/1‑151 (MQ=255)
                                                gATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCggg                                                                                                >  8:725609/1‑151 (MQ=255)
                                                 atatAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGt                                                                                               <  6:129602/151‑1 (MQ=255)
                                                  tataAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTg                                                                                              >  8:867960/1‑151 (MQ=255)
                                                   ataAAAGCAAAACGCCACGTAAACACGGGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTg                                                                                              >  1:178682/1‑150 (MQ=255)
                                                    taAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGa                                                                                             >  1:95113/1‑150 (MQ=255)
                                                     aaaaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAAt                                                                                           >  1:1352673/1‑151 (MQ=255)
                                                         gCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGGCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtt                                                                                       >  1:1184379/1‑151 (MQ=255)
                                                         gCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtt                                                                                       <  1:1002154/151‑1 (MQ=255)
                                                           aaaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtttt                                                                                     <  5:76760/151‑1 (MQ=255)
                                                           aaaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtttt                                                                                     >  7:272946/1‑151 (MQ=255)
                                                            aaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTGATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTg                                                                                    <  4:708370/151‑1 (MQ=255)
                                                            aaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAttttt                                                                                    >  3:322056/1‑150 (MQ=255)
                                                             aaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgc                                                                                   >  1:619729/1‑151 (MQ=255)
                                                              aCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcg                                                                                  <  5:1392362/151‑1 (MQ=255)
                                                               cGCCACGTAAACAGGTGTCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcgc                                                                                 <  7:1166426/151‑1 (MQ=255)
                                                                 ccACGTAAACACGTGGCGTTTTTGGCATAAGACAAAGTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAg                                                                               <  2:649377/151‑1 (MQ=255)
                                                                  cACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAgg                                                                              >  3:539711/1‑151 (MQ=255)
                                                                   aCGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGa                                                                             <  3:1521184/151‑1 (MQ=255)
                                                                       aaaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGAAGGATcc                                                                          >  7:1452638/1‑150 (MQ=255)
                                                                        aaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcac                                                                        >  3:1565134/1‑151 (MQ=255)
                                                                         acacGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcaca                                                                       <  3:864357/151‑1 (MQ=255)
                                                                          cacGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAt                                                                      <  4:1436658/151‑1 (MQ=255)
                                                                          cacGTCGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGACTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAt                                                                      >  1:1073996/1‑151 (MQ=255)
                                                                           acGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTAGTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAtt                                                                     >  4:379575/1‑151 (MQ=255)
                                                                            cGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAttt                                                                    <  2:1383251/151‑1 (MQ=255)
                                                                             gTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAtttt                                                                   <  6:201709/151‑1 (MQ=255)
                                                                               ggCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTc                                                                  <  4:1114556/150‑1 (MQ=255)
                                                                                 ggTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                               <  4:1285517/150‑1 (MQ=255)
                                                                                 cGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                               >  8:146058/1‑151 (MQ=255)
                                                                                     tttGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGt                                                            <  7:195507/150‑1 (MQ=255)
                                                                                        ggCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttct                                                         <  6:1252375/150‑1 (MQ=255)
                                                                                          cATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCTGGTGAATGATTTTGCGCAGGATCAACATTTTCTGCAGttcttc                                                       >  6:1512396/1‑150 (MQ=255)
                                                                                          cATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTcc                                                      >  4:1656485/1‑151 (MQ=255)
                                                                                             aaGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTg                                                    <  7:773538/150‑1 (MQ=255)
                                                                                             aaGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGa                                                   <  7:464601/151‑1 (MQ=255)
                                                                                               gACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGt                                                 <  1:1542166/151‑1 (MQ=255)
                                                                                               gACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGt                                                 <  2:744532/151‑1 (MQ=255)
                                                                                                aCAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTc                                                <  1:772570/151‑1 (MQ=255)
                                                                                                 aaaaTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCATTTTATTAATGAGGAATTCCATTGCATCGATTTCGCACATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCg                                               <  8:480990/150‑1 (MQ=255)
                                                                                                   aaTTTATGAGCGTTTCATCATTTCCAAGAAATCGTCATTGGACTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCATCCTGAGTCGTg                                             >  2:443611/1‑151 (MQ=255)
                                                                                                    aTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTg                                             >  8:57134/1‑150 (MQ=255)
                                                                                                    aTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGa                                            <  4:872464/151‑1 (MQ=255)
                                                                                                     tttATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCACATCCGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGACTCGTg                                             >  6:34773/1‑149 (MQ=255)
                                                                                                        aTGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagcac                                        >  7:889924/1‑150 (MQ=255)
                                                                                                        aTGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagca                                         <  4:124913/150‑1 (MQ=255)
                                                                                                         tGAGCGTTTCATCATTTCGAAGGAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagcagc                                       >  6:765793/1‑151 (MQ=255)
                                                                                                         tGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagcag                                        <  8:13219/150‑1 (MQ=255)
                                                                                                         tGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagcacc                                       >  8:1197847/1‑149 (MQ=255)
                                                                                                            gCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGGAGctct                                    >  3:1294917/1‑151 (MQ=255)
                                                                                                             cGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTTGCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCTGTTCTTCCTGAGTCGTGAGCAGCtctt                                   <  5:391596/151‑1 (MQ=255)
                                                                                                               tttCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCtcttct                                 >  2:188795/1‑151 (MQ=255)
                                                                                                                ttCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCtcttct                                 <  8:1680441/150‑1 (MQ=255)
                                                                                                                 tcatcaTTTCGAAGAAATCGTCATTTGTCTTGGGCATTGCGAGTTTATTTATGAGGAATTCCATTGCATCGATTTCGCGCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGGAGTTCTTCCTGAGTCGTGAGCAGCTcttctt                                <  5:1519899/150‑1 (MQ=255)
                                                                                                                   atcatTTCGTAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTAc                             >  4:48835/1‑151 (MQ=255)
                                                                                                                   atcatTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGAATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTATCTGAAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTa                              >  6:897105/1‑150 (MQ=255)
                                                                                                                     catTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACgg                           >  5:505246/1‑151 (MQ=255)
                                                                                                                     catTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTCACg                            <  5:1643323/150‑1 (MQ=255)
                                                                                                                      atTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACggg                          >  7:1009100/1‑151 (MQ=255)
                                                                                                                           gaagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTAcc                      >  1:464745/1‑150 (MQ=255)
                                                                                                                            aagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAg                    <  7:1317257/151‑1 (MQ=255)
                                                                                                                             agaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGa                   <  6:1064669/151‑1 (MQ=255)
                                                                                                                               aaaTCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAAc                 <  7:868909/151‑1 (MQ=255)
                                                                                                                               aaaTCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTAAGGGTACCAGAAc                 >  6:1154361/1‑151 (MQ=255)
                                                                                                                                   gttCATTGGTCTTTGGCATTGCCAGTTTATTAATTAGGAATTGCATTGCATCGATTTCGCCCATCGGTTTAATGATTTTGCGCAGGATCCAAATTTTCTGCAGTTCTTCCTGAGGCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGtt             <  3:683690/149‑1 (MQ=255)
                                                                                                                                          ggtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTc       <  6:599598/150‑1 (MQ=255)
                                                                                                                                               tgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATa   <  6:1645466/149‑1 (MQ=255)
                                                                                                                                               tgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAg  <  1:557035/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAG  >  DH1/4546692‑4546984

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: