Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,829 2 bp→AG pseudogene (907‑908/3090 nt) PA14_RS13085 → CusA/CzcA family heavy metal efflux RND transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,787,8290CA99.4% 1355.3 / NA 323pseudogene (907/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (142/178);  total (144/179)
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
*NC_0084632,787,8300TG99.7% 1254.4 / NA 323pseudogene (908/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (143/179);  minor base C (1/0);  total (144/179)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 323 total aligned reads displayed.

TGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCC‑AGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCAC  >  NC_008463/2787686‑2787974
                                                                                                                                                ||                                                                                                                                                 
tGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGa                                                                                                                                               <  2:3523054/150‑1 (MQ=255)
   gcACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCa                                                                                                                                            >  1:7951737/1‑150 (MQ=255)
    cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGACGTGGCCCAGGCGGTCGGCCAGCGCCGGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATccc                                                                                                                                           >  2:1716168/1‑148 (MQ=255)
      ccGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCTAGGCGATCACCg                                                                                                                                         >  1:11184697/1‑150 (MQ=255)
       cGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGt                                                                                                                                        <  2:7349523/150‑1 (MQ=255)
           gcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                                                                    <  2:6118505/150‑1 (MQ=255)
                cACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCg                                                                                                                               <  1:3256415/150‑1 (MQ=255)
                  ccGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGca                                                                                                                             >  2:948576/1‑150 (MQ=255)
                           gCCGCGAAGTGGTGCTCGGCACGGTGTACATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTgg                                                                                                                    <  2:5389626/150‑2 (MQ=255)
                             cgcgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtg                                                                                                                  >  1:8928244/1‑147 (MQ=255)
                              gcgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtgg                                                                                                                 >  2:1413854/1‑146 (MQ=255)
                                gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggac                                                                                                               >  2:3683019/1‑144 (MQ=255)
                                 aaGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggaca                                                                                                              >  2:7940786/1‑143 (MQ=255)
                                                        aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaag                                                                                      >  1:7693606/1‑150 (MQ=255)
                                                        aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaag                                                                                      <  2:9100128/150‑1 (MQ=255)
                                                          gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaagaa                                                                                    >  1:8177761/1‑150 (MQ=255)
                                                          gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaagaa                                                                                    >  2:10221531/1‑150 (MQ=255)
                                                           cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAc                                                                                   <  1:7541122/150‑1 (MQ=255)
                                                           cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAc                                                                                   <  1:3112603/150‑1 (MQ=255)
                                                             gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcct                                                                                 <  2:3854734/150‑2 (MQ=255)
                                                               tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctgg                                                                               <  1:8808582/150‑4 (MQ=255)
                                                                 ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtg                                                                            <  1:6187359/150‑1 (MQ=255)
                                                                  gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                           >  2:10476402/1‑150 (MQ=255)
                                                                  gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                           >  2:4101554/1‑150 (MQ=255)
                                                                  gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                           >  2:8931593/1‑150 (MQ=255)
                                                                    gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGaa                                                                         >  1:9786135/1‑150 (MQ=255)
                                                                    gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGaa                                                                         <  1:9246322/150‑1 (MQ=255)
                                                                     agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCNACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                        >  1:5866174/1‑150 (MQ=255)
                                                                      gtaCAGCCGGGAAGTGGCCCAGGCTGTCGGACAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAgg                                                                       <  1:489186/148‑1 (MQ=255)
                                                                      gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAgg                                                                       >  2:10081302/1‑150 (MQ=255)
                                                                      gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAgg                                                                       >  2:743183/1‑150 (MQ=255)
                                                                       aaCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggc                                                                      >  1:10677551/1‑150 (MQ=255)
                                                                        aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCCCG‑GTGAAGAAGAA‑CCTGGTGGAAggcg                                                                     >  1:6663886/1‑150 (MQ=255)
                                                                         cAGCCGGGAAGTGGCCCAGGCGGTCGGCCCGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgc                                                                    >  2:1143355/1‑150 (MQ=255)
                                                                          cgCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcg                                                                   <  2:10699883/149‑1 (MQ=255)
                                                                            ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcggc                                                                 >  2:8345008/1‑150 (MQ=255)
                                                                              gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcggcgc                                                               <  1:8492740/150‑1 (MQ=255)
                                                                                    tGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGa                                                         <  1:9547446/150‑1 (MQ=255)
                                                                                     ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGAt                                                        <  1:1108757/150‑1 (MQ=255)
                                                                                       cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                      <  1:5553090/150‑1 (MQ=255)
                                                                                       cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                      <  2:4462906/150‑1 (MQ=255)
                                                                                       cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                      <  2:4045708/150‑1 (MQ=255)
                                                                                          aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg                                                   <  1:205055/148‑1 (MQ=255)
                                                                                          aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg                                                   >  2:7674464/3‑150 (MQ=255)
                                                                                          aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg                                                   <  2:5991987/148‑1 (MQ=255)
                                                                                            gcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                 <  1:5020057/150‑1 (MQ=255)
                                                                                               gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                              >  1:7841201/1‑150 (MQ=255)
                                                                                               gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                              <  1:10304013/150‑1 (MQ=255)
                                                                                                tcggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctct                                             <  2:4970100/147‑1 (MQ=255)
                                                                                                tcggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctct                                             <  1:10806287/147‑1 (MQ=255)
                                                                                                 cggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCtt                                            >  1:9104765/3‑150 (MQ=255)
                                                                                                  ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTc                                           <  2:7012652/149‑1 (MQ=255)
                                                                                                   gCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGACAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGCAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                          <  2:9036963/150‑1 (MQ=255)
                                                                                                    ccAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                         <  1:5758551/150‑1 (MQ=255)
                                                                                                    ccAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGGGGACAAGG‑CCGTGGCCACGGAGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                         <  2:4569194/150‑1 (MQ=255)
                                                                                                      aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                       <  1:1409323/150‑1 (MQ=255)
                                                                                                      aGCGCCTGGAGGAGACCAACCGGCCCCGGCCCAAGGGGGTCAAGCCGATCACCGTCTACGACCGCACCACCCTGGTGCACAAGG‑CCGTGGCCACGGTGAAGCAGAACCC‑GGTGGAAGGCGCGGCGCTGGAGATCGCGGTGCTCTTCCTGt                                       <  2:922650/150‑1 (MQ=255)
                                                                                                       gcgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                      >  2:10578878/1‑150 (MQ=255)
                                                                                                       gcgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                      <  1:4375444/150‑1 (MQ=255)
                                                                                                         gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                    >  1:1000659/1‑150 (MQ=255)
                                                                                                         gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                    >  1:1938083/1‑150 (MQ=255)
                                                                                                         gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                    <  1:3458916/150‑1 (MQ=255)
                                                                                                           cTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTc                                  <  2:2881723/150‑1 (MQ=255)
                                                                                                            tGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                 <  1:3162/150‑1 (MQ=255)
                                                                                                            tGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTG‑GTGGCCAAGGCCGTGGCCACGGTGACGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                 >  2:2211963/1‑150 (MQ=255)
                                                                                                             ggaggaGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                <  2:3142769/150‑1 (MQ=255)
                                                                                                             ggaggaGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                <  2:3328062/150‑1 (MQ=255)
                                                                                                             ggaggaGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                <  1:7296969/150‑1 (MQ=255)
                                                                                                             ggaggaGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCGGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                <  2:5785577/150‑1 (MQ=255)
                                                                                                              gaggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                               <  2:2578230/150‑1 (MQ=255)
                                                                                                              gaggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                               <  1:6943023/150‑1 (MQ=255)
                                                                                                                ggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCaa                             <  1:4019808/150‑1 (MQ=255)
                                                                                                                ggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCaa                             <  1:4558202/150‑1 (MQ=255)
                                                                                                                 gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                            <  1:2952109/150‑1 (MQ=255)
                                                                                                                    aTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATc                         >  1:2198927/1‑150 (MQ=255)
                                                                                                                    aTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATc                         <  1:2617/150‑1 (MQ=255)
                                                                                                                     tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                        >  2:9819988/1‑150 (MQ=255)
                                                                                                                     tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                        >  2:4215132/1‑150 (MQ=255)
                                                                                                                     tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                        >  1:1383744/1‑150 (MQ=255)
                                                                                                                          cGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                   <  1:1750169/150‑1 (MQ=255)
                                                                                                                            gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                 >  2:6513870/1‑150 (MQ=255)
                                                                                                                              cccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGCAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGgcgc               <  1:5252367/150‑1 (MQ=255)
                                                                                                                              cccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGgcgc               >  2:1966740/1‑150 (MQ=255)
                                                                                                                               ccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCt              <  2:758181/150‑1 (MQ=255)
                                                                                                                                 tGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGa            <  2:9600673/150‑1 (MQ=255)
                                                                                                                                 tGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGa            <  1:4152441/150‑1 (MQ=255)
                                                                                                                                 tGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGa            >  2:8548027/1‑150 (MQ=255)
                                                                                                                                  gcccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGAt           <  1:5924740/145‑1 (MQ=255)
                                                                                                                                   cccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATc          >  2:5025682/5‑150 (MQ=255)
                                                                                                                                    ccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa         >  1:8436691/4‑150 (MQ=255)
                                                                                                                                     caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc        >  1:10003595/3‑150 (MQ=255)
                                                                                                                                     caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc        >  2:1633741/3‑150 (MQ=255)
                                                                                                                                      aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc       >  1:8128756/2‑150 (MQ=255)
                                                                                                                                       aGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg      <  2:10994810/150‑1 (MQ=255)
                                                                                                                                       aGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg      >  2:5974140/1‑150 (MQ=255)
                                                                                                                                        gggggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGc     >  1:633683/1‑150 (MQ=255)
                                                                                                                                        gggggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGc     >  1:6493140/1‑150 (MQ=255)
                                                                                                                                         ggggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGTCCACGGTGAAGAAGCA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGcc    <  2:8888979/150‑1 (MQ=255)
                                                                                                                                          gggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGcca   <  1:9940944/150‑1 (MQ=255)
                                                                                                                                           ggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccac  <  2:6595473/150‑1 (MQ=255)
                                                                                                                                                ||                                                                                                                                                 
TGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCC‑AGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCAC  >  NC_008463/2787686‑2787974

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