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,8380TG98.8% 1913.0 / 1.0 496T406P (ACC→CCC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (261/229);  total (264/232)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 523 total aligned reads displayed.

CAAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAG  >  DH1/4546691‑4546984
                                                                                                                                                   |                                                                                                                                                  
 aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                >  3:1596328/1‑151 (MQ=255)
 aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                >  2:306794/1‑151 (MQ=255)
  aaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTg                                                                                                                                               >  4:1680908/1‑151 (MQ=255)
  aaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTg                                                                                                                                               <  8:447014/151‑1 (MQ=255)
   aTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGc                                                                                                                                              >  5:88021/1‑151 (MQ=255)
            aTTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTttatt                                                                                                                                     >  8:1158742/1‑151 (MQ=255)
              tttCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTttatta                                                                                                                                    >  5:697073/1‑150 (MQ=255)
                   tAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAgg                                                                                                                              >  2:1850661/1‑151 (MQ=255)
                    aaCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAgg                                                                                                                              >  7:1351939/1‑150 (MQ=255)
                    aaCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAgg                                                                                                                              >  1:763849/1‑150 (MQ=255)
                     aCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGaa                                                                                                                            <  1:33550/151‑1 (MQ=255)
                     aCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGa                                                                                                                             >  4:1830939/1‑150 (MQ=255)
                       aTCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAAtt                                                                                                                          >  8:1612412/1‑151 (MQ=255)
                        tCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTc                                                                                                                         >  1:1522372/1‑151 (MQ=255)
                         cGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTc                                                                                                                         <  6:1948520/150‑1 (MQ=255)
                          gCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCa                                                                                                                       <  7:385326/151‑1 (MQ=255)
                           ccAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAAATTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCAt                                                                                                                      <  6:695008/151‑1 (MQ=255)
                              tatcgtcATTAATATTATAGATATAAAAGCAAAATGCCACGTAAACAAGTGTCGTTTTTTTTATAAGACAAATTTATGATCGTTTAATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTg                                                                                                                    <  4:726463/144‑1 (MQ=255)
                              gcgcGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGc                                                                                                                   >  6:157253/1‑151 (MQ=255)
                                gcgGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCAt                                                                                                                 <  3:69905/151‑1 (MQ=255)
                                gcgGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTCTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCAt                                                                                                                 <  3:85854/151‑1 (MQ=255)
                                 cgGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGGTTATTAATGAGGAATTCCATTGCATc                                                                                                                >  6:147845/1‑151 (MQ=255)
                                    gaTTAAGATTACAGATATAAAAGCAAAACGCCACGTAAAGACGTGGCGTTTTTTGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTTGTCTTGGGCATTGCCAGTTTATTAATGTGGAATTCCATTGCATCGAt                                                                                                             <  1:1053949/150‑1 (MQ=255)
                                    cATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAt                                                                                                             >  6:407554/1‑151 (MQ=255)
                                    cATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGAGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAt                                                                                                             <  8:1805975/151‑1 (MQ=255)
                                     aTTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGTTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAtt                                                                                                            <  8:1416612/151‑1 (MQ=255)
                                        aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                          <  8:557382/150‑1 (MQ=255)
                                          gATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGcc                                                                                                       >  2:1192613/1‑151 (MQ=255)
                                           aTTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGccc                                                                                                      <  6:789760/151‑1 (MQ=255)
                                             tACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCAt                                                                                                    >  8:77035/1‑151 (MQ=255)
                                              aCAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATc                                                                                                   <  4:1282219/151‑1 (MQ=255)
                                                aGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCgg                                                                                                 <  1:509053/151‑1 (MQ=255)
                                                   tataAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTg                                                                                              >  2:1108848/1‑151 (MQ=255)
                                                    ataAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGa                                                                                             <  6:1047994/151‑1 (MQ=255)
                                                       aaaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATg                                                                                          <  5:395804/151‑1 (MQ=255)
                                                       aaaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATg                                                                                          <  3:394911/151‑1 (MQ=255)
                                                        aaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATg                                                                                          >  5:1778713/1‑150 (MQ=255)
                                                             aaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTg                                                                                    >  4:513148/1‑151 (MQ=255)
                                                               aCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcg                                                                                  >  6:273478/1‑151 (MQ=255)
                                                               aCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcg                                                                                  <  7:542490/151‑1 (MQ=255)
                                                                 gCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcgc                                                                                 >  4:1571375/1‑150 (MQ=255)
                                                                 gCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCa                                                                                <  1:586006/151‑1 (MQ=255)
                                                                  ccACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcgc                                                                                 >  6:1325664/1‑149 (MQ=255)
                                                                  ccACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcgc                                                                                 >  2:403471/1‑149 (MQ=255)
                                                                   cACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCTgg                                                                              >  7:1215471/1‑151 (MQ=255)
                                                                    aCGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAgg                                                                              >  2:1551103/1‑150 (MQ=255)
                                                                       tAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATcc                                                                          >  1:903182/1‑151 (MQ=255)
                                                                         aaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcac                                                                        <  3:1448982/151‑1 (MQ=255)
                                                                          acacGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCcaca                                                                       <  3:1276582/151‑1 (MQ=255)
                                                                            acGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAtt                                                                     >  7:1653279/1‑151 (MQ=255)
                                                                            acGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAtt                                                                     <  1:678535/151‑1 (MQ=255)
                                                                             cGTGGCGTTTTTGGCATAAGACAAATTTATGATCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCACATCGGGTGAATGATTTTGCGCAGGATCCACAttt                                                                    <  6:1799150/151‑1 (MQ=255)
                                                                             cGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAttt                                                                    <  5:1135781/151‑1 (MQ=255)
                                                                                    tttttGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                               >  4:238648/1‑149 (MQ=255)
                                                                                      tttGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGtt                                                           >  5:1696888/1‑151 (MQ=255)
                                                                                         ggCATAAGACTAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttctt                                                        <  8:1325072/151‑1 (MQ=255)
                                                                                         ggCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttctt                                                        <  2:1675253/151‑1 (MQ=255)
                                                                                         ggCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttctt                                                        <  7:685148/151‑1 (MQ=255)
                                                                                            aTAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCt                                                     <  7:1078662/151‑1 (MQ=255)
                                                                                             tAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTg                                                    >  7:1833222/1‑151 (MQ=255)
                                                                                              aaGACAAATTTTTGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGGCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTg                                                    <  5:18486/150‑1 (MQ=255)
                                                                                              aaGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCtaa                                                   <  1:1894133/151‑3 (MQ=255)
                                                                                               aGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAg                                                  <  7:19547/151‑1 (MQ=255)
                                                                                                gACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGt                                                 <  2:1330120/151‑1 (MQ=255)
                                                                                                 aCAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTc                                                >  5:1591163/1‑151 (MQ=255)
                                                                                                  cAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCg                                               >  8:327992/1‑151 (MQ=255)
                                                                                                  cAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCg                                               >  4:406833/1‑151 (MQ=255)
                                                                                                   aaaTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATAGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGt                                              <  3:1869933/151‑1 (MQ=255)
                                                                                                     aTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGa                                            <  4:1167403/151‑1 (MQ=255)
                                                                                                       ttATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagc                                          >  7:305187/1‑151 (MQ=255)
                                                                                                       ttATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGag                                           <  1:1612328/150‑1 (MQ=255)
                                                                                                        tATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagc                                          <  5:1403928/150‑1 (MQ=255)
                                                                                                          tGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagcagc                                       >  3:547691/1‑151 (MQ=255)
                                                                                                          tGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagcagc                                       >  2:1803240/1‑151 (MQ=255)
                                                                                                           gAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGct                                      <  5:1216835/151‑1 (MQ=255)
                                                                                                                tttCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCtcttct                                 >  4:939453/1‑151 (MQ=255)
                                                                                                                 ttCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATAGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCtcttct                                 <  2:614191/150‑1 (MQ=255)
                                                                                                                  tcatcaTTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTcttctt                                <  8:250833/150‑1 (MQ=255)
                                                                                                                  tcatcaTTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTcttctt                                <  5:248652/150‑1 (MQ=255)
                                                                                                                   catcatTTCTAATAAATCGTCATTGTTCTTGGGCATTGCCATTTAATTAATGAGGAATTCAATTTCATCGATTTCGCCCATGGGTTGAATGATTTAGCGCAGGATCAACATTTTCTGAAGTTATCCCTGAGTCGTGAGCAGCTCTTGttt                               <  1:1793453/150‑1 (MQ=255)
                                                                                                                   catcatTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCttt                               <  6:8548/150‑1 (MQ=255)
                                                                                                                   catcatTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTa                              <  8:1067033/151‑1 (MQ=255)
                                                                                                                      catTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACgg                           <  1:1277747/151‑1 (MQ=255)
                                                                                                                      catTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCAATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTGCCTGAGTCGTGAGCAGCTCTTCTTTACgg                           >  7:1487991/1‑151 (MQ=255)
                                                                                                                            gaagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTAcc                      >  7:1134937/1‑150 (MQ=255)
                                                                                                                            gaagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCa                     <  3:1751626/151‑1 (MQ=255)
                                                                                                                             aagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTAcc                      >  2:170971/1‑149 (MQ=255)
                                                                                                                                aaaTCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAAc                 <  3:101480/151‑1 (MQ=255)
                                                                                                                                aaaTCGTCATTGGTCTTGGGCATAGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCCGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCAGAGTCGTGCGCACCTCATATATACGGGTACAAGaa                  >  2:850004/1‑150 (MQ=255)
                                                                                                                                 aaTCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACg                >  7:842707/1‑151 (MQ=255)
                                                                                                                                    cGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGtt             >  1:218533/1‑151 (MQ=255)
                                                                                                                                    cGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGtt             <  1:170519/151‑1 (MQ=255)
                                                                                                                                     gTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTg            >  2:1425045/1‑151 (MQ=255)
                                                                                                                                     gTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTg            >  6:1224944/1‑151 (MQ=255)
                                                                                                                                       cATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCCGCTCTTCTTTACGGGTACCAGAACGGTTGt           >  4:1418024/1‑150 (MQ=255)
                                                                                                                                        aTTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGt           >  8:1687747/1‑149 (MQ=255)
                                                                                                                                            gtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCg      >  1:1361140/1‑150 (MQ=255)
                                                                                                                                            gtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGa     <  7:756493/151‑1 (MQ=255)
                                                                                                                                                tgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATa   <  7:961980/149‑1 (MQ=255)
                                                                                                                                                tgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAg  <  4:341220/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAG  >  DH1/4546691‑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: