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.3% 2396.0 / 47.9 511pseudogene (662/957 nt)PA14_RS10890DUF58 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base . (237/260);  total (245/266)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 524 total aligned reads displayed.

GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848
                                                                                                                                                        |                                                                                                                                                  
 aGCTGGGTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCG‑GGCCGTTCCTCGCGCCGGGCGCGGCGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCGGGGACTGGGCGCCTAcc                                                                                                                                                   >  2:11322096/1‑151 (MQ=255)
   cTGGGTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccag                                                                                                                                                 <  1:4378959/151‑3 (MQ=255)
    tGGGTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCGCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagc                                                                                                                                                >  1:10869007/1‑148 (MQ=255)
     gggTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagcc                                                                                                                                               <  1:10769745/151‑5 (MQ=255)
      ggTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGTGCACGGGCGTCGACGTCTACCAGGGACTGCGCGCCTA‑CCAGCCg                                                                                                                                             >  1:10620192/1‑151 (MQ=255)
        tCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCggg                                                                                                                                           <  2:7002717/151‑1 (MQ=255)
        tCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACGGCGCGCCTA‑CCAGCCggg                                                                                                                                           >  2:8336645/1‑151 (MQ=255)
         cGACCTGCG‑CCAGCAGGCGCTGGTCTCCCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGc                                                                                                                                          <  2:9431034/151‑1 (MQ=255)
          gACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCg                                                                                                                                         >  1:531555/1‑151 (MQ=255)
             cTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACt                                                                                                                                      >  2:4586481/1‑151 (MQ=255)
               tgcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGCTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCC‑GCCGGGCGACTc                                                                                                                                     >  1:1558112/5‑151 (MQ=255)
                 cgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTcgc                                                                                                                                   <  1:5039207/149‑1 (MQ=255)
                  gccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTcgcg                                                                                                                                  >  1:3376094/2‑151 (MQ=255)
                  gccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTcgcg                                                                                                                                  <  2:5080486/150‑1 (MQ=255)
                     agcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCgcggcg                                                                                                                               >  1:8621852/1‑151 (MQ=255)
                     agcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCgcggcg                                                                                                                               >  2:6636308/1‑151 (MQ=255)
                     agcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCgcggcg                                                                                                                               <  2:3789573/151‑1 (MQ=255)
                     agcagGCGCTGGTCGACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCgcggcg                                                                                                                               <  1:11341592/151‑1 (MQ=255)
                           cgcTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCt                                                                                                                         <  1:5971935/151‑1 (MQ=255)
                             cTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGc                                                                                                                       <  2:2643579/151‑1 (MQ=255)
                             cTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGc                                                                                                                       >  1:2802967/1‑151 (MQ=255)
                                 tCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTg                                                                                                                   <  2:6091131/151‑1 (MQ=255)
                                      ccTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAgg                                                                                                              >  2:200755/1‑151 (MQ=255)
                                       cTCGGCCGCTGGCCGGCGAG‑CCGGTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCT‑CCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGc                                                                                                             <  1:10073877/151‑1 (MQ=255)
                                          ggCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGTAGGCGtt                                                                                                          >  1:5278027/1‑151 (MQ=255)
                                            ccGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtttt                                                                                                        <  2:9688762/151‑1 (MQ=255)
                                            ccGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtttt                                                                                                        >  1:10952873/1‑151 (MQ=255)
                                            ccGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtttt                                                                                                        <  2:9885267/151‑1 (MQ=255)
                                               cTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTccc                                                                                                     >  2:2219212/1‑151 (MQ=255)
                                               cTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTccc                                                                                                     >  1:2305108/1‑151 (MQ=255)
                                                tGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCCTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCg                                                                                                    >  1:956722/1‑151 (MQ=255)
                                                 ggCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                   <  2:2429129/151‑1 (MQ=255)
                                                 ggCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGt                                                                                                   <  1:9097804/151‑1 (MQ=255)
                                                  gCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTg                                                                                                  <  2:7637393/151‑1 (MQ=255)
                                                     ggCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGcc                                                                                               >  2:4254977/1‑151 (MQ=255)
                                                      gCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCa                                                                                              >  2:2313208/1‑151 (MQ=255)
                                                        cgagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCCGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAg                                                                                             <  2:9678386/147‑1 (MQ=255)
                                                         gagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAgg                                                                                            >  1:5493614/4‑151 (MQ=255)
                                                         gagccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAgg                                                                                            <  2:8995573/148‑1 (MQ=255)
                                                            ccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgc                                                                                         <  2:866040/151‑1 (MQ=255)
                                                              gTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctg                                                                                       <  2:5931128/151‑1 (MQ=255)
                                                                tCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgctgct                                                                                     >  2:4851143/1‑151 (MQ=255)
                                                                  cTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTgg                                                                                   >  1:10730361/1‑151 (MQ=255)
                                                                   tCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTggg                                                                                  >  2:3546300/1‑150 (MQ=255)
                                                                   tCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGt                                                                                  <  2:2968901/151‑1 (MQ=255)
                                                                   tCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGt                                                                                  >  1:3141155/1‑151 (MQ=255)
                                                                    cgcgCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGCGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCGA‑CCAGCCGGGCGCCTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTg                                                                                 <  2:1010893/151‑1 (MQ=255)
                                                                    cgcgCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTg                                                                                 >  1:6551983/1‑151 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGaa                                                                               >  1:10436418/1‑151 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGaa                                                                               <  2:6273385/151‑1 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGGCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGAGGCCAGGGGCGGCTGGTGaa                                                                               >  2:3751927/1‑151 (MQ=255)
                                                                        ccGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAgg                                                                             <  1:2886581/151‑1 (MQ=255)
                                                                        ccGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTGTCCCGTGGCCAGGGGCTGCTGGTGAAgg                                                                             >  1:3077180/1‑151 (MQ=255)
                                                                        ccGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCGA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAgg                                                                             <  2:3196816/151‑1 (MQ=255)
                                                                         cGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGa                                                                            >  2:571149/1‑151 (MQ=255)
                                                                           ggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACt                                                                          <  2:4422891/151‑1 (MQ=255)
                                                                            gcgcGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACtt                                                                         <  1:4147662/151‑1 (MQ=255)
                                                                               cgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgc                                                                      <  2:7846247/151‑1 (MQ=255)
                                                                               cgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgc                                                                      >  1:7961089/1‑151 (MQ=255)
                                                                               cgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgc                                                                      <  1:3953738/151‑1 (MQ=255)
                                                                                ggaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgcc                                                                     <  2:1073559/151‑1 (MQ=255)
                                                                                ggaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgcc                                                                     >  1:9629958/1‑151 (MQ=255)
                                                                                  aggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccgc                                                                   >  1:1723709/1‑151 (MQ=255)
                                                                                   ggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCcgcg                                                                  >  1:788667/1‑151 (MQ=255)
                                                                                   ggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCcgcg                                                                  <  2:9250296/151‑1 (MQ=255)
                                                                                   ggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCcgcg                                                                  <  1:1159103/151‑1 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                               <  2:1253566/151‑1 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            >  1:8179667/1‑151 (MQ=255)
                                                                                            ggACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCggg                                                         <  1:10479483/151‑1 (MQ=255)
                                                                                            ggACCAGGGGCGCGCGCTGGCGAGCACGGGCGACGACCACTACCAGGGACTGCGCGCCTA‑CCTGCCGGGCGACTCGCGGGGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGGGAAGGACTTCGCCGCGCTCGGCggg                                                         >  2:10996821/1‑151 (MQ=255)
                                                                                             gACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCgggg                                                        <  1:2604441/151‑1 (MQ=255)
                                                                                             gACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCgggg                                                        >  1:8366367/1‑151 (MQ=255)
                                                                                              aCCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGa                                                       >  1:5171730/1‑151 (MQ=255)
                                                                                               ccAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGaa                                                      >  1:6920607/1‑151 (MQ=255)
                                                                                                cAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTGTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAg                                                     <  2:7154309/151‑1 (MQ=255)
                                                                                                 aGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGa                                                    <  1:56236/151‑1 (MQ=255)
                                                                                                   gggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGAcc                                                  >  1:1019531/1‑151 (MQ=255)
                                                                                                    ggCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCt                                                 <  2:11248697/151‑1 (MQ=255)
                                                                                                         gcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGc                                            >  2:7318588/1‑151 (MQ=255)
                                                                                                           gcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTg                                          <  2:6838121/151‑1 (MQ=255)
                                                                                                            cTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTgg                                         >  2:531516/1‑151 (MQ=255)
                                                                                                            cTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTgg                                         <  1:7424981/151‑1 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                        <  2:1331954/151‑1 (MQ=255)
                                                                                                             tGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGa                                        >  1:584095/1‑151 (MQ=255)
                                                                                                               gCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCNACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAt                                      >  1:6118924/1‑151 (MQ=255)
                                                                                                               gCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGTAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAt                                      >  2:8547681/1‑151 (MQ=255)
                                                                                                               gCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAt                                      <  2:7420624/151‑1 (MQ=255)
                                                                                                                cGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAtt                                     <  1:7347770/151‑1 (MQ=255)
                                                                                                                  aGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCg                                   <  2:6743044/151‑1 (MQ=255)
                                                                                                                   gCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGa                                  >  1:983861/1‑151 (MQ=255)
                                                                                                                          cGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTc                           >  1:3791568/1‑151 (MQ=255)
                                                                                                                            tCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcgg                         >  1:8872818/1‑151 (MQ=255)
                                                                                                                             cgacgaCTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTGTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggc                        <  1:7593516/151‑1 (MQ=255)
                                                                                                                               acgacTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcgg                      >  1:665205/1‑151 (MQ=255)
                                                                                                                                 gacTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCgggg                    >  2:8136636/1‑151 (MQ=255)
                                                                                                                                  acTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGa                   <  1:2158944/151‑1 (MQ=255)
                                                                                                                                           gACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAg          <  2:10528172/151‑1 (MQ=255)
                                                                                                                                               gcgcgcCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGc      <  1:4729390/151‑1 (MQ=255)
                                                                                                                                                 gcgcCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCt    >  2:5465610/1‑151 (MQ=255)
                                                                                                                                                  cgcCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTg   >  2:2726453/1‑151 (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: ‑