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.5% 2306.4 / 17.9 480pseudogene (662/957 nt)PA14_RS10890DUF58 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (6/1);  new base . (239/233);  total (246/234)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 503 total aligned reads displayed.

GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGGGC  >  NC_008463/2341550‑2341850
                                                                                                                                                        |                                                                                                                                                    
gAGCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCGGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAc                                                                                                                                                      >  1:8567148/1‑151 (MQ=255)
  gCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAcca                                                                                                                                                    <  2:7360142/151‑2 (MQ=255)
  gCTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAcca                                                                                                                                                    <  1:4623591/151‑2 (MQ=255)
   cTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccag                                                                                                                                                   <  1:4466323/151‑3 (MQ=255)
   cTGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGAGGCCGAGGACCAGGGGCGCGCGCTGGCGGGCACGGGCGTCGACGACTACCAGGGCCTGCGCGCCTAccgg                                                                                                                                                   >  1:8015416/1‑149 (MQ=255)
    tGGGTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagc                                                                                                                                                  <  2:7780208/151‑4 (MQ=255)
     gggTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagcc                                                                                                                                                 <  1:7130461/151‑5 (MQ=255)
      ggTCGACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCg                                                                                                                                               <  1:7373608/151‑1 (MQ=255)
          gACCTGCGCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCg                                                                                                                                           <  1:5938792/151‑1 (MQ=255)
                    cagcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCTTCCTCGCGCCGGGCGCGCAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggc                                                                                                                                  <  1:1148062/151‑1 (MQ=255)
                    cagcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggc                                                                                                                                  <  1:6452034/151‑1 (MQ=255)
                      gcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcgg                                                                                                                                >  1:2230467/1‑151 (MQ=255)
                       cagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggc                                                                                                                               >  2:5132194/1‑151 (MQ=255)
                        agGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCACGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCgcggcggcg                                                                                                                              >  2:828770/1‑151 (MQ=255)
                          gcgcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGcc                                                                                                                            <  1:2880422/151‑1 (MQ=255)
                            gcTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTg                                                                                                                          <  2:6955272/151‑1 (MQ=255)
                             cTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGc                                                                                                                         >  1:7498104/1‑151 (MQ=255)
                               ggTCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCAc                                                                                                                       <  1:8938248/151‑1 (MQ=255)
                                 tCTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGTGGACCCGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCCGGGACTGCGCGCCTCCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTg                                                                                                                     >  2:1127354/1‑151 (MQ=255)
                                  cTACCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTgg                                                                                                                    <  1:1990137/151‑1 (MQ=255)
                                    aCCCTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGaa                                                                                                                  >  2:181177/1‑151 (MQ=255)
                                       cTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGc                                                                                                               >  2:2844970/1‑151 (MQ=255)
                                       cTCGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGc                                                                                                               >  1:9089542/1‑151 (MQ=255)
                                         cGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGt                                                                                                             <  1:7004523/151‑1 (MQ=255)
                                         cGGCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGt                                                                                                             <  2:3758141/151‑1 (MQ=255)
                                           gCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGttt                                                                                                           >  1:5491796/1‑151 (MQ=255)
                                           gCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGttt                                                                                                           >  2:5409874/1‑151 (MQ=255)
                                           gCCGCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGttt                                                                                                           <  2:5711090/151‑1 (MQ=255)
                                              gCTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTcc                                                                                                        >  2:3391003/1‑151 (MQ=255)
                                               cTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTccc                                                                                                       <  1:1000543/151‑1 (MQ=255)
                                               cTGGCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACGGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTccc                                                                                                       >  1:3823907/1‑151 (MQ=255)
                                                 ggCCGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                     <  2:2424752/151‑1 (MQ=255)
                                                   ccGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTgg                                                                                                   >  2:4636406/1‑151 (MQ=255)
                                                   ccGGCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTgg                                                                                                   <  2:750047/151‑1 (MQ=255)
                                                      gCGAGCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCa                                                                                                <  2:885534/151‑1 (MQ=255)
                                                        cgagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAg                                                                                               >  1:10816163/5‑151 (MQ=255)
                                                         gagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAgg                                                                                              >  1:5829792/4‑151 (MQ=255)
                                                            ccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgc                                                                                           >  2:3637427/1‑151 (MQ=255)
                                                            ccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgc                                                                                           >  2:7017580/1‑151 (MQ=255)
                                                            ccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgc                                                                                           >  2:2016987/1‑151 (MQ=255)
                                                             cGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgct                                                                                          <  2:8920416/151‑1 (MQ=255)
                                                              gTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctg                                                                                         >  1:6534422/1‑151 (MQ=255)
                                                               gtCCTCGCGCCGGGCGCGGAGGATGCCGCGGACCCGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTGGTCCCGTGGCCAGGGgctgc                                                                                        <  2:10010293/150‑1 (MQ=255)
                                                                tCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctgct                                                                                       <  2:9532198/151‑1 (MQ=255)
                                                                 ccTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGctgctg                                                                                      <  2:8531869/151‑1 (MQ=255)
                                                                     gcgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGa                                                                                  <  2:9311233/151‑1 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGaa                                                                                 >  1:6227550/1‑151 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGaa                                                                                 <  1:1434790/151‑1 (MQ=255)
                                                                           ggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACt                                                                            >  1:4753677/1‑151 (MQ=255)
                                                                              gcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcg                                                                         <  2:3963800/151‑1 (MQ=255)
                                                                              gcgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcg                                                                         >  2:10437330/1‑151 (MQ=255)
                                                                               cgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgc                                                                        <  1:4240195/151‑1 (MQ=255)
                                                                                ggaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGACGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgcc                                                                       >  1:8167862/1‑151 (MQ=255)
                                                                                     aTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCt                                                                  >  1:2540840/1‑151 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                                 >  1:8404446/1‑151 (MQ=255)
                                                                                       gCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcg                                                                <  2:4201098/151‑1 (MQ=255)
                                                                                        ccGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcgg                                                               >  1:10342280/1‑151 (MQ=255)
                                                                                            ggACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCggg                                                           >  1:5304549/1‑151 (MQ=255)
                                                                                              aCCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGa                                                         >  2:10103166/1‑151 (MQ=255)
                                                                                               ccAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGaa                                                        <  1:321096/151‑1 (MQ=255)
                                                                                                cAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAg                                                       >  1:10562766/1‑151 (MQ=255)
                                                                                                   gggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGAcc                                                    >  1:776530/1‑151 (MQ=255)
                                                                                                    ggCGCGCGCGGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCt                                                   <  1:4658948/151‑1 (MQ=255)
                                                                                                       gcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtgtg                                                >  1:1022217/1‑151 (MQ=255)
                                                                                                           gcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTg                                            >  2:4271395/1‑151 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                          >  1:2808582/1‑151 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                          >  2:10181309/1‑151 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                          >  1:3472906/1‑151 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCGGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                      >  1:2645414/1‑151 (MQ=255)
                                                                                                                   gCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGa                                    <  1:2368172/151‑1 (MQ=255)
                                                                                                                   gCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGa                                    <  2:9399647/151‑1 (MQ=255)
                                                                                                                    cACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAg                                   >  1:6669842/1‑151 (MQ=255)
                                                                                                                        ggCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGcc                               <  2:6145303/151‑1 (MQ=255)
                                                                                                                         gCGTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCt                              >  2:11088674/1‑151 (MQ=255)
                                                                                                                           gTCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcg                            <  1:9290859/151‑1 (MQ=255)
                                                                                                                            tCGACGACTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcgg                           <  1:9613738/151‑1 (MQ=255)
                                                                                                                             cgacgaCTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGTAGGCGTGTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggc                          <  1:4868678/151‑1 (MQ=255)
                                                                                                                              gacgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcg                         <  2:3030010/151‑1 (MQ=255)
                                                                                                                              gacgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcg                         >  2:6714247/1‑151 (MQ=255)
                                                                                                                              gacgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcg                         <  2:2163401/151‑1 (MQ=255)
                                                                                                                               acgacTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcgg                        <  2:1596274/151‑1 (MQ=255)
                                                                                                                                   cgaCCAGGGACCGCGCGCCCCCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAt                    <  2:7095118/149‑1 (MQ=255)
                                                                                                                                   cTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAt                    >  1:3539275/1‑151 (MQ=255)
                                                                                                                                   cTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAt                    >  1:989746/1‑151 (MQ=255)
                                                                                                                                   cTACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAt                    >  1:1636569/1‑151 (MQ=255)
                                                                                                                                    tACCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAtc                   >  2:10818172/1‑151 (MQ=255)
                                                                                                                                     aCCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAtct                  >  1:8773773/1‑151 (MQ=255)
                                                                                                                                     aCCAGGGACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAtct                  >  2:2536478/1‑151 (MQ=255)
                                                                                                                                          ggACTGCGCGCCTACCTGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaaa             <  2:4457080/151‑1 (MQ=255)
                                                                                                                                          ggACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaaa             <  1:5658545/151‑1 (MQ=255)
                                                                                                                                          ggACTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaaa             >  2:8369760/1‑151 (MQ=255)
                                                                                                                                            aCTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGt           >  1:6827424/1‑151 (MQ=255)
                                                                                                                                             cTGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTc          >  1:9489488/1‑151 (MQ=255)
                                                                                                                                              tGCGCGCCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCg         >  1:3992805/1‑151 (MQ=255)
                                                                                                                                               gcgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGc        >  2:8069866/1‑151 (MQ=255)
                                                                                                                                               gcgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGc        <  2:2359620/151‑1 (MQ=255)
                                                                                                                                               gcgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGc        <  2:9708706/151‑1 (MQ=255)
                                                                                                                                                cgcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAATGTCGcc       >  2:1509182/1‑151 (MQ=255)
                                                                                                                                                 gcgcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCt      >  1:5113422/1‑151 (MQ=255)
                                                                                                                                                   gcCTACCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTgg    >  2:4897380/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGGGC  >  NC_008463/2341550‑2341850

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: