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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,341,7020C.97.5% 2258.5 / 43.0 483pseudogene (662/957 nt)PA14_RS10890DUF58 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (5/7);  new base . (215/254);  total (220/263)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 503 total aligned reads displayed.

GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848
                                                                                                                                                        |                                                                                                                                                  
gAGCTGGGTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAc                                                                                                                                                    <  2:1594314/151‑1 (MQ=255)
  gCTGGGTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAcca                                                                                                                                                  >  1:7210817/1‑150 (MQ=255)
    tGGGTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTAccagc                                                                                                                                                >  1:6029504/1‑148 (MQ=255)
      ggTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGTGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGTCTGCGCGCCTA‑CCAGCCg                                                                                                                                             >  1:970926/1‑151 (MQ=255)
       gTCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCgg                                                                                                                                            >  1:2020432/1‑151 (MQ=255)
        tCGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCggg                                                                                                                                           <  2:2704237/151‑1 (MQ=255)
        tCGAACTGCG‑CCAGCAGGCGGTGGTCTACCCTCGGCCGCTTGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGTTGCCGAGGACAAGGGGCGCGCGCTGGCGAGCACGTGAGTCGACGACTTCAAGGGTCTGCGCTCCTA‑CAAGCCggg                                                                                                                                           <  1:4787594/151‑1 (MQ=255)
         cGACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGc                                                                                                                                          <  1:7112256/151‑1 (MQ=255)
          gACCTGCG‑CCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCg                                                                                                                                         <  1:2394667/151‑1 (MQ=255)
                gcgccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTcg                                                                                                                                    <  1:10768863/148‑1 (MQ=255)
                  gccAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTcgcg                                                                                                                                  >  2:7345037/2‑151 (MQ=255)
                     agcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCgcggcg                                                                                                                               >  1:388097/1‑151 (MQ=255)
                      gcagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCgcggcgg                                                                                                                              <  1:1738686/151‑1 (MQ=255)
                       cagGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCgcggcggc                                                                                                                             <  1:1811534/151‑1 (MQ=255)
                         ggCGCTTGTCTCCCCTCGGCCGCTGGCCGGCGA‑CCCGGACCTCGCGCCGGGCCCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCC‑GCCGGGCGACTCGCGGCGgcgc                                                                                                                           <  1:10587937/151‑1 (MQ=255)
                          gcgcTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGACGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGcc                                                                                                                          <  1:6172709/151‑1 (MQ=255)
                            gcTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTg                                                                                                                        <  1:7653467/151‑1 (MQ=255)
                            gcTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGAGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCTGGCGAATCGCGGCGGCGCCTg                                                                                                                        <  2:9374246/151‑1 (MQ=255)
                             cTGGTCTACCCTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGc                                                                                                                       <  1:7546581/151‑1 (MQ=255)
                                    aCCCTCGTCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGaa                                                                                                                <  2:2284920/151‑1 (MQ=255)
                                     ccctgcgCCGCTGTCCTGCGCG‑CCGTTCCTCTCGACCGGCCCGGCGGATGCCGCGGACCATGGGCGCGCGCTGGCTAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAg                                                                                                               <  2:10918885/145‑1 (MQ=255)
                                     cccTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAg                                                                                                               <  1:10340474/151‑1 (MQ=255)
                                     cccTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAg                                                                                                               >  2:4646282/1‑151 (MQ=255)
                                      ccTCGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAgg                                                                                                              <  1:4912230/151‑1 (MQ=255)
                                         cGGCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCCGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGt                                                                                                           >  1:3570708/1‑151 (MQ=255)
                                          ggCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtt                                                                                                          <  2:10433236/151‑1 (MQ=255)
                                          ggCCGCTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGtt                                                                                                          <  1:10414074/151‑1 (MQ=255)
                                               cTGGCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACCACTACCAGCGAATGCGCGCCTA‑CCAGACGGCAGTCTCGCGGCGGCGCCTGCACTGGAAGGCGttaaccc                                                                                                     >  2:9810026/1‑146 (MQ=255)
                                                  gCCGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTg                                                                                                  <  1:505011/151‑1 (MQ=255)
                                                   ccGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTgg                                                                                                 <  1:4250563/151‑1 (MQ=255)
                                                    cGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                >  1:4477164/1‑151 (MQ=255)
                                                    cGGCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGc                                                                                                <  2:8223390/151‑1 (MQ=255)
                                                     ggCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGcc                                                                                               >  2:8954319/1‑151 (MQ=255)
                                                      gCGAG‑CCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTGTCGCGTGGccg                                                                                              >  2:7820943/1‑150 (MQ=255)
                                                          agccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAggg                                                                                           >  1:6782105/3‑151 (MQ=255)
                                                          agccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAggg                                                                                           <  1:8429461/149‑1 (MQ=255)
                                                          agccGGTCCTCGCGCCGGGCGCGGAGGATGCCGTGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAggg                                                                                           <  2:1651195/149‑1 (MQ=255)
                                                            ccGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgc                                                                                         <  2:8431759/151‑1 (MQ=255)
                                                             cGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGgct                                                                                        >  1:1199815/1‑151 (MQ=255)
                                                                  cTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCggg                                                                                   >  1:10856558/1‑151 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGaa                                                                               <  1:3002893/151‑1 (MQ=255)
                                                                      cgcCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGaa                                                                               <  1:8628733/151‑1 (MQ=255)
                                                                         cGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGa                                                                            >  2:4243196/1‑151 (MQ=255)
                                                                         cGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGa                                                                            <  2:3189261/151‑1 (MQ=255)
                                                                         cGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGa                                                                            <  1:5674315/151‑1 (MQ=255)
                                                                          gggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGAc                                                                           <  2:6207773/151‑1 (MQ=255)
                                                                          gggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGAc                                                                           <  1:2841149/151‑1 (MQ=255)
                                                                           ggCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACt                                                                          >  2:5488736/1‑151 (MQ=255)
                                                                               cgGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgc                                                                      >  1:271437/1‑151 (MQ=255)
                                                                                ggaggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGTCTTcgcc                                                                     >  1:6959233/1‑151 (MQ=255)
                                                                                  aggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccgc                                                                   <  1:9904251/151‑1 (MQ=255)
                                                                                  aggaTGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTcgccgc                                                                   <  2:7064015/151‑1 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                               <  2:4542167/151‑1 (MQ=255)
                                                                                      tGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTc                                                               <  1:9227316/151‑1 (MQ=255)
                                                                                       gCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcg                                                              >  1:6674110/1‑151 (MQ=255)
                                                                                       gCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACGACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcg                                                              <  2:649421/151‑1 (MQ=255)
                                                                                        ccGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcgg                                                             <  1:10185968/151‑1 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            >  2:7137031/1‑151 (MQ=255)
                                                                                         cGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggc                                                            <  1:1585903/151‑1 (MQ=255)
                                                                                          gAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTcggcg                                                           >  1:526441/1‑151 (MQ=255)
                                                                                           aGGACCAGGGGCGCGCGCTTGCGAGCTCGGGCGTCGACGTCTACCAGGGACTGCGCTCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTTGAAGGCGTTTTCCCGTGGCCAGGTGCTGCTGGTGAAGGACTTCGCCGCGCTcggctg                                                          <  2:5585331/151‑3 (MQ=255)
                                                                                             gACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCgggg                                                        <  2:8830195/151‑1 (MQ=255)
                                                                                              aCCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGa                                                       >  1:4083690/1‑151 (MQ=255)
                                                                                                 aGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGAGACTCGCGGCGGCGCCTGCACTTTAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGaata                                                    <  2:1768750/151‑3 (MQ=255)
                                                                                                      cgcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtgt                                               >  2:3929954/1‑151 (MQ=255)
                                                                                                      cgcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtgt                                               >  2:3827755/1‑151 (MQ=255)
                                                                                                       gcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCtgtg                                              <  2:6643336/151‑1 (MQ=255)
                                                                                                       gcgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACATCGCCGCGCTCGGCGGGGAAGACCtgtg                                              <  1:6544769/151‑1 (MQ=255)
                                                                                                        cgcgcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTgg                                             >  2:4874239/1‑151 (MQ=255)
                                                                                                           gcTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTg                                          >  1:10087753/1‑151 (MQ=255)
                                                                                                              ggCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCC‑GCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGAGTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGaa                                       <  1:8815165/151‑1 (MQ=255)
                                                                                                                cGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAAtt                                     <  2:410685/151‑1 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                    >  1:4563820/1‑151 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                    <  2:5240047/151‑1 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                    <  1:7957048/151‑1 (MQ=255)
                                                                                                                 gAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTc                                    <  2:7712544/151‑1 (MQ=255)
                                                                                                                  aGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCg                                   <  2:7603429/151‑1 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGACCTGTGGCTGGAATTCGAgc                                >  2:3663959/1‑151 (MQ=255)
                                                                                                                     aCGGGCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgc                                <  1:5851815/151‑1 (MQ=255)
                                                                                                                       gggCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAgcgc                              >  2:2195534/1‑151 (MQ=255)
                                                                                                                         gCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGACGTCCTGGGGCTGGAATTCGAGCGCCt                            >  1:577955/1‑151 (MQ=255)
                                                                                                                         gCGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCt                            <  2:10527124/151‑1 (MQ=255)
                                                                                                                          cGTCGACGACTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTc                           <  1:8299873/151‑1 (MQ=255)
                                                                                                                             cgacgaCTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggc                        >  2:7466256/1‑151 (MQ=255)
                                                                                                                               acgacTACCAGGGACTGCGCGCCGACCC‑GCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTcggcgg                      <  2:5979167/151‑1 (MQ=255)
                                                                                                                                cgccTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCggg                     <  2:10249397/148‑1 (MQ=255)
                                                                                                                                 gacTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCgggg                    >  2:8043622/1‑151 (MQ=255)
                                                                                                                                 gacTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTGCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCgggg                    <  1:3419673/151‑1 (MQ=255)
                                                                                                                                  acTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGa                   >  2:9587183/1‑151 (MQ=255)
                                                                                                                                   cTACCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAt                  <  2:3474060/151‑1 (MQ=255)
                                                                                                                                     aCCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAtct                <  2:210033/151‑1 (MQ=255)
                                                                                                                                     aCCAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGAtct                >  2:261965/1‑151 (MQ=255)
                                                                                                                                       cAGGGACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCg              >  2:9160353/1‑151 (MQ=255)
                                                                                                                                         gggACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaa            <  1:6231356/151‑1 (MQ=255)
                                                                                                                                         gggACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGaa            >  2:11146400/1‑151 (MQ=255)
                                                                                                                                           gACTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAg          <  2:8665323/151‑1 (MQ=255)
                                                                                                                                             cTGCGCGCCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTc        <  2:1977581/151‑1 (MQ=255)
                                                                                                                                                 gcgcCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCt    <  2:2940887/151‑1 (MQ=255)
                                                                                                                                                   gcCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTgg  >  2:1513534/1‑151 (MQ=255)
                                                                                                                                                   gcCTA‑CCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTgg  >  2:1174629/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GAGCTGGGTCGACCTGCGCCCAGCAGGCGCTGGTCTACCCTCGGCCGCTGGCCGGCGAGCCCGTTCCTCGCGCCGGGCGCGGAGGATGCCGAGGACCAGGGGCGCGCGCTGGCGAGCACGGGCGTCGACGACTACCAGGGACTGCGCGCCTACCCAGCCGGGCGACTCGCGGCGGCGCCTGCACTGGAAGGCGTTTTCCCGTGGCCAGGGGCTGCTGGTGAAGGACTTCGCCGCGCTCGGCGGGGAAGACCTGTGGCTGGAATTCGAGCGCCTCGGCGGGGATCTCGAAAGTCGCCTGG  >  NC_008463/2341550‑2341848

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