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.99.5% 2017.8 / 2.6 434coding (2129/2244 nt)fhuAferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base . (221/210);  minor base C (0/2);  total (221/213)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 434 total aligned reads displayed.

GGCACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACC  >  DH1/3701620‑3701910
                                                                                                                                                 |                                                                                                                                                 
  cACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTg                                                                                                                                           >  3:1238453/1‑151 (MQ=255)
       aaaTCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGagcc                                                                                                                                      <  1:802853/151‑1 (MQ=255)
         atacGTGCCCCAAAAGAGAAATTAGAAAAGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCACAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGCTTGTTAACATGCAGCCCACGTTGGagcc                                                                                                                                      >  1:252105/4‑149 (MQ=255)
         aTCCGTGCCCCAAAAGAGAAATTAGAAAGGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGGAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCAGGATCGAACAGGTTGTTAACATGCAGCCCAGGTTGGagcg                                                                                                                                      >  4:248269/1‑148 (MQ=255)
           ccGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGagccagc                                                                                                                                   <  6:1783676/150‑1 (MQ=255)
             gTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCAt                                                                                                                                <  6:1465597/151‑1 (MQ=255)
              tGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATg                                                                                                                               >  2:929821/1‑151 (MQ=255)
               gCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGc                                                                                                                              >  3:1326424/1‑151 (MQ=255)
                ccccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGcc                                                                                                                             >  8:1244583/1‑151 (MQ=255)
                 cccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCg                                                                                                                            <  8:295204/151‑1 (MQ=255)
                   cAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGACGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGAc                                                                                                                          <  2:818699/151‑1 (MQ=255)
                    aaaaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACt                                                                                                                         >  5:1084016/1‑151 (MQ=255)
                    aaaaGAGAAATTAGAAACGGAAGGTTGCGGGTGCAACGACCTGACGGTCTGCGCCCCCGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACt                                                                                                                         <  8:614314/151‑1 (MQ=255)
                     aaaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTc                                                                                                                        >  1:1384935/1‑151 (MQ=255)
                      aaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTTAAGCAGCTGGCGCCGTATTCACGATCGCACAGGGTGTTACCATGCAGCCCACGTTGGCGCCCGCCATGCCGCCTcg                                                                                                                       >  2:651386/1‑151 (MQ=255)
                      aaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGCCTcg                                                                                                                       >  7:133305/1‑151 (MQ=255)
                      aaGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTcg                                                                                                                       >  4:1314948/1‑151 (MQ=255)
                       agagAAATTAGAAAGGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTcgc                                                                                                                      <  8:495719/151‑1 (MQ=255)
                       agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTcgc                                                                                                                      >  2:1171790/1‑151 (MQ=255)
                       agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTcgc                                                                                                                      >  5:1502445/1‑151 (MQ=255)
                        gagaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTcgcg                                                                                                                     <  5:1371599/151‑1 (MQ=255)
                         agaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCgcgc                                                                                                                    >  8:1638327/1‑151 (MQ=255)
                          gaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGcc                                                                                                                   <  5:1210521/151‑1 (MQ=255)
                            aaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGcc                                                                                                                   >  7:1668765/1‑149 (MQ=255)
                             aTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGcc                                                                                                                   >  6:1087840/1‑148 (MQ=255)
                             aTTAGAAACGGAAGGTTGCGGGTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGa                                                                                                                <  2:1266381/151‑1 (MQ=255)
                               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAAGAGGTTGTTAACATGCAGCCCAAGATGGAGCCAGCCATGCCGACTCGAGCCACATc                                                                                                              >  8:1078951/1‑151 (MQ=255)
                               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGAt                                                                                                               <  5:1004778/150‑1 (MQ=255)
                                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGCTc                                                                                                              >  7:29375/1‑150 (MQ=255)
                                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATc                                                                                                              >  5:155659/1‑150 (MQ=255)
                                 gAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCAt                                                                                                            <  7:1841358/151‑1 (MQ=255)
                                 caaaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCAt                                                                                                            <  2:1792202/150‑1 (MQ=255)
                                  aaaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCAt                                                                                                            >  6:967492/1‑150 (MQ=255)
                                    aCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAAc                                                                                                         <  7:1439227/151‑1 (MQ=255)
                                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACg                                                                                                        >  5:591730/1‑151 (MQ=255)
                                      ggAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGCCTCGCGCCAGATCATAACgg                                                                                                       >  2:528135/1‑150 (MQ=255)
                                      ggAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGt                                                                                                       >  2:1059510/1‑151 (MQ=255)
                                      ggAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGt                                                                                                       <  1:958648/151‑1 (MQ=255)
                                       gacggttgcggttgcAACGACCTGACGTTCTGCGCCCCCGAAGAAGCCATAAGTGTTACAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTa                                                                                                      <  4:636602/148‑1 (MQ=255)
                                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTa                                                                                                      <  4:1665909/151‑1 (MQ=255)
                                          ggttgcggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTa                                                                                                   <  3:332812/151‑1 (MQ=255)
                                           gttgcggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTaa                                                                                                  <  4:1468772/151‑1 (MQ=255)
                                             tgcggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAAc                                                                                                 <  1:1063346/150‑1 (MQ=255)
                                              gcggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACg                                                                                                <  2:252994/150‑1 (MQ=255)
                                                    gcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCAc                                                                                         <  7:471420/151‑1 (MQ=255)
                                                        cGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGAcc                                                                                     >  4:1168604/1‑151 (MQ=255)
                                                         gACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCg                                                                                    <  5:430286/151‑1 (MQ=255)
                                                         gACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCg                                                                                    >  3:18035/1‑151 (MQ=255)
                                                            cTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                 >  4:538834/1‑151 (MQ=255)
                                                              gACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                 <  4:1777749/149‑1 (MQ=255)
                                                              gACGTTCTGCGCCCCACACGCAGCAATAAGTGTTAACGCAGGTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCCGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                 <  8:780839/149‑1 (MQ=255)
                                                                 gTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACtt                                                                            >  6:214075/1‑151 (MQ=255)
                                                                 gTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAAAGCATCCACGACCGTATAACtt                                                                            >  2:478171/1‑151 (MQ=255)
                                                                    atGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTccc                                                                         <  1:1689034/150‑1 (MQ=255)
                                                                       cgcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACt                                                                      >  6:301876/1‑151 (MQ=255)
                                                                        gcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGTTTGTTAACATGCAGCCCACGTTGGAGCCAGCAATGCCGACTCGCCACAGATCATAACGTACTAACGCATCCACGACCGAATAAATACCCACtt                                                                     >  8:1761054/1‑151 (MQ=255)
                                                                        gcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGCCTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACtt                                                                     >  3:1544417/1‑151 (MQ=255)
                                                                           ccAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTaa                                                                  <  7:1775240/151‑1 (MQ=255)
                                                                            cAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTaaa                                                                 <  3:401968/151‑1 (MQ=255)
                                                                            cAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTaa                                                                  >  6:517465/1‑150 (MQ=255)
                                                                             aGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAg                                                                <  4:461345/151‑1 (MQ=255)
                                                                                 ccagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                            <  6:1342224/150‑1 (MQ=255)
                                                                                  cagcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                           <  7:1245486/151‑1 (MQ=255)
                                                                                   agcCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                           >  3:1493007/1‑150 (MQ=255)
                                                                                     ccATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGc                                                        >  1:694630/1‑151 (MQ=255)
                                                                                      cATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACGGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGcc                                                       >  5:785344/1‑151 (MQ=255)
                                                                                      cATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGcc                                                       >  1:54061/1‑151 (MQ=255)
                                                                                      cATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGcc                                                       >  4:1083843/1‑151 (MQ=255)
                                                                                          attgtTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGCCCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGAt                                                   >  7:61011/3‑151 (MQ=255)
                                                                                          aGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGAt                                                   <  5:244417/151‑1 (MQ=255)
                                                                                            tgtTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGAATc                                                  >  3:1286042/1‑150 (MQ=255)
                                                                                             gtTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAc                                                >  3:592498/1‑151 (MQ=255)
                                                                                              ttAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAcc                                               >  1:1523791/1‑151 (MQ=255)
                                                                                                  agccgcTGGCGACGTATTCAGGATGCAACAGGTTGTTAACATGCAGCCCACGTTGGCGCCCGACTTGCCGAATCGCGCCAGATCATAACGTACTAACGCACCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATaa                                           <  1:845592/151‑1 (MQ=255)
                                                                                                   gcagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAAc                                          <  5:898149/151‑1 (MQ=255)
                                                                                                    cagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCAAAGTTTTAGCAAGCCATGCCGACTCGATCCAGATCATAAAGTACTAAATCATCCAAGAACGTATAACTTCAAACTTTAAATGAGATATCCGGATCACCATAACt                                         >  3:1042708/1‑151 (MQ=255)
                                                                                                     agcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTg                                        >  3:852770/1‑151 (MQ=255)
                                                                                                     agcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTg                                        <  4:1544841/151‑1 (MQ=255)
                                                                                                         ggCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCAAGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGc                                    >  7:1394524/1‑151 (MQ=255)
                                                                                                          gCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTGTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                   >  3:758102/1‑151 (MQ=255)
                                                                                                             aCGTATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGt                                <  8:1347160/151‑1 (MQ=255)
                                                                                                             aCGTATTCAAGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCTCGCCAGATCATAACGTAATAACGCATCCACGACAGTATAAATTCCCTATTTAAAGGAGTTAGCCGGATAACCATAACTGGAGCCact                                >  1:59750/1‑149 (MQ=255)
                                                                                                                tATTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtat                              >  4:794560/1‑150 (MQ=255)
                                                                                                                 aTTCACGATCTAACATGTTGTTAACATGCAGCCCACGTTGGACCCAGCCATGCCGACCCGCGCCAGATCATCACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATaa                            <  3:1701337/151‑1 (MQ=255)
                                                                                                                 aTTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtata                             >  7:93017/1‑150 (MQ=255)
                                                                                                                 aTTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATaa                            <  6:864745/151‑1 (MQ=255)
                                                                                                                 aTTCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATaa                            <  1:410324/151‑1 (MQ=255)
                                                                                                                   tCACGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACg                          >  6:1407202/1‑151 (MQ=255)
                                                                                                                      cGATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGacc                       >  1:1228833/1‑151 (MQ=255)
                                                                                                                       gATCGAACAGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGacc                       >  2:1422407/1‑150 (MQ=255)
                                                                                                                           caaCAGGTTGTTAACATGCAGCCCTCGTTGGAGCCAGCCATGCCGACGCGCGCCACTTCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCgg                  <  7:1075523/150‑1 (MQ=255)
                                                                                                                               aGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGcc              >  5:5158/1‑151 (MQ=255)
                                                                                                                               aGGTTGTTAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGcc              >  7:1049368/1‑151 (MQ=255)
                                                                                                                                 gttgttAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGccc             >  5:1760039/1‑150 (MQ=255)
                                                                                                                                 gttgttAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCa            <  1:1169429/151‑1 (MQ=255)
                                                                                                                                   tgttAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAg           <  1:736258/150‑1 (MQ=255)
                                                                                                                                   tgttAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAg           <  1:1618356/150‑1 (MQ=255)
                                                                                                                                    gttAACATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGc          <  8:1475457/150‑1 (MQ=255)
                                                                                                                                       aaCATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCa      <  2:1609277/151‑1 (MQ=255)
                                                                                                                                         cATGCAGCCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGa    <  2:1019091/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGCACGGAAATCCGTGCCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACC  >  DH1/3701620‑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: