Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,895 Δ1 bp pseudogene (973/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,8930C.99.1% 1579.1 / 4.6 328pseudogene (971/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (166/157);  total (166/162)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

GCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCAGCGCCAACC  >  NC_008463/2787744‑2788039
                                                                                                                                                       |                                                                                                                                                  
    aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                                                                                                                                                <  1:2912357/151‑4 (MQ=255)
      cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtg                                                                                                                                             <  2:9843695/151‑1 (MQ=255)
      cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtg                                                                                                                                             <  2:3227269/151‑1 (MQ=255)
      cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtg                                                                                                                                             <  2:7280408/151‑1 (MQ=255)
       ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtgg                                                                                                                                            <  1:9885746/151‑1 (MQ=255)
        gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGa                                                                                                                                           >  1:6445310/1‑151 (MQ=255)
        gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGa                                                                                                                                           >  1:1674280/1‑151 (MQ=255)
         cGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGaa                                                                                                                                          <  2:1909622/151‑1 (MQ=255)
          gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                                                                                         <  2:7758201/151‑1 (MQ=255)
          gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAACTGGTGGAAg                                                                                                                                         >  2:1047396/1‑151 (MQ=255)
            gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggc                                                                                                                                       >  1:5148019/1‑151 (MQ=255)
             aaCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcg                                                                                                                                      <  1:616376/151‑1 (MQ=255)
              aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgc                                                                                                                                     >  2:5714930/1‑151 (MQ=255)
              aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgc                                                                                                                                     >  1:7263896/1‑151 (MQ=255)
                 gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggc                                                                                                                                  >  2:8858976/1‑151 (MQ=255)
                   cGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcgc                                                                                                                                >  1:1254461/1‑151 (MQ=255)
                    gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                                                                                               >  2:2215392/1‑151 (MQ=255)
                     ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTg                                                                                                                              >  2:2309023/1‑151 (MQ=255)
                      gAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTgg                                                                                                                             >  2:9676972/1‑151 (MQ=255)
                       aaGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGt                                                                                                                            <  1:5550969/151‑1 (MQ=255)
                        aGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTg                                                                                                                           >  2:7869741/1‑151 (MQ=255)
                           ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                                                                                        <  1:2726023/151‑1 (MQ=255)
                           ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                                                                                        >  1:5468598/1‑151 (MQ=255)
                           ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                                                                                        <  1:5002502/151‑1 (MQ=255)
                            gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                                                                                       <  1:10199469/151‑1 (MQ=255)
                            gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                                                                                       <  2:4466727/151‑1 (MQ=255)
                            gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                                                                                       <  2:7134621/151‑1 (MQ=255)
                            gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                                                                                       <  1:3229891/151‑1 (MQ=255)
                             cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgc                                                                                                                      <  1:8755340/151‑1 (MQ=255)
                             cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgc                                                                                                                      >  1:6032734/1‑151 (MQ=255)
                              ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                                                                                                     >  2:8386040/1‑151 (MQ=255)
                                    ggTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                                                                                               <  2:3654246/151‑1 (MQ=255)
                                    ggTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                                                                                               >  2:5279507/1‑151 (MQ=255)
                                     gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctct                                                                                                              >  1:4222237/1‑151 (MQ=255)
                                        ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCttac                                                                                                           >  2:10275384/2‑149 (MQ=255)
                                        ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                                                           >  2:2807953/2‑151 (MQ=255)
                                           cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                                                                                        <  2:6732076/151‑1 (MQ=255)
                                           cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                                                                                        >  1:4490187/1‑151 (MQ=255)
                                            aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                                                       >  2:3058034/1‑151 (MQ=255)
                                            aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCAGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                                                       >  2:10136882/1‑151 (MQ=255)
                                              cgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                                                                                     <  1:10685462/151‑1 (MQ=255)
                                                 cTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                                                                  >  1:147793/1‑151 (MQ=255)
                                                  tGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                                                                                 <  1:4840323/151‑1 (MQ=255)
                                                  tGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                                                                                 >  2:5475566/1‑151 (MQ=255)
                                                   ggaggaGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                                                                                                <  2:8594494/151‑1 (MQ=255)
                                                       gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCTACa                                                                                            >  2:7619811/1‑151 (MQ=255)
                                                        agaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACAt                                                                                           >  1:9225038/1‑151 (MQ=255)
                                                          aTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                                                                                         <  1:2319594/151‑1 (MQ=255)
                                                           tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCCCCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCGTCCTGTTCCTCGGCAACATCcg                                                                                        <  1:3737881/151‑1 (MQ=255)
                                                            cAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgc                                                                                       <  1:6199676/151‑1 (MQ=255)
                                                              aCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                                                                                     >  1:5541746/1‑151 (MQ=255)
                                                                cGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg                                                                                   <  2:6214833/151‑1 (MQ=255)
                                                                cGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg                                                                                   <  2:10011209/151‑1 (MQ=255)
                                                                  gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggcg                                                                                 <  2:3371330/151‑1 (MQ=255)
                                                                  gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggcg                                                                                 >  1:5224062/1‑151 (MQ=255)
                                                                  gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggcg                                                                                 <  1:194617/151‑1 (MQ=255)
                                                                         cccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa                                                                          >  1:7713372/5‑151 (MQ=255)
                                                                            aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                                                       >  1:2019661/2‑151 (MQ=255)
                                                                                   tcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGGTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcat                                                                <  1:531739/146‑1 (MQ=255)
                                                                                     aaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatca                                                              >  1:67327/4‑151 (MQ=255)
                                                                                        gCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATcc                                                           >  2:51213/1‑151 (MQ=255)
                                                                                         cGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATccc                                                          <  2:2084671/151‑1 (MQ=255)
                                                                                            tCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCt                                                       <  2:3268652/151‑1 (MQ=255)
                                                                                            tCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCt                                                       <  2:927963/151‑1 (MQ=255)
                                                                                              aCCGTCTACGACCGCACCACCCTTGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGt                                                     <  1:8150782/151‑1 (MQ=255)
                                                                                              aCCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGt                                                     <  2:3835340/151‑1 (MQ=255)
                                                                                               ccGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                                                    >  1:2546022/1‑151 (MQ=255)
                                                                                               ccGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                                                    >  1:12985/1‑151 (MQ=255)
                                                                                               ccGTCAACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                                                    <  1:566468/151‑1 (MQ=255)
                                                                                                cGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCg                                                   >  2:8231243/1‑151 (MQ=255)
                                                                                                 gTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGa                                                  >  1:8587934/1‑151 (MQ=255)
                                                                                                 gTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGGCGa                                                  >  2:5515403/1‑151 (MQ=255)
                                                                                                  tCTACGACCGCACCCCCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGACGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGAt                                                 <  1:1373992/151‑1 (MQ=255)
                                                                                                   cTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATg                                                >  2:4873649/1‑151 (MQ=255)
                                                                                                     aCGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCt                                              <  2:1602485/151‑1 (MQ=255)
                                                                                                      cGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg                                             <  2:4173221/151‑1 (MQ=255)
                                                                                                       gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGACGTAGATCCTGGTGGAAGGCGCGGCGCTGGTTATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGATGt                                            <  2:6934243/151‑1 (MQ=255)
                                                                                                       gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                                            >  1:7469807/1‑151 (MQ=255)
                                                                                                       gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                                            >  1:4686375/1‑151 (MQ=255)
                                                                                                        aCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGtt                                           >  1:9673448/1‑151 (MQ=255)
                                                                                                         ccGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttc                                          <  2:1113512/151‑1 (MQ=255)
                                                                                                         ccGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttc                                          <  1:6576876/151‑1 (MQ=255)
                                                                                                           gCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttcac                                        >  1:8318842/1‑151 (MQ=255)
                                                                                                           gCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGACACCATCATCCCGCTGTCGATGCTGttcac                                        >  1:5876919/1‑151 (MQ=255)
                                                                                                             accaccCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttcacct                                      >  1:3291865/1‑151 (MQ=255)
                                                                                                                            caaggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGGTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTc                        <  2:1146088/144‑1 (MQ=255)
                                                                                                                             aaggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCg                       >  1:7903733/7‑151 (MQ=255)
                                                                                                                             aaggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCg                       >  1:10428661/7‑151 (MQ=255)
                                                                                                                               ggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGc                     <  1:2442464/147‑1 (MQ=255)
                                                                                                                                  cgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCTCCTTCACCGGGATGGTCGGCAAc                  >  1:10537827/2‑151 (MQ=255)
                                                                                                                                  cgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAAc                  <  2:2584667/150‑1 (MQ=255)
                                                                                                                                   gTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGGTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACa                 <  2:5169522/151‑1 (MQ=255)
                                                                                                                                        cacggTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTc            >  1:2092330/5‑151 (MQ=255)
                                                                                                                                         acggTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAggggct           >  1:8946744/4‑148 (MQ=255)
                                                                                                                                         acggTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCa           >  1:7016582/4‑151 (MQ=255)
                                                                                                                                         acggTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCa           >  1:4042472/4‑151 (MQ=255)
                                                                                                                                            gTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCAgcg        <  2:6032734/151‑1 (MQ=255)
                                                                                                                                               aagaagaaCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCAGCGCCa     <  2:4650114/151‑1 (MQ=255)
                                                                                                                                               aagaagaaCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCAGCGCCa     <  1:10009046/151‑1 (MQ=255)
                                                                                                                                                  aagaaCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCAGCGCCAAcc  <  2:9543514/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
GCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAACAGGGTCAGCGCCAACC  >  NC_008463/2787744‑2788039

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: