Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,341,704 Δ1 bp pseudogene (664/957 nt) PA14_RS10890 → DUF58 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,341,7020C.97.5% 2270.3 / 48.2 487pseudogene (662/957 nt)PA14_RS10890DUF58 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (8/4);  new base . (229/243);  total (238/249)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 507 total aligned reads displayed.

GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGGG  >  NC_008463/2341550‑2341849
                                                                                                                                                        |                                                                                                                                                   
   cTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccag                                                                                                                                                  <  2:2992312/151‑3 (MQ=255)
    tGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTCCCAGGGACGGCGCGCCTAccagc                                                                                                                                                 >  2:2256228/1‑148 (MQ=255)
    tGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagc                                                                                                                                                 <  1:7375774/151‑4 (MQ=255)
    tGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagc                                                                                                                                                 >  1:6394217/1‑148 (MQ=255)
      ggTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCg                                                                                                                                              <  2:4912366/151‑1 (MQ=255)
      ggTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCg                                                                                                                                              <  1:4701120/151‑1 (MQ=255)
       gTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCgg                                                                                                                                             >  1:9352392/1‑151 (MQ=255)
          gACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCg                                                                                                                                          >  2:3389497/1‑151 (MQ=255)
             cTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACt                                                                                                                                       >  1:4437320/1‑151 (MQ=255)
               tgcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTc                                                                                                                                      <  1:1506285/147‑1 (MQ=255)
                 cgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGAGGAGGATGCCGCGGACCTGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTCCCAGGGACGGCGCGCCTACCAGCCGGGCGACTcgc                                                                                                                                    >  1:880270/3‑151 (MQ=255)
                    cagcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggc                                                                                                                                 <  1:10006481/151‑1 (MQ=255)
                       cagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggc                                                                                                                              <  1:1760075/151‑1 (MQ=255)
                           cgcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCt                                                                                                                          <  1:5559556/151‑1 (MQ=255)
                           cgcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCt                                                                                                                          <  1:11167847/151‑1 (MQ=255)
                           cgcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCt                                                                                                                          <  2:5751527/151‑1 (MQ=255)
                            gcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTg                                                                                                                         >  2:708267/1‑151 (MQ=255)
                            gcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTg                                                                                                                         >  2:11002164/1‑151 (MQ=255)
                             cTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGc                                                                                                                        >  2:9083000/1‑151 (MQ=255)
                                  cTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTgg                                                                                                                   >  2:5628747/1‑151 (MQ=255)
                                   tACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGa                                                                                                                  <  1:6328732/151‑1 (MQ=255)
                                      ccTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAgg                                                                                                               >  2:8357276/1‑151 (MQ=255)
                                      ccTCGGCCGCTGGCCGGCGAGCCGTTACTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGGACGGGCGTCGACGACTACCAGGGACGGAGCGCCGACCAGCCCGGCGACTCGCGGCGGCGCCTGCACTGGAAgg                                                                                                               >  2:4768921/1‑151 (MQ=255)
                                       cTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGc                                                                                                              >  1:7681046/1‑151 (MQ=255)
                                        tCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCg                                                                                                             >  1:2505687/1‑151 (MQ=255)
                                        tCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTCCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCg                                                                                                             >  1:11046344/1‑151 (MQ=255)
                                         cGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGGCGCCTACCAGGGACTGCGCGCCTCCCAGCCGGGCGCCTCGCGGCGGCGCCTGCACTGGAGGGCGt                                                                                                            >  2:1684284/1‑151 (MQ=255)
                                               cTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTccc                                                                                                      <  2:8231231/151‑1 (MQ=255)
                                                tGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCg                                                                                                     >  2:8836547/1‑151 (MQ=255)
                                                tGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCg                                                                                                     <  1:3390313/151‑1 (MQ=255)
                                                 ggCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                    >  1:1419936/1‑151 (MQ=255)
                                                    cGGCGAGCCGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                 <  1:4847705/151‑1 (MQ=255)
                                                    cGGCGAGCCGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                 <  2:3145336/151‑1 (MQ=255)
                                                     ggCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGcc                                                                                                <  1:4308017/151‑1 (MQ=255)
                                                      gCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCa                                                                                               <  1:9159793/151‑1 (MQ=255)
                                                         gagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCAAGGGCGTCGACGACTACCAGGGAATGCGCGCCTACCAGCCGGGAGACTCGCGGCGGCGCCTGAACTGGAAGGCGTTTTCCGGTGGCCAgg                                                                                             >  2:2165197/4‑151 (MQ=255)
                                                           gccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAgggg                                                                                           >  1:8036324/2‑151 (MQ=255)
                                                             cGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgct                                                                                         >  1:6044700/1‑151 (MQ=255)
                                                              gTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctg                                                                                        <  2:9376062/151‑1 (MQ=255)
                                                                tCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctgct                                                                                      <  1:10793968/151‑1 (MQ=255)
                                                                  cTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTgg                                                                                    >  1:3165599/1‑151 (MQ=255)
                                                                  cTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTgg                                                                                    <  2:4242426/151‑1 (MQ=255)
                                                                  cTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTgg                                                                                    <  1:2860127/151‑1 (MQ=255)
                                                                   tCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGt                                                                                   <  1:5021258/151‑1 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACCACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTAACCCGAGGCCAGGGGCTGCTGGTGaa                                                                                >  2:3560206/1‑151 (MQ=255)
                                                                       gcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAg                                                                               >  1:6852913/1‑151 (MQ=255)
                                                                       gcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAg                                                                               <  1:2748743/151‑1 (MQ=255)
                                                                        ccGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCGACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAgg                                                                              <  2:5151589/151‑1 (MQ=255)
                                                                          gggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGAc                                                                            <  2:3971128/151‑1 (MQ=255)
                                                                           ggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACCGGCGACGCCGACTAACAGGGAATCCGCGCCTACCAGCCGGGCGACTCGCGGAGGCGCCTGCACTGGAAGGCGTTAACCAGAGGCCAGAGGATGCTCGTGAAGGACt                                                                           >  2:4987931/1‑151 (MQ=255)
                                                                            gcgcGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACtt                                                                          <  1:9915393/151‑1 (MQ=255)
                                                                            gcgcGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACtt                                                                          <  1:2048532/151‑1 (MQ=255)
                                                                             cgcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTc                                                                         <  2:4488652/151‑1 (MQ=255)
                                                                             cgcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTc                                                                         <  1:2211662/151‑1 (MQ=255)
                                                                               cgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgc                                                                       >  2:1339480/1‑151 (MQ=255)
                                                                                 gaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccg                                                                     >  1:115965/1‑151 (MQ=255)
                                                                                 gaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccg                                                                     >  1:10388225/1‑151 (MQ=255)
                                                                                    gaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCgcgc                                                                  >  2:9727483/1‑151 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                                <  2:11147993/151‑1 (MQ=255)
                                                                                        ccGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcgg                                                              >  2:2018934/1‑151 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                             <  1:1274698/151‑1 (MQ=255)
                                                                                          gAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggcg                                                            >  2:2535884/1‑151 (MQ=255)
                                                                                           aGGACCAGGGGCGCGCGCTGGCGAGCAAGGGCGACGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTACCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggcgg                                                           >  2:10267299/1‑151 (MQ=255)
                                                                                            ggACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCggg                                                          >  2:4144197/1‑151 (MQ=255)
                                                                                              aCCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGa                                                        >  1:2346071/1‑151 (MQ=255)
                                                                                              aCCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTCCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGa                                                        <  1:4512490/151‑1 (MQ=255)
                                                                                                cAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAg                                                      <  1:7224043/151‑1 (MQ=255)
                                                                                                  ggggCGCGCGCTGGCGAGCACGGGCGACGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGAc                                                    >  1:9711545/1‑151 (MQ=255)
                                                                                                    ggCGCGCGCTGGCGAGCACGGGCGTCGTCGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGTACTTCGCCGCGCTCGGCGGGGAAGACCt                                                  <  1:5369500/151‑1 (MQ=255)
                                                                                                    ggCGCGCGCTGGCGAGCACGGGCGTCGCCCACTACCAGGGACGGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACGGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCGGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGAccc                                                  <  1:7961517/151‑2 (MQ=255)
                                                                                                       gcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTGCCGTGGCCGGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCtgtg                                               >  2:7348324/1‑151 (MQ=255)
                                                                                                       gcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtgtg                                               <  1:1868731/151‑1 (MQ=255)
                                                                                                        cgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTgg                                              <  2:1188746/151‑1 (MQ=255)
                                                                                                         gcgcTGGCGAGCACGGGCGTCGACGTCTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCGGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGc                                             >  1:6503432/1‑151 (MQ=255)
                                                                                                           gcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTg                                           <  1:4080505/151‑1 (MQ=255)
                                                                                                                cGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAtt                                      <  2:10582005/151‑1 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                     <  2:6105095/151‑1 (MQ=255)
                                                                                                                  aGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCg                                    <  1:2994360/151‑1 (MQ=255)
                                                                                                                  aGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCg                                    >  2:9205112/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCGGTGGCTGGAATTCGAgc                                 >  2:4518876/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                 >  2:2816683/1‑151 (MQ=255)
                                                                                                                            tCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcgg                          <  2:6374709/151‑1 (MQ=255)
                                                                                                                            tCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcgg                          <  1:5939309/151‑1 (MQ=255)
                                                                                                                            tCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcgg                          <  2:1548152/151‑1 (MQ=255)
                                                                                                                             cgacgaCTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggc                         >  1:9038188/1‑151 (MQ=255)
                                                                                                                             cgacgaCTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggc                         <  2:7781317/151‑1 (MQ=255)
                                                                                                                             cgacgaCTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggc                         <  1:8532729/151‑1 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                      >  1:6754967/1‑151 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTGTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                      <  2:265289/151‑1 (MQ=255)
                                                                                                                                  acTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCGGTGGCTGGAATTCGAGCGCCTCGGCGGGGa                    >  2:5985107/1‑151 (MQ=255)
                                                                                                                                   cTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAt                   <  2:8565748/151‑1 (MQ=255)
                                                                                                                                       cAGGGACTGCGCGCCGACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCg               <  1:2610195/151‑1 (MQ=255)
                                                                                                                                          ggACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaac            <  1:375098/151‑2 (MQ=255)
                                                                                                                                           gACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAg           <  1:6663142/151‑1 (MQ=255)
                                                                                                                                             cTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTc         <  1:8469473/151‑1 (MQ=255)
                                                                                                                                             cTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTc         <  2:8465970/151‑1 (MQ=255)
                                                                                                                                             cTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTc         >  1:7029053/1‑151 (MQ=255)
                                                                                                                                                cgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGACAGTCGcc      >  2:6590692/1‑151 (MQ=255)
                                                                                                                                                 gcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCt     <  2:7006345/151‑1 (MQ=255)
                                                                                                                                                   gcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTgg   >  1:10171361/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGGG  >  NC_008463/2341550‑2341849

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