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.8% 2730.2 / ‑2.3 550pseudogene (158/718 nt)PA14_RS15905DUF1294 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (199/350);  total (200/350)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 589 total aligned reads displayed.

CCGGCCGTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTACCCGTT  >  NC_008463/3489825‑3490121
                                                                                                                                                      |                                                                                                                                                  
    ccgTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCgggga                                                                                                                                                <  1:7490543/151‑2 (MQ=255)
     cgTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCggggac                                                                                                                                               >  2:8856312/1‑149 (MQ=255)
     cgTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCggggac                                                                                                                                               >  1:4680764/1‑149 (MQ=255)
      gTGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCggggacg                                                                                                                                              <  2:8228170/151‑4 (MQ=255)
       tGCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCggggacgg                                                                                                                                             <  1:2395794/151‑5 (MQ=255)
        gCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGAcggc                                                                                                                                           >  1:5336834/1‑151 (MQ=255)
        gCCTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGAcggc                                                                                                                                           <  1:850301/151‑1 (MQ=255)
          cTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGAcggcgg                                                                                                                                         >  1:455052/1‑151 (MQ=255)
          cTGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGAcggcgg                                                                                                                                         >  1:6019119/1‑151 (MQ=255)
           tGGGCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGt                                                                                                                                        <  2:7433708/151‑1 (MQ=255)
              gCCCGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCAc                                                                                                                                     >  2:5773556/1‑151 (MQ=255)
               cccGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCAcc                                                                                                                                    >  1:3684897/1‑151 (MQ=255)
               cccGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCAcc                                                                                                                                    <  2:4389352/151‑1 (MQ=255)
                ccGCGTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACcg                                                                                                                                   <  1:8432509/151‑1 (MQ=255)
                  gcgTCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACcgcg                                                                                                                                 <  2:3528445/151‑1 (MQ=255)
                     tCACCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGCGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGaa                                                                                                                              >  2:1046203/1‑151 (MQ=255)
                       aCCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCGGGCGGGGACGGCGGTCACCGCGGAAcg                                                                                                                            <  2:2209177/151‑1 (MQ=255)
                       aCCGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGAGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGTCGGCGGTCACCGCGGAAcg                                                                                                                            <  2:8031845/151‑1 (MQ=255)
                        ccGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAAcgc                                                                                                                           >  2:4764478/1‑151 (MQ=255)
                         cGGCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCGCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGTAcgcg                                                                                                                          >  2:6607326/1‑151 (MQ=255)
                          ggCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCgg                                                                                                                         <  1:784844/151‑1 (MQ=255)
                          ggCGCCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCGCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCgg                                                                                                                         >  2:7738616/1‑151 (MQ=255)
                             gcCTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGgaga                                                                                                                      >  2:8734215/1‑151 (MQ=255)
                               cTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGata                                                                                                                    <  1:8914634/151‑1 (MQ=255)
                               cTGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGata                                                                                                                    <  2:486365/151‑1 (MQ=255)
                                tGGGCGCCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGatat                                                                                                                   >  1:6873479/1‑151 (MQ=255)
                                     gcCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAgg                                                                                                              <  2:4235215/151‑1 (MQ=255)
                                     gcCACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAgg                                                                                                              <  2:9871000/151‑1 (MQ=255)
                                       cACGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGaa                                                                                                            <  1:10695911/151‑1 (MQ=255)
                                        aCGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAg                                                                                                           <  1:8297081/151‑1 (MQ=255)
                                        aCGGCGGACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAg                                                                                                           <  1:7970518/151‑1 (MQ=255)
                                          ggcggACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGAc                                                                                                         <  1:2178584/151‑1 (MQ=255)
                                           gcggACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGAct                                                                                                        >  1:2016436/1‑151 (MQ=255)
                                             ggACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtc                                                                                                      >  1:2786751/1‑151 (MQ=255)
                                              gccAGTCGGCGGGTCGAGGGGCAGCGCGGCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtct                                                                                                     <  2:6134330/149‑1 (MQ=255)
                                              gACAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtct                                                                                                     >  2:9517583/1‑151 (MQ=255)
                                               aCAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtctt                                                                                                    >  1:143159/1‑151 (MQ=255)
                                               aCAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGActtctt                                                                                                    <  1:6661893/151‑1 (MQ=255)
                                                cAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACttcttc                                                                                                   <  1:2612022/151‑1 (MQ=255)
                                                cAGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACttcttc                                                                                                   <  2:6330847/151‑1 (MQ=255)
                                                 cgccggcggGCCGCGGGGCAGCGCGTCCAGCCGCGCCCCCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCCGGCGAGGAACCACCCCCGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGCGCTCTGCAGGAACACTTCTTcc                                                                                                  <  2:11022133/148‑1 (MQ=255)
                                                 aGTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTcc                                                                                                  <  1:4332309/151‑1 (MQ=255)
                                                  gTCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTccc                                                                                                 <  1:1186080/151‑1 (MQ=255)
                                                   tCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTcccc                                                                                                <  1:1380343/151‑1 (MQ=255)
                                                   tCGGCGGGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTcccc                                                                                                <  1:5594898/151‑1 (MQ=255)
                                                     ggcggGTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCccgc                                                                                              >  2:7527359/1‑151 (MQ=255)
                                                          gTCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTc                                                                                         <  1:3134816/151‑1 (MQ=255)
                                                           tCGAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCt                                                                                        <  2:1272445/151‑1 (MQ=255)
                                                             gAGGGTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTc                                                                                      >  1:6133253/1‑151 (MQ=255)
                                                               gggTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCgg                                                                                    <  2:1834906/151‑1 (MQ=255)
                                                               gggTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCgg                                                                                    <  2:1513741/151‑1 (MQ=255)
                                                                ggTCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGc                                                                                   <  2:2102762/151‑1 (MQ=255)
                                                                ggTCAGCGCGTCCAGCCGCGCCCTCTCGGAGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGCCGGCGGTCACCGCGGAACGCGGAGATATGCAGGACGACTTCTTCCCCGCCGGTCTTCGGc                                                                                   >  1:3294853/1‑151 (MQ=255)
                                                                  tCAGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATGTGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAg                                                                                 >  1:2501165/1‑151 (MQ=255)
                                                                    aGCGCGTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAggt                                                                               >  1:501611/1‑150 (MQ=255)
                                                                      cgcgTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATg                                                                             >  1:523930/1‑151 (MQ=255)
                                                                       gcgTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATga                                                                            <  2:5501966/151‑1 (MQ=255)
                                                                       gcgTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATga                                                                            <  1:10213509/151‑1 (MQ=255)
                                                                         gTCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGACCGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaag                                                                          >  1:8119129/1‑151 (MQ=255)
                                                                          tCCAGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagc                                                                         <  1:3221575/151‑1 (MQ=255)
                                                                             aGCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagccga                                                                      >  1:8760780/1‑151 (MQ=255)
                                                                              gCCGCGCCCTCTCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATgaagccgaa                                                                     <  2:9068447/151‑1 (MQ=255)
                                                                                      ctctCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTgt                                                             <  1:7623065/151‑1 (MQ=255)
                                                                                      ctctCGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTgt                                                             >  2:340780/1‑151 (MQ=255)
                                                                                          cGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTgtcgtc                                                         >  2:2661469/1‑151 (MQ=255)
                                                                                          cGGCGCGCAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTgtcgtc                                                         <  1:7078256/151‑1 (MQ=255)
                                                                                            gcgcgcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTcgtcgt                                                       >  1:699200/1‑151 (MQ=255)
                                                                                             cgcgcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTa                                                      <  1:10306627/151‑1 (MQ=255)
                                                                                               cgcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTAcc                                                    >  1:6382905/1‑151 (MQ=255)
                                                                                                gcAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCa                                                   >  1:11187524/1‑151 (MQ=255)
                                                                                                 cAGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAg                                                  <  2:9025067/151‑1 (MQ=255)
                                                                                                  aGGCGGCCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGc                                                 >  1:9949199/1‑151 (MQ=255)
                                                                                                     cggcCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGc                                              <  2:2376441/151‑1 (MQ=255)
                                                                                                      ggcCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCt                                             >  1:10777405/1‑151 (MQ=255)
                                                                                                      ggcCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCt                                             <  2:3063659/151‑1 (MQ=255)
                                                                                                       gcCCTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTg                                            <  1:3418289/151‑1 (MQ=255)
                                                                                                        cccTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGa                                           >  1:7253879/1‑151 (MQ=255)
                                                                                                        cccTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGa                                           <  1:5780973/151‑1 (MQ=255)
                                                                                                          cTGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATc                                         <  1:3781259/151‑1 (MQ=255)
                                                                                                           tGGGCGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCa                                        >  2:3907344/1‑151 (MQ=255)
                                                                                                               cGTCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGcc                                    <  1:2931061/151‑1 (MQ=255)
                                                                                                                 tCCTGGCTGGCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCg                                  <  2:1652752/151‑1 (MQ=255)
                                                                                                                   ctggctggCGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGc                                <  1:1751474/151‑1 (MQ=255)
                                                                                                                    tggctggcGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCt                               >  1:10026928/1‑151 (MQ=255)
                                                                                                                         ggcGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTcc                          <  1:9597129/151‑1 (MQ=255)
                                                                                                                          gcGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCa                         <  1:9246725/151‑1 (MQ=255)
                                                                                                                          gcGAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCa                         <  2:6905488/151‑1 (MQ=255)
                                                                                                                            gAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCata                       >  1:3926751/1‑151 (MQ=255)
                                                                                                                            gAGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCata                       <  2:10129383/151‑1 (MQ=255)
                                                                                                                            caGGAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGAGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCAGGCCCGGCTTTTCCatc                       >  2:11193659/2‑150 (MQ=255)
                                                                                                                               gAACCACACCTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATcc                    >  2:10836494/1‑151 (MQ=255)
                                                                                                                                       ccTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCActcgctc            <  2:5131524/151‑1 (MQ=255)
                                                                                                                                       ccTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCActcgctc            <  2:7252230/151‑1 (MQ=255)
                                                                                                                                       ccTGGTCGCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCActcgctc            <  2:386779/151‑1 (MQ=255)
                                                                                                                                         tGGTCGCCCTGGCGGGGACGGCGGTCCCCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCAGATCCACtcgctcgc          >  2:3161836/1‑151 (MQ=255)
                                                                                                                                              gCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTAccc     >  1:4089731/1‑151 (MQ=255)
                                                                                                                                              gCCCTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGAGGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCGTGCCCGGCTTTTCCATCTCCACTCGCTCGCTAccc     >  1:4723234/1‑151 (MQ=255)
                                                                                                                                                ccTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATACCCACTCGCTCGCTACCCGt   <  2:10733354/151‑1 (MQ=255)
                                                                                                                                                 cTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCTTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTACCCGtt  >  2:4961295/1‑151 (MQ=255)
                                                                                                                                                 cTGGCGGGGACGGCGGTCACCGCGGAACGCGGAGATATGCAGGAAGACTTCTTCCCCGCCGGTCTTCGGCAGGATGAAGCCGAAGCCCGTGTCGTCGTACCAGCGGCTGATCATGCCCGGCTTTTCCATATCCACTCGCTCGCTACCCGtt  >  1:5121670/1‑151 (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: