Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,861 3 bp→GTG pseudogene (939‑941/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,8610CG99.4% 1291.7 / NA 338pseudogene (939/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base G (152/181);  minor base T (1/1);  total (154/184)
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,8620AT99.1% 1399.1 / ‑3.2 338pseudogene (940/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (152/183);  minor base G (2/1);  total (154/184)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084632,787,8630TG99.4% 1311.6 / ‑2.7 338pseudogene (941/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (153/183);  total (154/184)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

TCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCC‑AGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCT  >  NC_008463/2787740‑2788003
                                                                                                                          |||                                                                                                                                             
tCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaa                                                                                                                    >  1:8210883/1‑151 (MQ=255)
 cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaag                                                                                                                   <  1:1631331/151‑1 (MQ=255)
  aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaaga                                                                                                                  <  2:4358545/151‑1 (MQ=255)
   tGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaagaa                                                                                                                 <  1:2570417/151‑1 (MQ=255)
    gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAc                                                                                                                >  1:2092920/1‑151 (MQ=255)
    gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAc                                                                                                                <  1:2786110/151‑1 (MQ=255)
      tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcct                                                                                                              >  2:7689982/1‑150 (MQ=255)
       gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctg                                                                                                             <  1:684895/151‑3 (MQ=255)
       gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctg                                                                                                             >  2:6061516/1‑149 (MQ=255)
         tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctggt                                                                                                           <  2:2710776/151‑5 (MQ=255)
         tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctggt                                                                                                           <  1:8361682/151‑5 (MQ=255)
          cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGTA‑CCtggtg                                                                                                         <  1:10922642/151‑1 (MQ=255)
          cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtg                                                                                                         >  1:8689612/1‑151 (MQ=255)
           ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                                                        <  2:3127805/151‑1 (MQ=255)
           ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                                                        >  1:4559131/1‑151 (MQ=255)
             cGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGaa                                                                                                      >  2:1890061/1‑151 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCATGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGACGAAGAA‑CCTGGTGGacg                                                                                                     >  1:1116074/1‑149 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGCTCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                                                     <  2:5261776/151‑1 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                                                     <  2:5861574/151‑1 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                                                     >  2:2566392/1‑151 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                                                     <  1:6549408/151‑1 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                                                     <  2:4823881/151‑1 (MQ=255)
               agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAgg                                                                                                    <  1:8068840/151‑1 (MQ=255)
                gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggc                                                                                                   >  2:2811112/1‑151 (MQ=255)
                 aaCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcg                                                                                                  <  2:7411154/151‑1 (MQ=255)
                  aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgc                                                                                                 <  2:125481/151‑1 (MQ=255)
                    aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                                               <  2:927893/151‑1 (MQ=255)
                    aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                                               <  2:11018596/151‑1 (MQ=255)
                     gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcggc                                                                                              >  2:10769825/1‑151 (MQ=255)
                     gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcggc                                                                                              <  1:3090715/151‑1 (MQ=255)
                        gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCt                                                                                           <  2:9706027/151‑1 (MQ=255)
                          gAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTgg                                                                                         <  1:979528/151‑1 (MQ=255)
                           aaGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGt                                                                                        <  2:7233201/151‑1 (MQ=255)
                           aaGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGACCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCCCCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGCA‑CCTGGTGGAAGGCGCGGCGCTGGt                                                                                        <  2:3617983/151‑1 (MQ=255)
                             gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTG‑GTGGCCAAGGCCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGa                                                                                      >  2:6416765/1‑151 (MQ=255)
                              tGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGAt                                                                                     <  2:5755971/151‑1 (MQ=255)
                               ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATc                                                                                    >  1:9202615/1‑151 (MQ=255)
                                gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGAC‑ACGGCCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcg                                                                                   >  2:2645715/1‑151 (MQ=255)
                                    aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACCAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt                                                                               >  1:4009813/3‑151 (MQ=255)
                                     ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                                              <  2:4181898/150‑1 (MQ=255)
                                     ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                                              <  2:887067/150‑1 (MQ=255)
                                       cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                                                            >  2:8382593/1‑151 (MQ=255)
                                          tcggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCtt                                                                         <  1:4295899/148‑1 (MQ=255)
                                            ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                       >  2:3969723/2‑151 (MQ=255)
                                               cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                                                    <  1:4465710/151‑1 (MQ=255)
                                                aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                   >  2:231089/1‑151 (MQ=255)
                                                aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                   >  1:10518576/1‑151 (MQ=255)
                                                aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                   <  1:5826456/151‑1 (MQ=255)
                                                  cgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                                                 <  1:5393415/151‑1 (MQ=255)
                                                    ccTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTc                                                               <  1:7918606/151‑1 (MQ=255)
                                                    ccTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTc                                                               <  2:514061/151‑1 (MQ=255)
                                                     cTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                              <  1:9950723/151‑1 (MQ=255)
                                                     cTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                              >  1:5214175/1‑151 (MQ=255)
                                                         aggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCaa                                                          >  1:2850931/1‑151 (MQ=255)
                                                         aggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCaa                                                          <  1:1448091/151‑1 (MQ=255)
                                                           gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                                                        <  2:9559760/151‑1 (MQ=255)
                                                           gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                                                        <  1:5911943/151‑1 (MQ=255)
                                                            agaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACAt                                                       <  2:5357708/151‑1 (MQ=255)
                                                            agaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACAt                                                       <  2:1689735/151‑1 (MQ=255)
                                                             gaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATc                                                      >  1:5677830/1‑151 (MQ=255)
                                                              aTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                                                     <  1:3470136/151‑1 (MQ=255)
                                                                cAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgc                                                   >  2:8190918/1‑151 (MQ=255)
                                                                 aaCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                                                  <  2:4979418/151‑1 (MQ=255)
                                                                 aaCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGAGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                                                  <  1:6596004/151‑1 (MQ=255)
                                                                   ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                                                <  1:6069589/151‑1 (MQ=255)
                                                                   ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                                                >  2:7008702/1‑151 (MQ=255)
                                                                   ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                                                >  1:10705429/1‑151 (MQ=255)
                                                                   ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                                                <  1:7309180/151‑1 (MQ=255)
                                                                   ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCGGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                                                <  2:2368204/151‑1 (MQ=255)
                                                                     ggACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                              >  1:10011890/1‑151 (MQ=255)
                                                                      gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggcg                                             >  2:9388954/1‑151 (MQ=255)
                                                                            gcccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATc                                       <  2:8710535/146‑1 (MQ=255)
                                                                              ccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc                                     >  1:10197084/4‑151 (MQ=255)
                                                                               caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc                                    >  2:3270445/3‑151 (MQ=255)
                                                                               caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc                                    <  2:4657205/149‑1 (MQ=255)
                                                                                aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                   <  1:8985725/150‑1 (MQ=255)
                                                                                aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                   >  2:1402365/2‑151 (MQ=255)
                                                                                 aGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGc                                  >  1:406049/1‑151 (MQ=255)
                                                                                        caaggCGATCTCCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCAACGGTGTAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatc                           <  2:9029895/147‑1 (MQ=255)
                                                                                        caaggCGATCACCGTCTCCGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatc                           <  2:6477527/147‑1 (MQ=255)
                                                                                         aaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatca                          <  1:4732985/148‑1 (MQ=255)
                                                                                         aaggCGATCACCGTCTACGACAGCACAACCCTGTTGGACAAGG‑CCGTGGCCACGGAGAAGAATAA‑CAAGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTTTTCCTCGGCAACATCCGCTCGGCGCTGATCACCGCCACcatca                          >  1:3666472/4‑151 (MQ=255)
                                                                                          aggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatcat                         <  1:2371809/149‑1 (MQ=255)
                                                                                           ggCGATCACCGTCTACGACCGCACCACCCTTGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGAGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCatcatc                        <  2:8154970/150‑1 (MQ=255)
                                                                                                tCACCGTCTACGACCGCACCACCCTGGTGGAC‑ATGGCCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGTTCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCAGCATCCCGCt                   >  1:4620118/1‑151 (MQ=255)
                                                                                                   ccGTCTACGACCGCACCACCCTGGTGGACCAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGTTCGCGGTGCTCTTCCTGTTCCTCGGCAACAGCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                >  2:7906569/1‑151 (MQ=255)
                                                                                                   ccGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                <  2:10112422/151‑1 (MQ=255)
                                                                                                   ccGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                >  1:3160739/1‑151 (MQ=255)
                                                                                                   ccGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                <  1:8548528/151‑1 (MQ=255)
                                                                                                   ccGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                >  1:6339909/1‑151 (MQ=255)
                                                                                                    cGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCg               <  2:2149885/151‑1 (MQ=255)
                                                                                                     gTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGa              >  1:10098230/1‑151 (MQ=255)
                                                                                                     gTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGa              >  1:7168258/1‑151 (MQ=255)
                                                                                                     gTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGa              >  1:10744584/1‑151 (MQ=255)
                                                                                                      tCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGAt             <  2:3473166/151‑1 (MQ=255)
                                                                                                       cTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATg            <  1:1706371/151‑1 (MQ=255)
                                                                                                        tACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGc           >  2:5806597/1‑151 (MQ=255)
                                                                                                           gACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt        >  1:9553598/1‑151 (MQ=255)
                                                                                                            aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGtt       >  1:5989020/1‑151 (MQ=255)
                                                                                                             ccGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttc      <  2:6294535/151‑1 (MQ=255)
                                                                                                                          |||                                                                                                                                             
TCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCC‑AGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCT  >  NC_008463/2787740‑2788003

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