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.96.5% 1295.1 / 40.6 288coding (2108/2244 nt)fhuAferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base . (138/140);  total (145/143)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

CCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCAAAGAAGGTG  >  DH1/3701636‑3701935
                                                                                                                                                      |                                                                                                                                                     
ccccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGc                                                                                                                                                       <  5:697126/151‑1 (MQ=255)
 cccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGcg                                                                                                                                                      >  7:1045461/1‑150 (MQ=255)
 cccAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGcg                                                                                                                                                      >  4:1120703/1‑150 (MQ=255)
       agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTcgc                                                                                                                                               >  1:634545/1‑151 (MQ=255)
        gagaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTcgcg                                                                                                                                              >  5:228235/1‑151 (MQ=255)
         agaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTGACGATCGAACAGGTTGTTAACCTGCAGCGCCACGTTGGAGCCAGCCATGCGCCTCgcgc                                                                                                                                             >  2:725565/1‑151 (MQ=255)
           aaaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGcc                                                                                                                                            >  7:782344/1‑150 (MQ=255)
           aaaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGcc                                                                                                                                            >  3:899636/1‑150 (MQ=255)
           aaaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCa                                                                                                                                           <  1:461826/151‑1 (MQ=255)
               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGAt                                                                                                                                        <  2:848033/150‑1 (MQ=255)
               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGAt                                                                                                                                        <  7:88470/150‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCa                                                                                                                                      <  4:337534/151‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATCCCGACTCGCGcc                                                                                                                                            >  1:924619/1‑146 (MQ=255)
                  aaaCGGAAGGTTGCGGCTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATa                                                                                                                                    <  7:1252209/151‑1 (MQ=255)
                   aaCGGAAGGTTGCGGTTGCGCCGACCTGACCTTCTGCGCCCCAGAAGCCGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCCCGATCGAACACGGTGTTAACATGCAGCGCCACGTTGGAGCCAGCCCTGCGACTCGCGCCAGAGCATa                                                                                                                                    >  4:717126/1‑150 (MQ=255)
                   aaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATa                                                                                                                                    >  2:686042/1‑150 (MQ=255)
                    aCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAAc                                                                                                                                  >  1:747612/1‑151 (MQ=255)
                    aCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAAc                                                                                                                                  <  6:669988/151‑1 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAAGATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACg                                                                                                                                 >  5:794532/1‑151 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACg                                                                                                                                 <  3:527957/151‑1 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGt                                                                                                                                >  4:410350/1‑150 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTa                                                                                                                               <  3:1238119/151‑1 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTa                                                                                                                               <  3:76064/151‑1 (MQ=255)
                              gcggttgcAACGACCTGACGTTCTGCGGCCCAGAATCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGc                                                                                                                        <  2:171646/151‑1 (MQ=255)
                               cggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGc                                                                                                                        >  6:896930/1‑151 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCAt                                                                                                                      >  8:1219349/1‑151 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCAt                                                                                                                      <  5:462115/151‑1 (MQ=255)
                                  ttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGAGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                                                    >  8:1082707/1‑151 (MQ=255)
                                   tgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCCGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                                                    >  7:917521/1‑150 (MQ=255)
                                    gcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCa                                                                                                                   <  2:966179/150‑1 (MQ=255)
                                      aaCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCGTCCACgc                                                                                                                >  3:199504/1‑150 (MQ=255)
                                      aaCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGa                                                                                                                <  4:1209546/151‑1 (MQ=255)
                                         gACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGCTCATAACGTACTAACGCATCCACGACCg                                                                                                             >  5:636943/1‑151 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                            <  4:882676/151‑1 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                            >  1:1029314/1‑151 (MQ=255)
                                              gACGTTCTGCGCCCCATAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtata                                                                                                         <  8:535451/150‑1 (MQ=255)
                                              gACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                                          <  6:1071070/149‑1 (MQ=255)
                                                      gcgcCCCTTAAGCAGCCATAAGTGTTAAAGCAGCTGGCGAGGTATTCACGATCGAACAGGTTGTTAACATGCAGGGCCACGTTGGAGCCAGCCATGGGACTCGCGCCAGATCATAAGTTACTAACGCATCCACGACCGTATAACTTccc                                                                                                  <  1:1100289/149‑1 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCATAAGTGTTAACGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCa                                                                                                 <  7:1195108/150‑1 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCa                                                                                                 <  5:371744/150‑1 (MQ=255)
                                                       cgcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACt                                                                                               >  7:469127/1‑151 (MQ=255)
                                                       cgcCCCAGAAGCAGCCATAAGTGTAAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTAAACATGCCGCGCCAGGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCAGCCACGACCGTATAACTTCCCACt                                                                                               <  7:71619/151‑1 (MQ=255)
                                                        gcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCTACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACtt                                                                                              <  8:582328/151‑1 (MQ=255)
                                                        gcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCAcgt                                                                                              >  6:810740/1‑149 (MQ=255)
                                                                 gcagcCATAAGTGTTAAAGCAGCTGGCTACTTATTCACTATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCAATGCGAGTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                     <  1:336239/151‑1 (MQ=255)
                                                                  cagcCATAAGTGTTAAAGCAGCTGGGGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGCGCCAACCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                                                    <  5:847522/151‑1 (MQ=255)
                                                                   agcCATAAGTGTTAAAGCACCTGGCGACGAATTCACGATCGAACAGGTTGTAAAGATGAAGCGCCACGTTGGAGCCAGCGATCCGACTGGGGCCAGATCTTAAGGAACTAACCCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTa                                                                                   <  3:163891/151‑1 (MQ=255)
                                                                     ccATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGc                                                                                 >  6:463951/1‑151 (MQ=255)
                                                                       aTAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCg                                                                               >  5:1215128/1‑151 (MQ=255)
                                                                         aaGTGTTAAACCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCgg                                                                              <  8:1136805/150‑1 (MQ=255)
                                                                           gtgtTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATc                                                                           <  8:964655/151‑1 (MQ=255)
                                                                             gtTAAAGCAGCTGGCGACGTATTCACGATCGTACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAc                                                                         <  8:1167129/151‑1 (MQ=255)
                                                                              ttAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAcc                                                                        >  7:585269/1‑151 (MQ=255)
                                                                                aaaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGACGGATCAACAt                                                                      >  4:202069/1‑151 (MQ=255)
                                                                                 atggagcTGGCTTCGTATTTACGTTCGAACATGTTGTTAACATGCAGCGCCACGTTGTAGCGAGCCATGCGACTCGCGCCATATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                      <  2:1120856/146‑1 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                      >  1:211725/1‑150 (MQ=255)
                                                                                    cagcTGGCGACGTATTCAGGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACGCGCGCCAGATCATAACGTACTAACGCATGCCCGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCCCCATaagt                                                                  >  6:1175285/1‑149 (MQ=255)
                                                                                    cagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACt                                                                  >  8:261200/1‑151 (MQ=255)
                                                                                     cgcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTg                                                                 <  8:801608/150‑1 (MQ=255)
                                                                                      gcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                                >  5:257084/1‑151 (MQ=255)
                                                                                      gcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                                <  2:134074/151‑1 (MQ=255)
                                                                                       cTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                                >  4:1265156/1‑150 (MQ=255)
                                                                                         ggCGACGTATTCCCGAGCGAACAGGGTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAg                                                              <  4:587732/150‑1 (MQ=255)
                                                                                         ggCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGa                                                               <  6:124259/149‑1 (MQ=255)
                                                                                         ggCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAg                                                              <  7:998192/150‑1 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACCGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                           <  2:1149930/151‑1 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                           <  1:1212802/151‑1 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATACCGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                                            >  6:778026/1‑149 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                                            >  1:1049108/1‑149 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGGTAGCCGGATCACCATAACTGGAGcc                                                            >  8:1242886/1‑149 (MQ=255)
                                                                                              cGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGta                                                        <  4:995691/151‑1 (MQ=255)
                                                                                              cGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGt                                                         >  8:388125/1‑150 (MQ=255)
                                                                                               gTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtat                                                       >  7:668012/1‑151 (MQ=255)
                                                                                                 aTTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATaa                                                     <  1:794086/151‑1 (MQ=255)
                                                                                                    cACGATGGAACAGGTTGTTAACATGCAGCGCCAGGTTGCAGCCAGCCATGCGACTCGCGCCAGATAATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACg                                                   <  3:607570/150‑1 (MQ=255)
                                                                                                    cACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCAAGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACg                                                   >  6:584448/1‑150 (MQ=255)
                                                                                                     aCGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGCCCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGac                                                 >  5:53729/1‑151 (MQ=255)
                                                                                                      cGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGacc                                                <  3:657395/151‑1 (MQ=255)
                                                                                                        aTCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGcccac                                              >  5:95772/1‑151 (MQ=255)
                                                                                                           gAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCgg                                           <  4:510687/151‑1 (MQ=255)
                                                                                                            aaCAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCAAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGt                                          <  8:196688/151‑1 (MQ=255)
                                                                                                              cAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGc                                        <  8:227057/151‑1 (MQ=255)
                                                                                                                 gttgttAACATGCAGCGCCACGTTGTAGCCAGCCATGCGACTCGCGCCACATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCa                                     <  8:471266/151‑1 (MQ=255)
                                                                                                                 gttgttAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGccc                                      >  5:111815/1‑150 (MQ=255)
                                                                                                                      tAACATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                                >  5:1055528/1‑151 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCa                               <  6:367834/151‑1 (MQ=255)
                                                                                                                         cATGCAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                                >  4:94680/1‑148 (MQ=255)
                                                                                                                             cAGCGCCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTg                         <  6:521909/151‑1 (MQ=255)
                                                                                                                                cgcCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGaaa                      <  7:176947/151‑1 (MQ=255)
                                                                                                                                 gcCACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     <  8:1288817/151‑1 (MQ=255)
                                                                                                                                 gcCACGTGGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     <  7:238840/151‑1 (MQ=255)
                                                                                                                                 cccACGGTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     <  8:856312/150‑1 (MQ=255)
                                                                                                                                  ccACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGc                    >  1:600645/1‑151 (MQ=255)
                                                                                                                                  ccACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGc                    >  7:319527/1‑151 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     >  2:419829/1‑149 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg                   >  3:969731/1‑151 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg                   <  6:247694/151‑1 (MQ=255)
                                                                                                                                         ggAGCCAGCCATGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCg              <  2:142153/150‑1 (MQ=255)
                                                                                                                                             ccagccaTGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaaa         <  2:226723/151‑1 (MQ=255)
                                                                                                                                             ccagccaTGCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaaa         <  7:1037167/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCCCAAAAGAGAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCAAAGAAGGTG  >  DH1/3701636‑3701935

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