Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 3,701,787 Δ1 bp coding (2107/2244 nt) fhuA ← ferrichrome‑iron receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,701,7860C.94.4% 1807.7 / 117.9 426coding (2108/2244 nt)fhuAferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base C (13/11);  new base . (206/196);  total (219/207)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 443 total aligned reads displayed.

CCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCAAAGAA  >  DH1/3701636‑3701931
                                                                                                                                                      |                                                                                                                                                 
ccccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGc                                                                                                                                                   >  3:1493455/1‑151 (MQ=255)
ccccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTAGTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGc                                                                                                                                                   >  3:1216001/1‑151 (MQ=255)
  ccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGcg                                                                                                                                                  >  5:1150394/1‑149 (MQ=255)
     aaaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGcgcctc                                                                                                                                              >  3:1406112/1‑146 (MQ=255)
     aaaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTc                                                                                                                                             >  5:1640048/1‑151 (MQ=255)
       agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTcgc                                                                                                                                           >  5:1212016/1‑151 (MQ=255)
        gagaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCCCGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTcgcg                                                                                                                                          >  5:982086/1‑151 (MQ=255)
        gagaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTcgcg                                                                                                                                          <  1:1541707/151‑1 (MQ=255)
         agaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCgcgc                                                                                                                                         <  6:564685/151‑1 (MQ=255)
             aTTAGAAACGGAAGGTTGCGGTTTCAACGTACTGTAGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGa                                                                                                                                     <  8:550643/151‑1 (MQ=255)
             aTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGcc                                                                                                                                        >  7:1440369/1‑148 (MQ=255)
               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGAt                                                                                                                                    <  2:1609671/150‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCa                                                                                                                                  <  7:276680/151‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCa                                                                                                                                  <  1:1184395/151‑1 (MQ=255)
                   aaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATaa                                                                                                                               <  7:1534276/151‑1 (MQ=255)
                    aCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAAc                                                                                                                              <  1:548827/151‑1 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACg                                                                                                                             >  5:1542954/1‑151 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACg                                                                                                                             >  4:1623379/1‑151 (MQ=255)
                      ggAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGt                                                                                                                            >  6:381327/1‑151 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGt                                                                                                                            >  7:1183451/1‑150 (MQ=255)
                          ggttgcggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACt                                                                                                                         >  3:410943/1‑150 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCAt                                                                                                                  >  5:1343781/1‑151 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGAGTTAAACCAGCTCGCGACGTATTCACGATCGAACAGGTTGATAACATGCAGCGCCACGTTCCAGCCAGCCATGCCGACACGCGCCAGATCCTAACGTACAAAccc                                                                                                                    >  2:80327/1‑148 (MQ=255)
                                   tgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                                                >  3:1803121/1‑150 (MQ=255)
                                       aCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGAc                                                                                                           <  1:1157502/151‑1 (MQ=255)
                                       aCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGAc                                                                                                           >  1:1392567/1‑151 (MQ=255)
                                       aCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGAAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGAc                                                                                                           >  8:247200/1‑151 (MQ=255)
                                         gACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCg                                                                                                         >  5:871097/1‑151 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                        <  3:355115/151‑1 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCACAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACCATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                        <  8:1144035/151‑1 (MQ=255)
                                           ccTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                        >  2:1518335/1‑150 (MQ=255)
                                              gACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACTATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                                      <  7:1265065/149‑1 (MQ=255)
                                               aCGTTCTGCGCCCCAGAAGCAGCCATAAGTTTTAAAGCAGCTGGCGACTTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATaa                                                                                                    <  7:693841/151‑1 (MQ=255)
                                                  ttCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTc                                                                                                >  1:1104017/1‑151 (MQ=255)
                                                     tGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTccc                                                                                              >  5:311839/1‑150 (MQ=255)
                                                     tGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGGATAACTTccc                                                                                              >  8:1454540/1‑150 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCCTAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCAc                                                                                            <  6:1418655/151‑1 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCATAAGTGTTAAAGGAGCTGGCGACGTATTCACGATCGAACTGGTTGTTAACATGCAGCGCCAGGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCAc                                                                                            <  4:1753834/151‑1 (MQ=255)
                                                         ccccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACttt                                                                                         >  8:715932/1‑151 (MQ=255)
                                                          cccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTa                                                                                        <  3:1659913/151‑1 (MQ=255)
                                                                agcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAgg                                                                                    >  6:914769/1‑149 (MQ=255)
                                                                 gcagcCATAAGTTTTAAAGAAGCTGGCGACGTATTCCCGATCGAACAGGTTGTTAACATGCCGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGTTCATAACGTACTAACGAATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  7:1798158/151‑1 (MQ=255)
                                                                 gcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  4:1491369/151‑1 (MQ=255)
                                                                 gcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  2:596651/151‑1 (MQ=255)
                                                                  cagcCATAAGTGTTAAAGCAGCTGGCGACTTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  5:44838/151‑1 (MQ=255)
                                                                  cagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                                                >  8:1637746/1‑151 (MQ=255)
                                                                   agcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                                                >  5:141345/1‑150 (MQ=255)
                                                                     ccATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGc                                                                             <  1:61424/151‑1 (MQ=255)
                                                                     caaTAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCAGAACGTACTAACGCATCCACGCCAGTATAACTTCCCACTTTAAAGGAGTTatc                                                                             >  7:1610083/3‑149 (MQ=255)
                                                                         aaGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCgg                                                                          <  1:307833/151‑1 (MQ=255)
                                                                         aaGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCgg                                                                          <  6:201235/150‑1 (MQ=255)
                                                                               tAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAcc                                                                    <  8:812356/150‑1 (MQ=255)
                                                                                aaaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCa                                                                   <  2:171311/151‑1 (MQ=255)
                                                                                aaaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                  <  5:1743588/151‑1 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                  >  3:99315/1‑150 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATa                                                                 <  4:1803689/151‑1 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTAAAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                  >  6:1101161/1‑150 (MQ=255)
                                                                                  agcagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATaa                                                                <  4:348062/151‑1 (MQ=255)
                                                                                   gcagcTGGCGACGTATTCCCGCTCGAACAGGTTGTTACCATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCCTCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGACGTCTCACCATCAc                                                               >  6:1283821/1‑151 (MQ=255)
                                                                                    cagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACt                                                              >  4:1049704/1‑151 (MQ=255)
                                                                                      gcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            <  3:508650/151‑1 (MQ=255)
                                                                                       cTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  2:251636/1‑150 (MQ=255)
                                                                                       cTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  7:1119380/1‑150 (MQ=255)
                                                                                        tGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  5:1847820/1‑149 (MQ=255)
                                                                                        tGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  6:3341/1‑149 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                       <  7:1371141/151‑1 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                       <  3:815385/151‑1 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                                        >  2:890852/1‑149 (MQ=255)
                                                                                             aCGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGt                                                     <  5:712379/151‑1 (MQ=255)
                                                                                                 aTTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtat                                                   >  7:322810/1‑150 (MQ=255)
                                                                                                  ttCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAAc                                                >  6:784588/1‑151 (MQ=255)
                                                                                                   tCACGATCTAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAAc                                                <  1:1703242/150‑1 (MQ=255)
                                                                                                    cACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACgc                                              >  2:454178/1‑150 (MQ=255)
                                                                                                     aCGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGac                                             >  4:560706/1‑151 (MQ=255)
                                                                                                      cGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGacc                                            >  4:239553/1‑151 (MQ=255)
                                                                                                        aTCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGaccac                                          >  1:809774/1‑151 (MQ=255)
                                                                                                         tCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGaccacc                                         >  5:633371/1‑151 (MQ=255)
                                                                                                          cGACCAGGTTGTTAACATGCAGCGCCACGTTGCAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAATTGTAGCCAGTATAACGACCACCg                                        <  4:1269171/151‑1 (MQ=255)
                                                                                                          cGAACAGGTTTTTAACATGCAGCGCCACGTTGGATCCTGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCAGCCACGACCGTATAACTTCCGACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGaccacc                                         <  4:958747/150‑1 (MQ=255)
                                                                                                                ggTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGccc                                  <  8:1358107/151‑1 (MQ=255)
                                                                                                                    gttAACATGCCGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGTTGCCCAGc                               <  4:1403145/150‑1 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                            >  3:1117322/1‑150 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                            >  8:1034825/1‑150 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCa                           <  6:1777368/151‑1 (MQ=255)
                                                                                                                        aCATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAg                          <  2:1353221/151‑1 (MQ=255)
                                                                                                                        aCATGCAGCGCCACGTTGGAGCCACCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAg                          <  4:1209828/151‑1 (MQ=255)
                                                                                                                            gCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCt                      <  1:214706/151‑1 (MQ=255)
                                                                                                                                 gcCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGaa                   >  8:1781752/1‑149 (MQ=255)
                                                                                                                                  ccAGGTTGGAGTTAGCCTTGCGACTCTCGCGAGATGATAACGTATTAACGCATCCACGACCGTATAACTTCAAACTTTAAAGGTTTTAGCCGGATCAGCATAAGTGGATCCTCTATAACGACGACAGGTGCACAGCGTCATTCCTGAAAGc                <  4:773726/151‑1 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATGCGACTCTCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg               <  5:757973/151‑1 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATGCGACTCGCGCCAGCTCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGCCCACCGGTGCGCAGCGTCAGACCTGAAAGCg               >  2:28961/1‑151 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg               <  6:706928/151‑1 (MQ=255)
                                                                                                                                    aCGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCgg              <  4:412838/151‑1 (MQ=255)
                                                                                                                                     cGTTGGAGCAAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTAAAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAATTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGa             <  6:1582561/151‑1 (MQ=255)
                                                                                                                                        tgtagccagccATGCGACTCGCGCCAGATCATAACGTACTAACGATTCCACGCCCGTATAACTTCCCACTTTAAAGGAGTTAGCCTGATCACCATAACTGGAGCCAGTATAACGACCACCGGTTCCCAGCGTCAGACCTGAAAGCGGAcc           <  6:1713671/147‑1 (MQ=255)
                                                                                                                                           agccagccATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCa       <  6:471334/151‑1 (MQ=255)
                                                                                                                                            gccagccaTGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCa       >  1:1608567/1‑150 (MQ=255)
                                                                                                                                              cagccaTGCGACTCGCGCCAGATCATAACGTACTAACGCATACACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaa      >  8:1227895/1‑149 (MQ=255)
                                                                                                                                               agccaTGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaa      >  7:858460/1‑148 (MQ=255)
                                                                                                                                                gccaTGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCAaagaa  <  2:1331169/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCAAAGAA  >  DH1/3701636‑3701931

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: