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)
     aaaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTc                                                                                                                                             >  5:1640048/1‑151 (MQ=255)
       agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTcgc                                                                                                                                           >  5:1212016/1‑151 (MQ=255)
        gagaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCCCGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTcgcg                                                                                                                                          >  5:982086/1‑151 (MQ=255)
        gagaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTcgcg                                                                                                                                          <  1:1541707/151‑1 (MQ=255)
         agaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCgcgc                                                                                                                                         <  6:564685/151‑1 (MQ=255)
             aTTAGAAACGGAAGGTTGCGGTTTCAACGTACTGTAGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGa                                                                                                                                     <  8:550643/151‑1 (MQ=255)
             aTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGcc                                                                                                                                        >  7:1440369/1‑148 (MQ=255)
               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGAt                                                                                                                                    <  2:1609671/150‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCa                                                                                                                                  <  7:276680/151‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCa                                                                                                                                  <  1:1184395/151‑1 (MQ=255)
                   aaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATaa                                                                                                                               <  7:1534276/151‑1 (MQ=255)
                    aCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAAc                                                                                                                              <  1:548827/151‑1 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACg                                                                                                                             >  5:1542954/1‑151 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACg                                                                                                                             >  4:1623379/1‑151 (MQ=255)
                      ggAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGt                                                                                                                            >  6:381327/1‑151 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGt                                                                                                                            >  7:1183451/1‑150 (MQ=255)
                          ggttgcggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACt                                                                                                                         >  3:410943/1‑150 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCAt                                                                                                                  >  5:1343781/1‑151 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGAGTTAAACCAGCTCGCGACGTATTCACGATCGAACAGGTTGATAACATGCAGCGCCACGTTCCAGCCAGCCATGCCGACACGCGCCAGATCCTAACGTACAAAccc                                                                                                                    >  2:80327/1‑148 (MQ=255)
                                   tgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                                                >  3:1803121/1‑150 (MQ=255)
                                       aCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGAc                                                                                                           <  1:1157502/151‑1 (MQ=255)
                                       aCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGAc                                                                                                           >  1:1392567/1‑151 (MQ=255)
                                       aCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGAAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGAc                                                                                                           >  8:247200/1‑151 (MQ=255)
                                         gACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCg                                                                                                         >  5:871097/1‑151 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                        <  3:355115/151‑1 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCACAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACCATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                        <  8:1144035/151‑1 (MQ=255)
                                           ccTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                        >  2:1518335/1‑150 (MQ=255)
                                              gACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACTATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                                      <  7:1265065/149‑1 (MQ=255)
                                               aCGTTCTGCGCCCCAGAAGCAGCCATAAGTTTTAAAGCAGCTGGCGACTTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATaa                                                                                                    <  7:693841/151‑1 (MQ=255)
                                                  ttCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTc                                                                                                >  1:1104017/1‑151 (MQ=255)
                                                     tGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTccc                                                                                              >  5:311839/1‑150 (MQ=255)
                                                     tGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGGATAACTTccc                                                                                              >  8:1454540/1‑150 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCCTAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCAc                                                                                            <  6:1418655/151‑1 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCATAAGTGTTAAAGGAGCTGGCGACGTATTCACGATCGAACTGGTTGTTAACATGCAGCGCCAGGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCAc                                                                                            <  4:1753834/151‑1 (MQ=255)
                                                         ccccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACttt                                                                                         >  8:715932/1‑151 (MQ=255)
                                                          cccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTa                                                                                        <  3:1659913/151‑1 (MQ=255)
                                                                agcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAgg                                                                                    >  6:914769/1‑149 (MQ=255)
                                                                 gcagcCATAAGTTTTAAAGAAGCTGGCGACGTATTCCCGATCGAACAGGTTGTTAACATGCCGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGTTCATAACGTACTAACGAATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  7:1798158/151‑1 (MQ=255)
                                                                 gcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  4:1491369/151‑1 (MQ=255)
                                                                 gcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  2:596651/151‑1 (MQ=255)
                                                                  cagcCATAAGTGTTAAAGCAGCTGGCGACTTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                 <  5:44838/151‑1 (MQ=255)
                                                                  cagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCAT‑CCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                                                >  8:1637746/1‑151 (MQ=255)
                                                                   agcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                                                >  5:141345/1‑150 (MQ=255)
                                                                     ccATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGc                                                                             <  1:61424/151‑1 (MQ=255)
                                                                     caaTAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCAGAACGTACTAACGCATCCACGCCAGTATAACTTCCCACTTTAAAGGAGTTatc                                                                             >  7:1610083/3‑149 (MQ=255)
                                                                         aaGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCgg                                                                          <  1:307833/151‑1 (MQ=255)
                                                                         aaGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCgg                                                                          <  6:201235/150‑1 (MQ=255)
                                                                               tAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAcc                                                                    <  8:812356/150‑1 (MQ=255)
                                                                                aaaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCa                                                                   <  2:171311/151‑1 (MQ=255)
                                                                                aaaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                  <  5:1743588/151‑1 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                  >  3:99315/1‑150 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATa                                                                 <  4:1803689/151‑1 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTAAAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                  >  6:1101161/1‑150 (MQ=255)
                                                                                  agcagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATaa                                                                <  4:348062/151‑1 (MQ=255)
                                                                                   gcagcTGGCGACGTATTCCCGCTCGAACAGGTTGTTACCATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCCTCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGACGTCTCACCATCAc                                                               >  6:1283821/1‑151 (MQ=255)
                                                                                    cagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACt                                                              >  4:1049704/1‑151 (MQ=255)
                                                                                      gcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            <  3:508650/151‑1 (MQ=255)
                                                                                       cTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  2:251636/1‑150 (MQ=255)
                                                                                       cTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  7:1119380/1‑150 (MQ=255)
                                                                                        tGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  5:1847820/1‑149 (MQ=255)
                                                                                        tGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                            >  6:3341/1‑149 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                       <  7:1371141/151‑1 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                       <  3:815385/151‑1 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                                        >  2:890852/1‑149 (MQ=255)
                                                                                             aCGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGt                                                     <  5:712379/151‑1 (MQ=255)
                                                                                                 aTTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtat                                                   >  7:322810/1‑150 (MQ=255)
                                                                                                  ttCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAAc                                                >  6:784588/1‑151 (MQ=255)
                                                                                                   tCACGATCTAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAAc                                                <  1:1703242/150‑1 (MQ=255)
                                                                                                    cACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACgc                                              >  2:454178/1‑150 (MQ=255)
                                                                                                     aCGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGac                                             >  4:560706/1‑151 (MQ=255)
                                                                                                      cGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGacc                                            >  4:239553/1‑151 (MQ=255)
                                                                                                        aTCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGaccac                                          >  1:809774/1‑151 (MQ=255)
                                                                                                         tCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGaccacc                                         >  5:633371/1‑151 (MQ=255)
                                                                                                          cGACCAGGTTGTTAACATGCAGCGCCACGTTGCAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAATTGTAGCCAGTATAACGACCACCg                                        <  4:1269171/151‑1 (MQ=255)
                                                                                                          cGAACAGGTTTTTAACATGCAGCGCCACGTTGGATCCTGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCAGCCACGACCGTATAACTTCCGACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGaccacc                                         <  4:958747/150‑1 (MQ=255)
                                                                                                                ggTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGccc                                  <  8:1358107/151‑1 (MQ=255)
                                                                                                                    gttAACATGCCGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGTTGCCCAGc                               <  4:1403145/150‑1 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                            >  3:1117322/1‑150 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                            >  8:1034825/1‑150 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCa                           <  6:1777368/151‑1 (MQ=255)
                                                                                                                        aCATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAg                          <  2:1353221/151‑1 (MQ=255)
                                                                                                                        aCATGCAGCGCCACGTTGGAGCCACCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAg                          <  4:1209828/151‑1 (MQ=255)
                                                                                                                            gCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCt                      <  1:214706/151‑1 (MQ=255)
                                                                                                                                 gcCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGaa                   >  8:1781752/1‑149 (MQ=255)
                                                                                                                                  ccAGGTTGGAGTTAGCCTTG‑CGACTCTCGCGAGATGATAACGTATTAACGCATCCACGACCGTATAACTTCAAACTTTAAAGGTTTTAGCCGGATCAGCATAAGTGGATCCTCTATAACGACGACAGGTGCACAGCGTCATTCCTGAAAGc                <  4:773726/151‑1 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATG‑CGACTCTCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg               <  5:757973/151‑1 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGCTCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGCCCACCGGTGCGCAGCGTCAGACCTGAAAGCg               >  2:28961/1‑151 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg               <  6:706928/151‑1 (MQ=255)
                                                                                                                                    aCGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCgg              <  4:412838/151‑1 (MQ=255)
                                                                                                                                     cGTTGGAGCAAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTAAAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAATTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGa             <  6:1582561/151‑1 (MQ=255)
                                                                                                                                        tgtagccagccATG‑CGACTCGCGCCAGATCATAACGTACTAACGATTCCACGCCCGTATAACTTCCCACTTTAAAGGAGTTAGCCTGATCACCATAACTGGAGCCAGTATAACGACCACCGGTTCCCAGCGTCAGACCTGAAAGCGGAcc           <  6:1713671/147‑1 (MQ=255)
                                                                                                                                           agccagccATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCa       <  6:471334/151‑1 (MQ=255)
                                                                                                                                            gccagccaTG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCa       >  1:1608567/1‑150 (MQ=255)
                                                                                                                                              cagccaTG‑CGACTCGCGCCAGATCATAACGTACTAACGCATACACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaa      >  8:1227895/1‑149 (MQ=255)
                                                                                                                                               agccaTG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaa      >  7:858460/1‑148 (MQ=255)
                                                                                                                                                gccaTG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCAaagaa  <  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: ‑