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,8380TG94.0% 1029.5 / 33.8 284T406P (ACC→CCC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (11/6);  new base G (123/144);  total (134/150)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

Only 100 of 310 total aligned reads displayed.

GCAAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAGC  >  DH1/4546690‑4546985
                                                                                                                                                    |                                                                                                                                                   
 cAAATTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAt                                                                                                                                                  <  6:250144/151‑1 (MQ=255)
  aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACATATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                 <  4:594845/151‑1 (MQ=255)
  aaaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCAtt                                                                                                                                                 >  4:835875/1‑151 (MQ=255)
   aaTTCCAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTg                                                                                                                                                >  3:716183/1‑151 (MQ=255)
       ccATGAATTTTCCTAACATCGCCATCGAGGCATTAAGATTACAGATTTAAAAGGAAATCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCa                                                                                                                                             <  2:43989/150‑1 (MQ=255)
       ccAGGAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCa                                                                                                                                             <  5:695226/150‑1 (MQ=255)
           gAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGttt                                                                                                                                         >  1:738149/1‑150 (MQ=255)
           gAATTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGttt                                                                                                                                         >  5:2526/1‑150 (MQ=255)
             aTTTTCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTATTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGCTCTTCGGCATTGCCAGTttatt                                                                                                                                      >  6:842164/1‑151 (MQ=255)
              ttttCCTAACATCGCCATCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACTTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCTAAGAAATCGTCATTTGTCTTGGGCAATGACAGTttatt                                                                                                                                      >  8:809422/1‑150 (MQ=255)
                ttCCTAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAAt                                                                                                                                   >  7:142716/1‑151 (MQ=255)
                    tAACATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAgg                                                                                                                               >  1:254588/1‑151 (MQ=255)
                      aCATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGa                                                                                                                              >  5:395016/1‑150 (MQ=255)
                       cATCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAAt                                                                                                                            <  5:544761/151‑1 (MQ=255)
                        aTCGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAAtt                                                                                                                           >  6:372797/1‑151 (MQ=255)
                          cGCCAGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTcc                                                                                                                         >  3:30707/1‑151 (MQ=255)
                             cAGCGCGGCATTAAGATTACAGATATAATAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCAtt                                                                                                                      <  2:569897/151‑1 (MQ=255)
                              aGCGCGGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTg                                                                                                                     >  2:337910/1‑151 (MQ=255)
                                 gcgGCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTCCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCAt                                                                                                                  <  1:348549/151‑1 (MQ=255)
                                   ggCATTAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCg                                                                                                                <  4:516565/151‑1 (MQ=255)
                                        tAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAttt                                                                                                            <  4:16469/150‑1 (MQ=255)
                                        tAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGAttt                                                                                                            <  7:431443/150‑1 (MQ=255)
                                        tAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                           >  5:881564/1‑151 (MQ=255)
                                        tAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                           >  6:680265/1‑151 (MQ=255)
                                        tAAGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                           >  8:155403/1‑151 (MQ=255)
                                         aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                           <  7:257667/150‑1 (MQ=255)
                                         aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                           <  4:302250/150‑1 (MQ=255)
                                         aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                           <  2:208906/150‑1 (MQ=255)
                                         aaGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTc                                                                                                           <  7:640651/150‑1 (MQ=255)
                                          aGATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGc                                                                                                         <  4:1041336/151‑1 (MQ=255)
                                           gATTACAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGcc                                                                                                        >  3:115122/1‑151 (MQ=255)
                                                cAGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCg                                                                                                   >  4:916782/1‑151 (MQ=255)
                                                 aGATATAAAATCAAAACGCCACGTAAAGACGTTGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCgg                                                                                                  <  6:318235/151‑1 (MQ=255)
                                                 aGATATAAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCgg                                                                                                  <  3:204623/151‑1 (MQ=255)
                                                 aGATATAAAAGCAAAACGCCACGT‑‑ACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGt                                                                                                >  8:1052757/1‑151 (MQ=255)
                                                     ataAAAGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTg                                                                                               <  7:286518/150‑1 (MQ=255)
                                                       aaaaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAAt                                                                                            >  8:210171/1‑151 (MQ=255)
                                                       aaaaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAAt                                                                                            <  1:504022/151‑1 (MQ=255)
                                                        aaaGCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATg                                                                                           >  7:304431/1‑151 (MQ=255)
                                                           gCAAAACGCCACGTAAACACGTGGCTTTTTTTGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtt                                                                                        <  6:537879/151‑1 (MQ=255)
                                                           gCAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtt                                                                                        <  1:926836/151‑1 (MQ=255)
                                                            cAAAACGCCACGTAAACACGTTGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAttt                                                                                       <  8:270490/151‑1 (MQ=255)
                                                            cAAAACGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAttt                                                                                       <  2:810844/151‑1 (MQ=255)
                                                             aaaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtttt                                                                                      >  5:1063267/1‑151 (MQ=255)
                                                             aaaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGAtttt                                                                                      <  8:869458/151‑1 (MQ=255)
                                                              aaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTGATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGGTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATGGGGTGAATGATTTTg                                                                                     >  5:979156/1‑151 (MQ=255)
                                                              aaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTTATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTg                                                                                     <  7:979259/151‑1 (MQ=255)
                                                              aaaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTg                                                                                     >  2:982346/1‑151 (MQ=255)
                                                               aaCGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgc                                                                                    <  7:857664/151‑1 (MQ=255)
                                                                 cGCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTgcgc                                                                                  >  7:169705/1‑151 (MQ=255)
                                                                  gCCACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCa                                                                                 <  7:927908/151‑1 (MQ=255)
                                                                   ccACGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGGCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAg                                                                                <  7:793389/151‑1 (MQ=255)
                                                                     aCGTAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAgg                                                                               >  2:536599/1‑150 (MQ=255)
                                                                        tAAACACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATcc                                                                           >  5:280832/1‑151 (MQ=255)
                                                                         aaaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATcc                                                                           >  6:996389/1‑150 (MQ=255)
                                                                         aaaCACGTGGCGTTTTTGGCATAAGACAAATTTATGAGAGTTTCATGATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTGCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCca                                                                          <  8:744922/151‑1 (MQ=255)
                                                                                tGGCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACAtttt                                                                    >  3:564294/1‑150 (MQ=255)
                                                                                 ggCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCt                                                                  <  2:786496/151‑1 (MQ=255)
                                                                                 ggCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCt                                                                  >  8:562433/1‑151 (MQ=255)
                                                                                 ggCGTTTTTGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAG‑TTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTg                                                                 <  1:604991/151‑1 (MQ=255)
                                                                                   cGTTTTTGGCATAAGACAAATTTATGAGCGATTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCTGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGc                                                                <  3:412155/151‑1 (MQ=255)
                                                                                      ttttGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGTAATTCCATTGCATCGATTTCGCCCATCGGGTTAATGATTTTGCGCAGGATCCACATTTTCTGCAg                                                              <  2:1039205/150‑1 (MQ=255)
                                                                                        ttGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttc                                                           >  5:826044/1‑151 (MQ=255)
                                                                                        ttGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGtt                                                            <  2:397126/150‑1 (MQ=255)
                                                                                         tGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttct                                                          >  2:651714/1‑151 (MQ=255)
                                                                                         tGGCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttc                                                           <  8:601637/150‑1 (MQ=255)
                                                                                          ggCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttctt                                                         <  7:322808/151‑1 (MQ=255)
                                                                                          ggCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttct                                                          <  4:649710/150‑1 (MQ=255)
                                                                                           gCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttcttc                                                        >  6:892247/1‑151 (MQ=255)
                                                                                           gCATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGttctt                                                         <  7:336757/150‑1 (MQ=255)
                                                                                            cATAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTcc                                                       >  7:509500/1‑151 (MQ=255)
                                                                                             aTAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCt                                                      >  2:936702/1‑151 (MQ=255)
                                                                                              tAAGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCTGGTGaa                                                                                             <  4:43356/111‑1 (MQ=255)
                                                                                               aaGACAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTg                                                     <  5:670208/150‑1 (MQ=255)
                                                                                                aGACATATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAg                                                   <  4:255308/151‑1 (MQ=255)
                                                                                                  aCAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGTCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTc                                                 >  1:403285/1‑151 (MQ=255)
                                                                                                  aCAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTc                                                 <  7:205163/151‑1 (MQ=255)
                                                                                                  aCAAATTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTc                                                 >  1:879689/1‑151 (MQ=255)
                                                                                                    aaaTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGt                                               >  5:133591/1‑151 (MQ=255)
                                                                                                     aaTTTATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTg                                              >  3:420750/1‑151 (MQ=255)
                                                                                                         tATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagc                                           <  5:332576/150‑1 (MQ=255)
                                                                                                         tATGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagc                                           >  3:933990/1‑150 (MQ=255)
                                                                                                          aTGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAGTTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagca                                          <  3:674133/150‑1 (MQ=255)
                                                                                                          aTGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagca                                          <  2:254322/150‑1 (MQ=255)
                                                                                                          aTGAGCGTTTCATCATTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGagca                                          <  4:762199/150‑1 (MQ=255)
                                                                                                                       catTTCGAAGAAATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACgg                            <  7:210373/151‑1 (MQ=255)
                                                                                                                             gaagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTAcc                       >  8:638819/1‑150 (MQ=255)
                                                                                                                             gaagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCa                      <  3:83263/151‑1 (MQ=255)
                                                                                                                              aagaaATCGTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCac                     >  3:659568/1‑150 (MQ=255)
                                                                                                                               agaaATCGTCATTGGTCTTGGGCATTGCCGGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTAcc                       >  3:234517/1‑148 (MQ=255)
                                                                                                                                 aaaTCGTCATTGGTCTTGGGCAGTGCCAGTTTATTAATGAGGAATGCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAAc                  <  6:462457/151‑1 (MQ=255)
                                                                                                                                      gTCATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTg             <  7:670634/151‑1 (MQ=255)
                                                                                                                                        cATTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGt            <  8:419263/150‑1 (MQ=255)
                                                                                                                                         aTTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAg          <  2:185653/151‑1 (MQ=255)
                                                                                                                                         aTTGGTCTTGGGCATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAg          <  6:191728/151‑1 (MQ=255)
                                                                                                                                           tggtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGt         <  1:688822/150‑1 (MQ=255)
                                                                                                                                             gtcttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGa      <  5:933707/151‑1 (MQ=255)
                                                                                                                                               gttggtcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATa    <  4:361841/150‑1 (MQ=255)
                                                                                                                                               cttgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATATCGCCCATCGGGTCAATGATTTTGCGCAGGATCCACATTTTCTGCACATCTTCCTGAGTCGAGAGCACCTGATGTTTATCGCTAGCAGATAGGTTGtttttcat     >  7:995281/1‑143 (MQ=255)
                                                                                                                                                 tgggcATTGCCAGTTTATTAATGAGGAATTCCATTGCATCGATTTCGCCCATCGGGTGAATGATTTTGCGCAGGATCCACATTTTCTGCAGTTCTTCCTGAGTCGTGAGCAGCTCTTCTTTACGGGTACCAGAACGGTTGTAGTCGATAg   <  1:1064775/150‑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: ‑