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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,701,7650G.97.0% 1477.0 / 49.8 339coding (2129/2244 nt)fhuAferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base . (148/176);  total (154/185)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 342 total aligned reads displayed.

ATGGGCACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACC  >  DH1/3701617‑3701910
                                                                                                                                                    |                                                                                                                                                 
  gggCACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACg                                                                                                                                              <  6:316875/151‑1 (MQ=255)
    gCACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGtt                                                                                                                                            <  4:340839/151‑1 (MQ=255)
    gCACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGtt                                                                                                                                            <  3:782682/151‑1 (MQ=255)
      aCGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTgg                                                                                                                                          >  3:1055307/1‑151 (MQ=255)
      aCGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTgg                                                                                                                                          >  8:216940/1‑151 (MQ=255)
        ggAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTgg                                                                                                                                          >  4:1067461/1‑149 (MQ=255)
          aaaTCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGagcc                                                                                                                                      >  5:745003/1‑151 (MQ=255)
           aaTCCGTGCCCCATAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGagcc                                                                                                                                      >  2:1291100/1‑150 (MQ=255)
           aaTCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGagcca                                                                                                                                     <  5:622485/151‑1 (MQ=255)
             tCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGagccagc                                                                                                                                   >  7:1034406/1‑151 (MQ=255)
              ccGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCATCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGagccagc                                                                                                                                   <  3:277854/150‑1 (MQ=255)
               cGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAgccagcca                                                                                                                                 <  5:280802/151‑1 (MQ=255)
                   ccccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGcc                                                                                                                             <  7:1154121/151‑1 (MQ=255)
                      cAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGACAGCCATGCCGAc                                                                                                                          >  2:488768/1‑151 (MQ=255)
                       aaaaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACt                                                                                                                         >  5:208108/1‑151 (MQ=255)
                         aaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTcg                                                                                                                       >  7:230406/1‑151 (MQ=255)
                         aaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTcg                                                                                                                       <  4:1000377/151‑1 (MQ=255)
                          agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCCGCCATGCCGACTcgc                                                                                                                      >  6:946831/1‑151 (MQ=255)
                          agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTcgc                                                                                                                      >  5:68567/1‑151 (MQ=255)
                            agaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCgcgc                                                                                                                    >  1:929534/1‑151 (MQ=255)
                             gaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGCTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGcc                                                                                                                   >  8:1197684/1‑151 (MQ=255)
                              aaaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCAATAAGTGTTAAAGCAGCTGGCGACGTATTAACGATCGAACAGGTTGTTAACATGAAGC‑CAACGTTGGAGCCAGCCATGCCGACTCgcgca                                                                                                                   >  4:1107289/1‑149 (MQ=255)
                               aaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGCAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGcc                                                                                                                   >  1:1193695/1‑149 (MQ=255)
                               aaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGcc                                                                                                                   >  7:1296850/1‑149 (MQ=255)
                               aaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAg                                                                                                                 <  6:1082568/151‑1 (MQ=255)
                                aTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGa                                                                                                                <  7:1213398/151‑1 (MQ=255)
                                aTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGa                                                                                                                <  1:303059/151‑1 (MQ=255)
                                  tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGcc                                                                                                                   >  3:1183814/1‑146 (MQ=255)
                                     aaaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATa                                                                                                           <  7:1240401/151‑1 (MQ=255)
                                      aaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATaa                                                                                                          <  8:1077176/151‑1 (MQ=255)
                                           aaGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTAc                                                                                                     <  7:591314/151‑1 (MQ=255)
                                              gttgcggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTa                                                                                                   >  5:181442/1‑150 (MQ=255)
                                                  cggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGc                                                                                               >  5:1189425/1‑150 (MQ=255)
                                                     ttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                           >  3:1070411/1‑151 (MQ=255)
                                                      tgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                           <  4:540498/150‑1 (MQ=255)
                                                      tgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                           <  5:427587/150‑1 (MQ=255)
                                                       gcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCa                                                                                          <  5:1351646/150‑1 (MQ=255)
                                                       gcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCa                                                                                          <  6:783902/150‑1 (MQ=255)
                                                       gcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCAc                                                                                         <  7:950891/151‑1 (MQ=255)
                                                            gACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCg                                                                                    >  2:898782/1‑151 (MQ=255)
                                                             aCCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCTGC‑CCACGTTGTAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                   <  3:79126/151‑1 (MQ=255)
                                                              ccTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                   >  1:630707/1‑151 (MQ=255)
                                                               gtGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGGGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATCCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                 <  8:1250903/150‑1 (MQ=255)
                                                                tGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATAATAACGTACTAACGCATCCACGACCGtat                                                                                 >  1:447126/1‑150 (MQ=255)
                                                                  aCGTTCTGCGCCCCAGAAGCAGCCATAAGTTTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAAc                                                                              <  7:1114809/151‑1 (MQ=255)
                                                                    gTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACtt                                                                            >  1:1363899/1‑151 (MQ=255)
                                                                        tGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTccc                                                                         >  2:427661/1‑150 (MQ=255)
                                                                         gcgcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTccc                                                                         <  4:679856/149‑1 (MQ=255)
                                                                         gcgcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTcc                                                                          <  1:25770/148‑1 (MQ=255)
                                                                            ccccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACttt                                                                    >  3:384474/1‑151 (MQ=255)
                                                                             cccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACCGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTa                                                                   <  7:363892/151‑1 (MQ=255)
                                                                             cccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACttt                                                                    >  5:577639/1‑150 (MQ=255)
                                                                               caaaaGCACCCATAAGTGTAAAGCCCGGTGCCGACTTATTCACGATCTAACAGGTTGTTACCATGCAGC‑CGAGTTTGGAGCCAGCCTTGCCGACGCGCGCCATATAAAAACGTACAACCCCCCCCACGACCGTATAACTTCCCACTTTaaa                                                                 <  7:472545/148‑1 (MQ=255)
                                                                                 gAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAgg                                                               >  2:778497/1‑151 (MQ=255)
                                                                                   agcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAgg                                                               >  8:620043/1‑149 (MQ=255)
                                                                                   agcagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAg                                                             <  6:1224750/151‑1 (MQ=255)
                                                                                      agcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCATATCATAACGTACTAAAGCATCCACGACCGTATAACGTCCCACTTTAAAGGAGtt                                                           >  1:725762/1‑150 (MQ=255)
                                                                                      agcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                           >  5:890015/1‑150 (MQ=255)
                                                                                       gcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAg                                                         <  2:1052784/151‑1 (MQ=255)
                                                                                        ccATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGc                                                        >  2:643201/1‑151 (MQ=255)
                                                                                        ccATAAGTGTTAAAGCAGCTGGCGACGTATTCACGAACGAACAGGTTGTTAACATGCAGC‑CCACGTTCGAGCCAGCCATGCCGACTCGCGCCACATCATAACCTACTAACGCATCCACGACCGTATAACATCCCACTTTAAAGGACTTAGc                                                        >  6:165595/1‑151 (MQ=255)
                                                                                         cATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGTGTTAGcc                                                       >  4:180305/1‑151 (MQ=255)
                                                                                          aTAAGTGTTAAAGGAGCTGTCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCg                                                      <  2:713256/151‑1 (MQ=255)
                                                                                          aTAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCg                                                      <  3:185697/151‑1 (MQ=255)
                                                                                          aTAAGTCTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCg                                                      <  8:432741/151‑1 (MQ=255)
                                                                                               ggtTAAAGCAGCGGGCGACGTATGCAGGAGCGAACAGGTTGTTAACAGGCAGC‑CCCGGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCa                                                 <  3:1300573/150‑1 (MQ=255)
                                                                                                   aaaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                             <  3:69956/151‑1 (MQ=255)
                                                                                                    aaGCAGCTGGCGACGTATTCACGATCGAACCGGTTGTTAACATGCCGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATa                                            <  6:595702/151‑1 (MQ=255)
                                                                                                     agcagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATaa                                           <  3:1228876/151‑1 (MQ=255)
                                                                                                       cagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACt                                         <  4:846011/151‑1 (MQ=255)
                                                                                                       cagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACt                                         <  1:898484/151‑1 (MQ=255)
                                                                                                        cgcTGGCGACGTATTCACGATCCAACAGGTTGTTAACATGCA‑CCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTg                                        <  3:432620/150‑1 (MQ=255)
                                                                                                         gcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGACAGATCATAACGTACTAACGCATGAACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCTGATCACCATAACTgg                                       <  1:1156369/151‑1 (MQ=255)
                                                                                                          cTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                       >  3:501994/1‑150 (MQ=255)
                                                                                                            ggCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCAGGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAAGGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGa                                      <  8:855938/149‑1 (MQ=255)
                                                                                                            ggCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAg                                     <  1:502864/150‑1 (MQ=255)
                                                                                                             gCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                   <  1:344950/151‑1 (MQ=255)
                                                                                                              ccaCGTATTCACGATCGAACAGGTTGATAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                  <  4:1136085/149‑1 (MQ=255)
                                                                                                                   tATTCACGATCGAACAGGTTGTTAACATGCA‑CCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCACATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCCGATCACCATAACAGGAGCCAGtat                              >  6:705955/1‑150 (MQ=255)
                                                                                                                    aTTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtata                             <  5:1233247/150‑1 (MQ=255)
                                                                                                                    aTTCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtata                             <  7:755393/150‑1 (MQ=255)
                                                                                                                    aTTCACGAGCGAAAAGGTTGTTAACATGCAGC‑CCCCGTTGTATCCAGCAATGCATACTCGCTCCAGATCATAACTTACTAACGCATCAACGACAGTATAACTTCCCACTTTAAATGAGTTATCCGTAGCACCATAACTGGAGCCAGTATaa                            <  4:123688/151‑1 (MQ=255)
                                                                                                                      tCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAAc                           <  5:1148206/150‑1 (MQ=255)
                                                                                                                      tCACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACg                          >  8:54448/1‑151 (MQ=255)
                                                                                                                       cACGATCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACg                          >  8:734801/1‑150 (MQ=255)
                                                                                                                           aTCGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGaccac                     >  1:168477/1‑151 (MQ=255)
                                                                                                                             cGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCCCCATAACTGGAGCCAGTATAACGACCACCg                   >  8:330315/1‑151 (MQ=255)
                                                                                                                             cGAACAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCg                   >  2:1348516/1‑151 (MQ=255)
                                                                                                                               aaCAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGt                 <  1:1005494/151‑1 (MQ=255)
                                                                                                                                aCAGTTTGTTAACATGCAGC‑CCACGTTGGAGCCATCCATGCGGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTg                <  3:747120/151‑1 (MQ=255)
                                                                                                                                aCAGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGCCCACCGGTg                >  4:788218/1‑151 (MQ=255)
                                                                                                                                  aGGTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGCTGcc              >  8:467150/1‑151 (MQ=255)
                                                                                                                                   ggTTGTTAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAAAACTGGAGCCAGTATAACGACCACCGGTGccc             <  3:1125719/151‑1 (MQ=255)
                                                                                                                                    gttgttAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCa            <  8:14622/151‑1 (MQ=255)
                                                                                                                                       gttACCATGCAGC‑CCACGTTGGCCCCAGCCATGCCGACTGGCGCAAGATCTAAACGTACTAACGAATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCg         <  2:870472/151‑1 (MQ=255)
                                                                                                                                         tAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTACCGCCTCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGCCCACCGGTGCCCAGCGTc       >  5:1013338/1‑151 (MQ=255)
                                                                                                                                         tAACATGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc       >  7:873586/1‑151 (MQ=255)
                                                                                                                                           aCATGCAGC‑CCCCGGTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTCCTAACGCATCCACGCCCGGATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAg     <  6:1374157/151‑1 (MQ=255)
                                                                                                                                            cATGCAGC‑CCACGTTGGAGCCCGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCCCCGGTGCCCAGCGTCAgc    >  6:24473/1‑150 (MQ=255)
                                                                                                                                             aTGCAGC‑CCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGAc   <  7:282689/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATGGGCACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACC  >  DH1/3701617‑3701910

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