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.4% 1344.6 / 18.9 358T406P (ACC→CCC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (9/4);  new base G (165/180);  total (174/184)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 390 total aligned reads displayed.

GCAAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAGC  >  DH1/4546690‑4546985
                                                                                                                                                    |                                                                                                                                                   
 cAAATTGCAGGAATTATCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAt                                                                                                                                                  <  8:1078830/151‑1 (MQ=255)
 cAAATTCCAGGAATTTTCCTAACATCTCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCa                                                                                                                                                   <  7:570083/150‑1 (MQ=255)
  aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                 >  8:116671/1‑151 (MQ=255)
  aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                 >  1:582950/1‑151 (MQ=255)
    aTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGc                                                                                                                                               <  5:592029/151‑1 (MQ=255)
         aGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGtt                                                                                                                                          <  3:134788/151‑1 (MQ=255)
         aGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGtt                                                                                                                                          <  7:1074120/151‑1 (MQ=255)
         aGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGtt                                                                                                                                          <  5:513242/151‑1 (MQ=255)
            aaTTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTttat                                                                                                                                       <  2:882861/151‑1 (MQ=255)
             aTTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAATAAATCGTCATTGGTCTTGGGCATTGCCAGTttatt                                                                                                                                      <  1:1081325/151‑1 (MQ=255)
             aTTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTttatt                                                                                                                                      <  4:69567/151‑1 (MQ=255)
              ttttCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTttatt                                                                                                                                      >  2:803084/1‑150 (MQ=255)
               tttCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTttatta                                                                                                                                     <  7:19226/150‑1 (MQ=255)
                ttCCTAACATCGCCAGAGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTaa                                                                                                                                    <  8:644172/150‑1 (MQ=255)
                  ccTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTTGGCATTGCCAGTTTATTAATg                                                                                                                                  <  5:1023505/150‑1 (MQ=255)
                    tAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAgg                                                                                                                               >  6:276686/1‑151 (MQ=255)
                    tAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTAGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAgg                                                                                                                               >  5:647310/1‑151 (MQ=255)
                    aaaCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGACAGTTTATTAATGAgg                                                                                                                               >  3:1189748/2‑151 (MQ=255)
                     aaCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGGCAGTTTATTAATGAgg                                                                                                                               >  5:1256082/1‑150 (MQ=255)
                     aaCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGa                                                                                                                              <  8:1150952/151‑1 (MQ=255)
                      aCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGa                                                                                                                              >  3:660836/1‑150 (MQ=255)
                           gCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTcc                                                                                                                         >  1:1030375/1‑150 (MQ=255)
                             cAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCAtt                                                                                                                      >  7:1271364/1‑151 (MQ=255)
                             cAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCAtt                                                                                                                      >  8:273765/1‑151 (MQ=255)
                              aGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGGATTGCCAGTTTATTAATGAGGAATTCCATTg                                                                                                                     <  2:173771/151‑1 (MQ=255)
                               gcgcGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGc                                                                                                                    <  5:1090905/151‑1 (MQ=255)
                                  cgGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATc                                                                                                                 <  7:308893/151‑1 (MQ=255)
                                    gCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGa                                                                                                               <  1:319876/151‑1 (MQ=255)
                                     cATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAt                                                                                                              >  2:1009420/1‑151 (MQ=255)
                                     cATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAt                                                                                                              <  4:427897/151‑1 (MQ=255)
                                     cATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAt                                                                                                              >  4:282277/1‑151 (MQ=255)
                                         aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCg                                                                                                          <  3:274210/151‑1 (MQ=255)
                                         aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCg                                                                                                          >  6:371144/1‑151 (MQ=255)
                                          aGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGc                                                                                                         >  4:362871/1‑151 (MQ=255)
                                            aTTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGccc                                                                                                       <  4:516131/151‑1 (MQ=255)
                                              tACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCAt                                                                                                     >  2:589684/1‑151 (MQ=255)
                                                   atatAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGt                                                                                                <  2:1087416/151‑1 (MQ=255)
                                                     ataAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCTGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGa                                                                                              <  5:92299/151‑1 (MQ=255)
                                                        aaaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATg                                                                                           <  6:553842/151‑1 (MQ=255)
                                                          aGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAt                                                                                         <  8:1286719/151‑1 (MQ=255)
                                                           gCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtt                                                                                        <  1:860917/151‑1 (MQ=255)
                                                            cAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAttt                                                                                       <  3:85703/151‑1 (MQ=255)
                                                             aaaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtttt                                                                                      <  7:371384/151‑1 (MQ=255)
                                                              aaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTg                                                                                     >  2:256448/1‑151 (MQ=255)
                                                               aaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgc                                                                                    >  6:441852/1‑151 (MQ=255)
                                                                     aCGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAgg                                                                               >  4:106690/1‑150 (MQ=255)
                                                                       gTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATc                                                                            >  5:1134534/1‑151 (MQ=255)
                                                                       gTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATc                                                                            <  1:1361338/151‑1 (MQ=255)
                                                                        tAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATcc                                                                           >  5:866813/1‑151 (MQ=255)
                                                                        tAAAAACGTGGCGTTTTTGGAATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATcc                                                                           >  8:468387/1‑151 (MQ=255)
                                                                         aaaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATcc                                                                           >  2:683974/1‑150 (MQ=255)
                                                                         aaaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATcc                                                                           >  3:1095687/1‑150 (MQ=255)
                                                                          aaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcac                                                                         >  6:1323673/1‑151 (MQ=255)
                                                                          aaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcac                                                                         >  2:135713/1‑151 (MQ=255)
                                                                           acacGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcac                                                                         >  4:400490/1‑150 (MQ=255)
                                                                           acacGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcac                                                                         >  2:1248512/1‑150 (MQ=255)
                                                                             acGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAtt                                                                      <  7:1150758/151‑1 (MQ=255)
                                                                                tGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTc                                                                   >  3:1246036/1‑151 (MQ=255)
                                                                                  gCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTg                                                                 <  7:747042/151‑1 (MQ=255)
                                                                                   cGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                                <  3:159247/151‑1 (MQ=255)
                                                                                   cGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                                <  6:834506/151‑1 (MQ=255)
                                                                                   cGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                                >  6:329543/1‑151 (MQ=255)
                                                                                    gTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                                >  3:472883/1‑150 (MQ=255)
                                                                                          ggCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttct                                                          <  8:549940/150‑1 (MQ=255)
                                                                                          ggCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttct                                                          <  4:1362475/150‑1 (MQ=255)
                                                                                             aTAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCt                                                      >  4:776470/1‑151 (MQ=255)
                                                                                             aTAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCt                                                      <  5:1369371/151‑1 (MQ=255)
                                                                                              tAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTg                                                     >  6:724526/1‑151 (MQ=255)
                                                                                              tAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTg                                                     >  7:1247346/1‑151 (MQ=255)
                                                                                              tAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTg                                                     >  1:914525/1‑151 (MQ=255)
                                                                                                   cAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCg                                                <  4:1241895/151‑1 (MQ=255)
                                                                                                    aaaTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGt                                               >  4:908925/1‑151 (MQ=255)
                                                                                                    aaaTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGt                                               <  7:676809/151‑1 (MQ=255)
                                                                                                      aTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTg                                              >  7:977228/1‑150 (MQ=255)
                                                                                                       tttATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTg                                              >  5:724271/1‑149 (MQ=255)
                                                                                                          aTGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagca                                          <  4:119569/150‑1 (MQ=255)
                                                                                                           tGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagcagc                                        >  3:1316463/1‑151 (MQ=255)
                                                                                                            gAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGct                                       <  8:733468/151‑1 (MQ=255)
                                                                                                            gAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGct                                       >  1:124065/1‑151 (MQ=255)
                                                                                                                    catcatTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCttt                                >  7:43959/1‑150 (MQ=255)
                                                                                                                    catcatTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCttt                                <  1:683350/150‑1 (MQ=255)
                                                                                                                      tcatTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACg                             >  5:951201/1‑151 (MQ=255)
                                                                                                                        atTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACggg                           >  7:236957/1‑151 (MQ=255)
                                                                                                                        atTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTAAggg                           >  5:1263980/1‑151 (MQ=255)
                                                                                                                               agaaATCGTCATTGGTGTTGGGGATTGCCAGTTTATTAATGATGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGa                    <  8:1370255/151‑1 (MQ=255)
                                                                                                                                gaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGaa                   <  2:310369/151‑1 (MQ=255)
                                                                                                                                gaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGa                    >  8:1043420/1‑150 (MQ=255)
                                                                                                                                 aaaTCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAAc                  >  4:82751/1‑151 (MQ=255)
                                                                                                                                 aaaTCGTAATTGTTCTTGGGGATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTGTCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCACTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAAc                  <  6:646841/151‑1 (MQ=255)
                                                                                                                                  aaTCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACg                 <  6:1195550/151‑1 (MQ=255)
                                                                                                                                      gTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGGAGCTCTTCTTTACGGGTACCAGAACGGTTg             <  4:471280/151‑1 (MQ=255)
                                                                                                                                       tCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGt            >  7:133614/1‑151 (MQ=255)
                                                                                                                                        cATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGt            <  8:31086/150‑1 (MQ=255)
                                                                                                                                           tggtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGt         <  5:400702/150‑1 (MQ=255)
                                                                                                                                           tggtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTc        >  1:745182/1‑151 (MQ=255)
                                                                                                                                           tggtcttgggcATTGCCAGTTTATTAAAGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTAACAGAACGGTTGTAGt         >  2:513780/1‑150 (MQ=255)
                                                                                                                                             gtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAg                                                   >  5:267527/1‑106 (MQ=255)
                                                                                                                                              tcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTGTTTACGGGTACCAGAACGGTTGTAGTCGAt     >  7:483272/1‑151 (MQ=255)
                                                                                                                                                ttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGGTGTAGTCGAt     >  6:149723/1‑149 (MQ=255)
                                                                                                                                                 tggtcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATa    <  1:244208/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAGC  >  DH1/4546690‑4546985

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