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.98.4% 2368.8 / 25.7 497pseudogene (662/957 nt)PA14_RS10890DUF58 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base . (244/242);  total (252/245)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 511 total aligned reads displayed.

GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848
                                                                                                                                                        |                                                                                                                                                  
gAGCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAc                                                                                                                                                    >  2:10276532/1‑151 (MQ=255)
  gCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAcca                                                                                                                                                  <  1:6108157/151‑2 (MQ=255)
       gTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCgg                                                                                                                                            <  1:9697850/151‑1 (MQ=255)
        tCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCggg                                                                                                                                           <  2:6424379/151‑1 (MQ=255)
             cTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACt                                                                                                                                      <  1:1761308/151‑1 (MQ=255)
                gcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTcg                                                                                                                                    >  2:4494455/4‑151 (MQ=255)
                 cgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTcgc                                                                                                                                   >  2:8507429/3‑151 (MQ=255)
                 cgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTcgc                                                                                                                                   >  1:222235/3‑151 (MQ=255)
                    cagcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggc                                                                                                                                <  2:8434374/151‑1 (MQ=255)
                     agcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcg                                                                                                                               <  2:10725128/151‑1 (MQ=255)
                     agcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcg                                                                                                                               >  2:5914737/1‑151 (MQ=255)
                      gcagGCGCTGGTCTACCCTCGGCCGCTGGACGGCGAGCCGTTCCTCGCGCCTGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCTCGCCTACCAGCCGGGCGACTCgcggcgg                                                                                                                              <  2:6663332/151‑1 (MQ=255)
                       cagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggc                                                                                                                             <  1:1481574/151‑1 (MQ=255)
                       cagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggc                                                                                                                             <  1:2821915/151‑1 (MQ=255)
                        agGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggcg                                                                                                                            >  2:7210749/1‑151 (MQ=255)
                          gcgcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGcc                                                                                                                          <  1:7806993/151‑1 (MQ=255)
                             cTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACCGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGc                                                                                                                       >  2:9147360/1‑151 (MQ=255)
                              tGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCa                                                                                                                      >  2:5157651/1‑151 (MQ=255)
                               ggTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCAc                                                                                                                     <  2:10344968/151‑1 (MQ=255)
                               cgTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGAGGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCGCCGGGCGACTCGCGGCGGCGCCTGCAc                                                                                                                     <  2:9662610/150‑1 (MQ=255)
                                   tACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGCCTCGCGGCGGCGCCTGCACTGGa                                                                                                                 >  2:72732/1‑151 (MQ=255)
                                    aCCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGaa                                                                                                                >  1:10548560/1‑151 (MQ=255)
                                      ccTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAgg                                                                                                              <  1:8235567/151‑1 (MQ=255)
                                      ccTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAgg                                                                                                              <  1:6563087/151‑1 (MQ=255)
                                        tCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCg                                                                                                            >  1:8467276/1‑151 (MQ=255)
                                        tCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCg                                                                                                            >  1:813169/1‑151 (MQ=255)
                                             cGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTc                                                                                                       >  2:11028197/1‑151 (MQ=255)
                                             cGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTc                                                                                                       >  1:3323668/1‑151 (MQ=255)
                                              gCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGAGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTcc                                                                                                      >  2:9561007/1‑151 (MQ=255)
                                              gCTGGCCGGCGAGCCGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTcc                                                                                                      <  1:3669948/151‑1 (MQ=255)
                                                tGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCg                                                                                                    >  1:1952437/1‑151 (MQ=255)
                                                tGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCg                                                                                                    >  1:767581/1‑151 (MQ=255)
                                                   ccGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTgg                                                                                                 <  2:9280550/151‑1 (MQ=255)
                                                    cGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                >  2:1991629/1‑151 (MQ=255)
                                                        cgagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAg                                                                                             <  2:1447174/147‑1 (MQ=255)
                                                         gagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAgg                                                                                            >  1:7268845/4‑151 (MQ=255)
                                                          agccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAggg                                                                                           >  1:10780192/3‑151 (MQ=255)
                                                            ccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgc                                                                                         <  1:1621924/151‑1 (MQ=255)
                                                             cGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgct                                                                                        >  1:8689539/1‑151 (MQ=255)
                                                               ttCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctgc                                                                                      >  2:1323861/1‑151 (MQ=255)
                                                                 ccTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGctgctg                                                                                    >  2:2193393/1‑151 (MQ=255)
                                                                 ccTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGctgctg                                                                                    >  1:10936789/1‑151 (MQ=255)
                                                                          gggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGAc                                                                           >  2:5720946/1‑151 (MQ=255)
                                                                            gcgcGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACtt                                                                         <  2:1500024/151‑1 (MQ=255)
                                                                              gcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcg                                                                       >  2:6144042/1‑151 (MQ=255)
                                                                              gcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcg                                                                       <  2:10851453/151‑1 (MQ=255)
                                                                                ggaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgcc                                                                     >  1:1267588/1‑151 (MQ=255)
                                                                                 gaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccg                                                                    <  1:10105492/151‑1 (MQ=255)
                                                                                 gaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccg                                                                    <  1:10396033/151‑1 (MQ=255)
                                                                                   ggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCcgcg                                                                  <  2:3635769/151‑1 (MQ=255)
                                                                                         cGAGGTCCAGGGGCGCGCGCTGGCGATCACGGGCGTCGTCGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGACTGCACTGGAAGGCGTTTTCCCGTGGACAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            <  1:7517356/151‑1 (MQ=255)
                                                                                          gAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGACGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggcg                                                           <  2:4258927/151‑1 (MQ=255)
                                                                                           aGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGTCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggcgg                                                          <  1:5035448/151‑1 (MQ=255)
                                                                                             gACCAGGGGCGCGCGCGGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCgggg                                                        <  2:4432831/151‑1 (MQ=255)
                                                                                              aCCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGa                                                       <  1:9461204/151‑1 (MQ=255)
                                                                                               ccAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGaa                                                      <  1:3120479/151‑1 (MQ=255)
                                                                                               ccACGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGaa                                                      <  2:1060493/151‑1 (MQ=255)
                                                                                                 aGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGa                                                    >  2:743405/1‑151 (MQ=255)
                                                                                                 aGGGGCGCGCGCTCGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGa                                                    <  2:9066460/151‑1 (MQ=255)
                                                                                                  ggggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGAc                                                   <  2:10090939/151‑1 (MQ=255)
                                                                                                  ggggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGAc                                                   <  1:3957821/151‑1 (MQ=255)
                                                                                                    ggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCGACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCt                                                 <  1:5752506/151‑1 (MQ=255)
                                                                                                     gcgcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtg                                                <  1:10603563/151‑1 (MQ=255)
                                                                                                       gcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtgtg                                              <  1:2490057/151‑1 (MQ=255)
                                                                                                        cgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTgg                                             <  1:422264/151‑1 (MQ=255)
                                                                                                        cgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTgg                                             >  2:8622624/1‑151 (MQ=255)
                                                                                                         gcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGc                                            <  2:4871041/151‑1 (MQ=255)
                                                                                                           gcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTg                                          >  2:7636500/1‑151 (MQ=255)
                                                                                                           gcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTGTCCCGGGGCCAGGGGCTGCTGGTGAAGGACTGCGCCGCGCTCGGCGGGGACGACCTGTGGCTg                                          >  2:808635/1‑151 (MQ=255)
                                                                                                            cTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTgg                                         >  1:11113366/1‑151 (MQ=255)
                                                                                                              ggCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGaa                                       >  2:2544247/1‑151 (MQ=255)
                                                                                                               gCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAt                                      >  1:6321332/1‑151 (MQ=255)
                                                                                                                cGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAtt                                     <  2:9775508/151‑1 (MQ=255)
                                                                                                                cGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAtt                                     <  2:9384816/151‑1 (MQ=255)
                                                                                                                 gAGCACGGGAGTCGACGACTACCAGGGACTTCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGTAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAGTc                                    <  1:6979616/151‑1 (MQ=255)
                                                                                                                    cACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCTGTGGCTGGAATTCGAg                                 >  2:2872887/1‑151 (MQ=255)
                                                                                                                    cACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAg                                 >  1:929430/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                >  2:5996366/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                >  2:3724595/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                <  2:2807529/151‑1 (MQ=255)
                                                                                                                      cGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgcg                               >  2:4968368/1‑151 (MQ=255)
                                                                                                                      cGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgcg                               >  1:10274411/1‑151 (MQ=255)
                                                                                                                         gCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCt                            <  2:8332803/151‑1 (MQ=255)
                                                                                                                          cGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTc                           >  2:4065796/1‑151 (MQ=255)
                                                                                                                          cGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGGTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGACTTCGAGCGCCTc                           >  1:9065291/1‑151 (MQ=255)
                                                                                                                           gTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcg                          >  2:6941290/1‑151 (MQ=255)
                                                                                                                            tCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcgg                         >  1:4616804/1‑151 (MQ=255)
                                                                                                                             cgacgaCTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggc                        <  2:3128362/151‑1 (MQ=255)
                                                                                                                              gacgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcg                       >  2:7528160/1‑151 (MQ=255)
                                                                                                                               acgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcgg                      >  2:119012/1‑151 (MQ=255)
                                                                                                                                cgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                     >  1:13924/1‑151 (MQ=255)
                                                                                                                                 gacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCgggg                    >  2:3998718/1‑151 (MQ=255)
                                                                                                                                      ccAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAtctc               >  2:3382365/1‑151 (MQ=255)
                                                                                                                                        aGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGa             >  2:3092491/1‑151 (MQ=255)
                                                                                                                                         gggACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaa            >  2:9163680/1‑151 (MQ=255)
                                                                                                                                         gggACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaa            >  1:5471104/1‑151 (MQ=255)
                                                                                                                                             cTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTc        <  1:6851064/151‑1 (MQ=255)
                                                                                                                                              tGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCg       <  2:459654/151‑1 (MQ=255)
                                                                                                                                               gcgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGc      >  2:10483885/1‑151 (MQ=255)
                                                                                                                                               gcgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGc      <  2:2080403/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848

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