Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,489,979 (G)5→4 pseudogene (154/718 nt) PA14_RS15905 ← DUF1294 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,489,9750G.99.7% 3004.2 / 0.7 607pseudogene (158/718 nt)PA14_RS15905DUF1294 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (231/374);  total (233/374)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 639 total aligned reads displayed.

CCGGCCGTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTACCCGTT  >  NC_008463/3489825‑3490121
                                                                                                                                                      |                                                                                                                                                  
ccggccgTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCg                                                                                                                                                    >  2:10185652/1‑151 (MQ=255)
ccggccgTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCg                                                                                                                                                    <  2:4443323/151‑1 (MQ=255)
    ccgTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCgggga                                                                                                                                                <  1:10411288/151‑2 (MQ=255)
     cgTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCggggac                                                                                                                                               >  2:5060347/1‑149 (MQ=255)
         ccTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGAcggcg                                                                                                                                          <  2:753769/151‑1 (MQ=255)
           tGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGt                                                                                                                                        <  1:2574151/151‑1 (MQ=255)
             ggCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCa                                                                                                                                      >  1:1854460/1‑151 (MQ=255)
             ggCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCa                                                                                                                                      <  2:379764/151‑1 (MQ=255)
              gCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCAc                                                                                                                                     <  1:9880576/151‑1 (MQ=255)
                 cgcgTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGACCCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACcgc                                                                                                                                  >  1:9098831/1‑151 (MQ=255)
                  gcgTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACcgcg                                                                                                                                 <  1:1192314/151‑1 (MQ=255)
                  gcgTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACcgcg                                                                                                                                 <  1:2821243/151‑1 (MQ=255)
                   cgTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCgg                                                                                                                                <  1:4961322/151‑1 (MQ=255)
                    gTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGa                                                                                                                               >  2:4258273/1‑151 (MQ=255)
                    gTCACCGGCGCCTGGGCGCCACGGCGGACAGACGGCGGTTAGAGGGTCAGCGCGTCCAGCCGCGCCCACACGTCGCGCAGGCGGCCCAGGGCGTCCTGGCAGGCGATGAACCACACCTGGACGCCCTGGC‑GCGGACGCCGCTCACagcgcc                                                                                                                               >  1:5524756/1‑149 (MQ=255)
                      cACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAAc                                                                                                                             <  1:9197375/151‑1 (MQ=255)
                      cACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAAc                                                                                                                             <  1:4423874/151‑1 (MQ=255)
                       aCCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAAcg                                                                                                                            >  2:10444052/1‑151 (MQ=255)
                         cGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAAcgcg                                                                                                                          >  2:5386488/1‑151 (MQ=255)
                         cGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAAcgcg                                                                                                                          <  1:2153443/151‑1 (MQ=255)
                          ggCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCgg                                                                                                                         <  1:6780042/151‑1 (MQ=255)
                                tGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGatat                                                                                                                   <  2:2188879/151‑1 (MQ=255)
                                tGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGatat                                                                                                                   <  2:3425391/151‑1 (MQ=255)
                                  ggCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGc                                                                                                                 >  1:3565436/1‑151 (MQ=255)
                                        aCGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAg                                                                                                           <  1:8460848/151‑1 (MQ=255)
                                         cggcggACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGa                                                                                                          <  2:10575199/151‑1 (MQ=255)
                                             ggACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtc                                                                                                      <  2:582484/151‑1 (MQ=255)
                                             ggACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtc                                                                                                      <  1:5770333/151‑1 (MQ=255)
                                              gcccggcggcggGTCGAGGGTCAGCGCGTCCAGCCGCGCCCACTCGGCGCGCAGGCGGCCCGGGGCGTCCTGGCTGGCGCGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGACACCGCGGAACGCGGAGAGATGCAGGAAGActtct                                                                                                     <  2:10028611/145‑1 (MQ=255)
                                              gACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtct                                                                                                     <  2:1985032/151‑1 (MQ=255)
                                               aCAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtctt                                                                                                    <  1:8427937/151‑1 (MQ=255)
                                                cAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACttcttc                                                                                                   >  2:2474558/1‑151 (MQ=255)
                                                 aGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTcc                                                                                                  >  2:9307267/1‑151 (MQ=255)
                                                 aGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTcc                                                                                                  <  1:6558929/151‑1 (MQ=255)
                                                 aGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTcc                                                                                                  <  2:10879946/151‑1 (MQ=255)
                                                  gTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTccc                                                                                                 <  2:5508113/151‑1 (MQ=255)
                                                  gTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTccc                                                                                                 <  2:3628804/151‑1 (MQ=255)
                                                     ggcggGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCccgc                                                                                              >  1:3760092/1‑151 (MQ=255)
                                                      gcggGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCccgcc                                                                                             <  2:8322197/151‑1 (MQ=255)
                                                         ggTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGt                                                                                          >  1:3099961/1‑151 (MQ=255)
                                                          gTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTc                                                                                         <  1:6894249/151‑1 (MQ=255)
                                                          gTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATAGGCAGGAAGACGTCTTCCCCGCCGGTc                                                                                         <  2:8515413/151‑1 (MQ=255)
                                                           tCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCt                                                                                        <  1:10864958/151‑1 (MQ=255)
                                                             gAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTc                                                                                      <  1:1322231/151‑1 (MQ=255)
                                                             gAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTc                                                                                      <  2:7654329/151‑1 (MQ=255)
                                                              aGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCGTCCCCGCCGGTCTTCg                                                                                     <  2:10587010/151‑1 (MQ=255)
                                                               gggTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑TGGGACGGCGGTCACCGCGGACCGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCgg                                                                                    >  1:2427007/1‑151 (MQ=255)
                                                               gggTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCgg                                                                                    <  2:7252836/151‑1 (MQ=255)
                                                                   cAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCCGGCGGCCCTGGGCGTCCTGGCTGGCGAGGACCCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAgg                                                                                >  2:10150956/1‑151 (MQ=255)
                                                                   cAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAgg                                                                                <  2:1806706/151‑1 (MQ=255)
                                                                     gcgcGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGAt                                                                              <  2:4906262/151‑1 (MQ=255)
                                                                          tCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagc                                                                         <  1:206017/151‑1 (MQ=255)
                                                                            cAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagccg                                                                       <  2:2358783/151‑1 (MQ=255)
                                                                            cAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagccg                                                                       >  1:324647/1‑151 (MQ=255)
                                                                            cAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagccg                                                                       <  2:9518834/151‑1 (MQ=255)
                                                                             aGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTAACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagccga                                                                      <  1:10721818/151‑1 (MQ=255)
                                                                                cgcgCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagccgaagc                                                                   <  1:9445449/151‑1 (MQ=255)
                                                                                  cgcCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGccc                                                                 >  1:9385538/1‑151 (MQ=255)
                                                                                   gcCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCt                                                                <  2:7992815/151‑1 (MQ=255)
                                                                                     ccTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTg                                                              <  1:7214029/151‑1 (MQ=255)
                                                                                      ctctCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTgt                                                             <  1:10158734/151‑1 (MQ=255)
                                                                                        ctcGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTgtcg                                                           >  1:1761865/1‑151 (MQ=255)
                                                                                        ctcGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTGCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTgtcg                                                           >  1:8047942/1‑151 (MQ=255)
                                                                                             cgcgcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTa                                                      <  1:9654366/151‑1 (MQ=255)
                                                                                             cgcgcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTa                                                      >  1:8281766/1‑151 (MQ=255)
                                                                                             cgcgcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTa                                                      <  2:9235118/151‑1 (MQ=255)
                                                                                               cgcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTAcc                                                    <  1:6108504/151‑1 (MQ=255)
                                                                                                gcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCa                                                   <  1:2935988/151‑1 (MQ=255)
                                                                                                gcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCa                                                   <  1:3665297/151‑1 (MQ=255)
                                                                                                  aGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGc                                                 <  2:2808313/151‑1 (MQ=255)
                                                                                                  aGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGc                                                 >  2:10739917/1‑151 (MQ=255)
                                                                                                      ggcCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCt                                             <  2:6648057/151‑1 (MQ=255)
                                                                                                         ccTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGAt                                          >  1:871868/1‑151 (MQ=255)
                                                                                                          cTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATc                                         >  2:4668411/1‑151 (MQ=255)
                                                                                                           tGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCa                                        <  2:1707370/151‑1 (MQ=255)
                                                                                                           tGGGCGTCCGGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCa                                        <  1:5395968/151‑1 (MQ=255)
                                                                                                               cGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGcc                                    <  1:5956764/151‑1 (MQ=255)
                                                                                                                gTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGccc                                   >  1:7015644/1‑151 (MQ=255)
                                                                                                                gTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCGTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGccc                                   <  2:1303051/151‑1 (MQ=255)
                                                                                                                 tCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCTGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCg                                  <  1:395376/151‑1 (MQ=255)
                                                                                                                 tCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCg                                  <  1:10938273/151‑1 (MQ=255)
                                                                                                                 tCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCg                                  >  2:7585427/1‑151 (MQ=255)
                                                                                                                     ggctggcGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCtt                              >  1:4788474/1‑151 (MQ=255)
                                                                                                                      gctggcGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCttt                             >  1:7734014/1‑151 (MQ=255)
                                                                                                                          gcGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCa                         <  1:5587376/151‑1 (MQ=255)
                                                                                                                           cGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCat                        >  1:6402076/1‑151 (MQ=255)
                                                                                                                           cGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCat                        <  1:7468451/151‑1 (MQ=255)
                                                                                                                           cGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCat                        >  1:7868572/1‑151 (MQ=255)
                                                                                                                           cGAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCat                        <  2:1426131/151‑1 (MQ=255)
                                                                                                                            gAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCata                       <  1:5089925/151‑1 (MQ=255)
                                                                                                                            gAGGAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCata                       >  2:8822102/1‑151 (MQ=255)
                                                                                                                              ggAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATc                     >  2:6394122/1‑151 (MQ=255)
                                                                                                                               gAACCACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATcc                    <  2:6959049/151‑1 (MQ=255)
                                                                                                                                  ccACACCTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCAct                 <  2:3187797/151‑1 (MQ=255)
                                                                                                                                        cTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCActcgctcg           >  1:5164805/1‑151 (MQ=255)
                                                                                                                                        cTGGTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCActcgctcg           >  1:7652138/1‑151 (MQ=255)
                                                                                                                                           gTCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTa        <  2:6111539/151‑1 (MQ=255)
                                                                                                                                            tCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTAc       <  1:1874323/151‑1 (MQ=255)
                                                                                                                                            tCGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTAc       <  2:9738052/151‑1 (MQ=255)
                                                                                                                                             cGCCCTGGC‑GGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTAcc      <  2:5312559/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCGGCCGTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTACCCGTT  >  NC_008463/3489825‑3490121

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: ‑