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% 1429.9 / NA 343pseudogene (907/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (163/178);  total (165/178)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084632,787,8300TG100.0% 1294.4 / NA 343pseudogene (908/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (165/178);  total (165/178)

Only 100 of 343 total aligned reads displayed.

CTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACC  >  NC_008463/2787685‑2787975
                                                                                                                                                   ||                                                                                                                                                  
cTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGa                                                                                                                                                >  1:4541786/1‑151 (MQ=255)
   cgcACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCa                                                                                                                                             <  1:7268088/151‑1 (MQ=255)
       ccGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGt                                                                                                                                         <  2:2001519/151‑1 (MQ=255)
         ggCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCt                                                                                                                                       <  2:8413632/151‑1 (MQ=255)
          gcgcGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTa                                                                                                                                      <  2:5182363/151‑1 (MQ=255)
          gcgcGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTa                                                                                                                                      >  2:3823353/1‑151 (MQ=255)
           cgcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                                                                     <  1:990003/151‑1 (MQ=255)
           cgcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                                                                     <  1:2236606/151‑1 (MQ=255)
              ggCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAc                                                                                                                                  <  1:4270379/151‑1 (MQ=255)
              ggCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAc                                                                                                                                  >  2:8182871/1‑151 (MQ=255)
              ggCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAc                                                                                                                                  <  2:8805779/151‑1 (MQ=255)
               gCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAcg                                                                                                                                 <  2:9146055/151‑2 (MQ=255)
               gCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAcc                                                                                                                                 >  2:9397040/1‑151 (MQ=255)
                ccACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCg                                                                                                                                <  2:10192230/151‑1 (MQ=255)
                 cACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGc                                                                                                                               >  1:9335678/1‑151 (MQ=255)
                  aCCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGca                                                                                                                              >  1:6412702/1‑151 (MQ=255)
                    cGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacc                                                                                                                            >  1:6268227/1‑151 (MQ=255)
                     gagaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacca                                                                                                                           <  1:5401992/151‑1 (MQ=255)
                     gagaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacca                                                                                                                           <  2:2474524/151‑1 (MQ=255)
                       gaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCaccacc                                                                                                                         >  2:2885454/1‑151 (MQ=255)
                        aaCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCAccc                                                                                                                        >  1:8978308/1‑151 (MQ=255)
                        aaCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCAccc                                                                                                                        >  2:7739355/1‑151 (MQ=255)
                          cGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTg                                                                                                                      <  2:178833/151‑1 (MQ=255)
                           ggCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTgg                                                                                                                     >  2:2184316/1‑150 (MQ=255)
                            gCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCACGGCGATCACCGTCTACGACCGCACCACCCTggt                                                                                                                    >  2:5849650/1‑149 (MQ=255)
                             ccGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtg                                                                                                                   <  2:10754147/151‑4 (MQ=255)
                             ccGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtg                                                                                                                   <  2:7543656/151‑4 (MQ=255)
                               gcgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtgga                                                                                                                 >  2:10950619/1‑146 (MQ=255)
                               gcgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtgga                                                                                                                 >  2:5385930/1‑146 (MQ=255)
                                 gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGAGCACCGTCTACGACCGCACCACCCTggtggaca                                                                                                               >  2:8574055/1‑144 (MQ=255)
                                                        cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                                       <  1:6136281/151‑1 (MQ=255)
                                                        cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                                       >  2:960287/1‑151 (MQ=255)
                                                        cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                                       <  1:1120719/151‑1 (MQ=255)
                                                        cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                                       >  1:1814197/1‑151 (MQ=255)
                                                            cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcc                                                                                   >  2:7917644/1‑151 (MQ=255)
                                                               aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                                                                                >  1:4352412/1‑148 (MQ=255)
                                                                 cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtg                                                                             >  1:5217398/1‑151 (MQ=255)
                                                                  ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtgg                                                                            >  1:9094748/1‑151 (MQ=255)
                                                                   gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGa                                                                           <  1:1382276/151‑1 (MQ=255)
                                                                     gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                         >  1:6103383/1‑151 (MQ=255)
                                                                          cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCATGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcg                                                                    >  1:3833765/1‑151 (MQ=255)
                                                                            gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggc                                                                  <  2:3106613/151‑1 (MQ=255)
                                                                               gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                               <  2:1593751/151‑1 (MQ=255)
                                                                                ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTg                                                              <  1:6561554/151‑1 (MQ=255)
                                                                                ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTg                                                              <  2:3959988/151‑1 (MQ=255)
                                                                                    gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGa                                                          >  1:10273902/1‑151 (MQ=255)
                                                                                    gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGa                                                          <  1:2051914/151‑1 (MQ=255)
                                                                                     tGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGAt                                                         <  1:6680684/151‑1 (MQ=255)
                                                                                      ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGTGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGATCCTGGTGGAAGGCGCGGCGCTGGTGATc                                                        <  2:834368/151‑1 (MQ=255)
                                                                                      ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGCACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                        <  1:2548384/151‑1 (MQ=255)
                                                                                      ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                        >  1:10889769/1‑151 (MQ=255)
                                                                                      ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                        <  1:3433228/151‑1 (MQ=255)
                                                                                      ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                        <  1:3307536/151‑1 (MQ=255)
                                                                                      ggCCCAGGCGGTCGGCCAGCGCCGGGAGGAGATCAACCGGACCCTGCCCACGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc                                                        <  2:5431717/151‑1 (MQ=255)
                                                                                        cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGACGGCGCGGCGCTGGTGATcgc                                                      >  2:4408272/1‑151 (MQ=255)
                                                                                         ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                                     <  2:7384484/151‑1 (MQ=255)
                                                                                           aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt                                                   <  2:426832/149‑1 (MQ=255)
                                                                                            ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                  <  1:7621284/150‑1 (MQ=255)
                                                                                              cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                                >  1:6985711/1‑151 (MQ=255)
                                                                                               ggTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                               >  2:9779930/1‑151 (MQ=255)
                                                                                                gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctct                                              <  2:6463150/151‑1 (MQ=255)
                                                                                                   ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCCAGGGGGACAAGGCGATCACCGTCGACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGACGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                           <  2:6828515/150‑1 (MQ=255)
                                                                                                   ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                           >  1:2283112/2‑151 (MQ=255)
                                                                                                   ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                           <  1:950662/150‑1 (MQ=255)
                                                                                                      cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                        >  1:10038679/1‑151 (MQ=255)
                                                                                                       aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                       <  1:7368427/151‑1 (MQ=255)
                                                                                                           ccTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTc                                   >  2:11213616/1‑151 (MQ=255)
                                                                                                             tGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGANCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                 <  1:5885928/151‑1 (MQ=255)
                                                                                                              ggaggaGATTAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                                <  1:2943285/151‑1 (MQ=255)
                                                                                                                aggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCaa                              <  1:8447556/151‑1 (MQ=255)
                                                                                                                 ggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                             <  1:8133530/151‑1 (MQ=255)
                                                                                                                 ggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                             >  1:4600982/1‑151 (MQ=255)
                                                                                                                  gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                            <  1:7763764/151‑1 (MQ=255)
                                                                                                                      tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                        >  1:10223544/1‑151 (MQ=255)
                                                                                                                      tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                        >  2:4889192/1‑151 (MQ=255)
                                                                                                                      tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                        <  2:5924594/151‑1 (MQ=255)
                                                                                                                       cAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgc                       <  2:6157448/151‑1 (MQ=255)
                                                                                                                        aaCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                      >  1:7423345/1‑151 (MQ=255)
                                                                                                                         aCCGGACCCTGCCCAAGGTGGTCAAGGCGATAACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                     <  2:2574347/151‑1 (MQ=255)
                                                                                                                         aCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                     >  2:7255518/1‑151 (MQ=255)
                                                                                                                         aCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                     <  2:8447055/151‑1 (MQ=255)
                                                                                                                           cGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCACCATCCGCgcgg                   >  2:8028759/1‑151 (MQ=255)
                                                                                                                            ggACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCGACGACCGCACCACCCTGGTGGACAAGGACGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                  <  2:504240/151‑1 (MQ=255)
                                                                                                                             gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggcg                 <  2:3235610/151‑1 (MQ=255)
                                                                                                                              aCCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGgcgc                <  1:4761517/151‑1 (MQ=255)
                                                                                                                               cccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCt               <  2:2648740/151‑1 (MQ=255)
                                                                                                                                ccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg              >  2:6967023/1‑151 (MQ=255)
                                                                                                                                ccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg              <  1:2844023/151‑1 (MQ=255)
                                                                                                                                ccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg              <  1:3197618/151‑1 (MQ=255)
                                                                                                                                     cccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa          >  2:3566775/5‑151 (MQ=255)
                                                                                                                                     cccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa          >  1:9450467/5‑151 (MQ=255)
                                                                                                                                      ccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc         >  1:3113403/4‑151 (MQ=255)
                                                                                                                                      ccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc         >  1:1642554/4‑151 (MQ=255)
                                                                                                                                       caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc        >  1:7830690/3‑151 (MQ=255)
                                                                                                                                       caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc        <  2:5666404/149‑1 (MQ=255)
                                                                                                                                        aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg       >  2:7755990/2‑151 (MQ=255)
                                                                                                                                           gggggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGcc     >  1:8232925/1‑151 (MQ=255)
                                                                                                                                            ggggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGcca    <  1:11082121/151‑1 (MQ=255)
                                                                                                                                              ggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacc  >  1:4471033/1‑151 (MQ=255)
                                                                                                                                              ggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacc  >  2:9264295/1‑151 (MQ=255)
                                                                                                                                                   ||                                                                                                                                                  
CTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACC  >  NC_008463/2787685‑2787975

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: