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)
       agagAAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTcgc                                                                                                                                               >  1:634545/1‑151 (MQ=255)
        gagaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTcgcg                                                                                                                                              >  5:228235/1‑151 (MQ=255)
         agaAATTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTGACGATCGAACAGGTTGTTAACCTGCAGCGCCACGTTGGAGCCAGCCATG‑CGCCTCgcgc                                                                                                                                             >  2:725565/1‑151 (MQ=255)
           aaaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGcc                                                                                                                                            >  7:782344/1‑150 (MQ=255)
           aaaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGcc                                                                                                                                            >  3:899636/1‑150 (MQ=255)
           aaaTTAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCa                                                                                                                                           <  1:461826/151‑1 (MQ=255)
               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGAt                                                                                                                                        <  2:848033/150‑1 (MQ=255)
               tAGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGAt                                                                                                                                        <  7:88470/150‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCa                                                                                                                                      <  4:337534/151‑1 (MQ=255)
                aGAAACGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATCCCGACTCGCGcc                                                                                                                                            >  1:924619/1‑146 (MQ=255)
                  aaaCGGAAGGTTGCGGCTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATa                                                                                                                                    <  7:1252209/151‑1 (MQ=255)
                   aaCGGAAGGTTGCGGTTGCGCCGACCTGACCTTCTGCGCCCCAGAAGCCGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCCCGATCGAACACGGTGTTAACATGCAGCGCCACGTTGGAGCCAGCCCTG‑CGACTCGCGCCAGAGCATa                                                                                                                                    >  4:717126/1‑150 (MQ=255)
                   aaCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATa                                                                                                                                    >  2:686042/1‑150 (MQ=255)
                    aCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAAc                                                                                                                                  >  1:747612/1‑151 (MQ=255)
                    aCGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAAc                                                                                                                                  <  6:669988/151‑1 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAAGATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACg                                                                                                                                 >  5:794532/1‑151 (MQ=255)
                     cGGAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACg                                                                                                                                 <  3:527957/151‑1 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGt                                                                                                                                >  4:410350/1‑150 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTa                                                                                                                               <  3:1238119/151‑1 (MQ=255)
                       gAAGGTTGCGGTTGCAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTa                                                                                                                               <  3:76064/151‑1 (MQ=255)
                              gcggttgcAACGACCTGACGTTCTGCGGCCCAGAATCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGc                                                                                                                        <  2:171646/151‑1 (MQ=255)
                               cggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATGCCGACTCGCGCCAGATCATAACGTACTAACGc                                                                                                                        >  6:896930/1‑151 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCAt                                                                                                                      >  8:1219349/1‑151 (MQ=255)
                                ggttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCAt                                                                                                                      <  5:462115/151‑1 (MQ=255)
                                  ttgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGAGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                                                    >  8:1082707/1‑151 (MQ=255)
                                   tgcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCCGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATcc                                                                                                                    >  7:917521/1‑150 (MQ=255)
                                    gcAACGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCa                                                                                                                   <  2:966179/150‑1 (MQ=255)
                                      aaCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCGTCCACgc                                                                                                                >  3:199504/1‑150 (MQ=255)
                                      aaCGACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGa                                                                                                                <  4:1209546/151‑1 (MQ=255)
                                         gACCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGCTCATAACGTACTAACGCATCCACGACCg                                                                                                             >  5:636943/1‑151 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                            <  4:882676/151‑1 (MQ=255)
                                          aCCTGACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGt                                                                                                            >  1:1029314/1‑151 (MQ=255)
                                              gACGTTCTGCGCCCCATAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtata                                                                                                         <  8:535451/150‑1 (MQ=255)
                                              gACGTTCTGCGCCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGtat                                                                                                          <  6:1071070/149‑1 (MQ=255)
                                                      gcgcCCCTTAAGCAGCCATAAGTGTTAAAGCAGCTGGCGAGGTATTCACGATCGAACAGGTTGTTAACATGCAGGGCCACGTTGGAGCCAGCCATG‑GGACTCGCGCCAGATCATAAGTTACTAACGCATCCACGACCGTATAACTTccc                                                                                                  <  1:1100289/149‑1 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCATAAGTGTTAACGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCa                                                                                                 <  7:1195108/150‑1 (MQ=255)
                                                      gcgcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCa                                                                                                 <  5:371744/150‑1 (MQ=255)
                                                       cgcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACt                                                                                               >  7:469127/1‑151 (MQ=255)
                                                       cgcCCCAGAAGCAGCCATAAGTGTAAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTAAACATGCCGCGCCAGGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCAGCCACGACCGTATAACTTCCCACt                                                                                               <  7:71619/151‑1 (MQ=255)
                                                        gcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CTACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACtt                                                                                              <  8:582328/151‑1 (MQ=255)
                                                        gcCCCAGAAGCAGCCATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCAcgt                                                                                              >  6:810740/1‑149 (MQ=255)
                                                                 gcagcCATAAGTGTTAAAGCAGCTGGCTACTTATTCACTATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCAATG‑CGAGTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGt                                                                                     <  1:336239/151‑1 (MQ=255)
                                                                  cagcCATAAGTGTTAAAGCAGCTGGGGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGCGCCAACCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGtt                                                                                    <  5:847522/151‑1 (MQ=255)
                                                                   agcCATAAGTGTTAAAGCACCTGGCGACGAATTCACGATCGAACAGGTTGTAAAGATGAAGCGCCACGTTGGAGCCAGCGAT‑CCGACTGGGGCCAGATCTTAAGGAACTAACCCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTa                                                                                   <  3:163891/151‑1 (MQ=255)
                                                                     ccATAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGc                                                                                 >  6:463951/1‑151 (MQ=255)
                                                                       aTAAGTGTTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCg                                                                               >  5:1215128/1‑151 (MQ=255)
                                                                         aaGTGTTAAACCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCgg                                                                              <  8:1136805/150‑1 (MQ=255)
                                                                           gtgtTAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATc                                                                           <  8:964655/151‑1 (MQ=255)
                                                                             gtTAAAGCAGCTGGCGACGTATTCACGATCGTACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAc                                                                         <  8:1167129/151‑1 (MQ=255)
                                                                              ttAAAGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCAcc                                                                        >  7:585269/1‑151 (MQ=255)
                                                                                aaaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGACGGATCAACAt                                                                      >  4:202069/1‑151 (MQ=255)
                                                                                 atggagcTGGCTTCGTATTTACGTTCGAACATGTTGTTAACATGCAGCGCCACGTTGTAGCGAGCCATG‑CGACTCGCGCCATATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                      <  2:1120856/146‑1 (MQ=255)
                                                                                 aaGCAGCTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCAt                                                                      >  1:211725/1‑150 (MQ=255)
                                                                                    cagcTGGCGACGTATTCAGGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACGCGCGCCAGATCATAACGTACTAACGCATGCCCGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCCCCATaagt                                                                  >  6:1175285/1‑149 (MQ=255)
                                                                                    cagcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACt                                                                  >  8:261200/1‑151 (MQ=255)
                                                                                     cgcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTg                                                                 <  8:801608/150‑1 (MQ=255)
                                                                                      gcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                                >  5:257084/1‑151 (MQ=255)
                                                                                      gcTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                                <  2:134074/151‑1 (MQ=255)
                                                                                       cTGGCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTgg                                                                >  4:1265156/1‑150 (MQ=255)
                                                                                         ggCGACGTATTCCCGAGCGAACAGGGTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAg                                                              <  4:587732/150‑1 (MQ=255)
                                                                                         ggCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGa                                                               <  6:124259/149‑1 (MQ=255)
                                                                                         ggCGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAg                                                              <  7:998192/150‑1 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACCGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                           <  2:1149930/151‑1 (MQ=255)
                                                                                           cGACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCa                                                           <  1:1212802/151‑1 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATACCGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                                            >  6:778026/1‑149 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGcc                                                            >  1:1049108/1‑149 (MQ=255)
                                                                                            gACGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGGTAGCCGGATCACCATAACTGGAGcc                                                            >  8:1242886/1‑149 (MQ=255)
                                                                                              cGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGta                                                        <  4:995691/151‑1 (MQ=255)
                                                                                              cGTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGt                                                         >  8:388125/1‑150 (MQ=255)
                                                                                               gTATTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGtat                                                       >  7:668012/1‑151 (MQ=255)
                                                                                                 aTTCACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATaa                                                     <  1:794086/151‑1 (MQ=255)
                                                                                                    cACGATGGAACAGGTTGTTAACATGCAGCGCCAGGTTGCAGCCAGCCATG‑CGACTCGCGCCAGATAATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACg                                                   <  3:607570/150‑1 (MQ=255)
                                                                                                    cACGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCAAG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACg                                                   >  6:584448/1‑150 (MQ=255)
                                                                                                     aCGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGCCCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGac                                                 >  5:53729/1‑151 (MQ=255)
                                                                                                      cGATCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGacc                                                <  3:657395/151‑1 (MQ=255)
                                                                                                        aTCGAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGcccac                                              >  5:95772/1‑151 (MQ=255)
                                                                                                           gAACAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCgg                                           <  4:510687/151‑1 (MQ=255)
                                                                                                            aaCAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCAAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGt                                          <  8:196688/151‑1 (MQ=255)
                                                                                                              cAGGTTGTTAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGc                                        <  8:227057/151‑1 (MQ=255)
                                                                                                                 gttgttAACATGCAGCGCCACGTTGTAGCCAGCCATG‑CGACTCGCGCCACATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCa                                     <  8:471266/151‑1 (MQ=255)
                                                                                                                 gttgttAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGccc                                      >  5:111815/1‑150 (MQ=255)
                                                                                                                      tAACATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                                >  5:1055528/1‑151 (MQ=255)
                                                                                                                       aaCATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCa                               <  6:367834/151‑1 (MQ=255)
                                                                                                                         cATGCAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTc                                >  4:94680/1‑148 (MQ=255)
                                                                                                                             cAGCGCCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTg                         <  6:521909/151‑1 (MQ=255)
                                                                                                                                cgcCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGaaa                      <  7:176947/151‑1 (MQ=255)
                                                                                                                                 gcCACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     <  8:1288817/151‑1 (MQ=255)
                                                                                                                                 gcCACGTGGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     <  7:238840/151‑1 (MQ=255)
                                                                                                                                 cccACGGTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     <  8:856312/150‑1 (MQ=255)
                                                                                                                                  ccACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGc                    >  1:600645/1‑151 (MQ=255)
                                                                                                                                  ccACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGc                    >  7:319527/1‑151 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAg                     >  2:419829/1‑149 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg                   >  3:969731/1‑151 (MQ=255)
                                                                                                                                   cACGTTGGAGCCAGCCAT‑CCGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCg                   <  6:247694/151‑1 (MQ=255)
                                                                                                                                         ggAGCCAGCCATG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCg              <  2:142153/150‑1 (MQ=255)
                                                                                                                                             ccagccaTG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaaa         <  2:226723/151‑1 (MQ=255)
                                                                                                                                             ccagccaTG‑CGACTCGCGCCAGATCATAACGTACTAACGCATCCACGACCGTATAACTTCCCACTTTAAAGGAGTTAGCCGGATCACCATAACTGGAGCCAGTATAACGACCACCGGTGCCCAGCGTCAGACCTGAAAGCGGACCGTCaaa         <  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: ‑