Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,784 T→G pseudogene (862/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,7840TG99.7% 1186.0 / NA 310pseudogene (862/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (145/164);  minor base C (1/0);  total (146/164)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

CGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCT  >  NC_008463/2787636‑2787931
                                                                                                                                                    |                                                                                                                                                      
cGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGACCAGCCGGGAAGTGGCCCAGGCGGTCGGc                                                                                                                                                      >  1:10085104/1‑151 (MQ=255)
   gCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAg                                                                                                                                                   <  1:7973730/151‑1 (MQ=255)
    ccGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAgc                                                                                                                                                  >  1:8764503/1‑151 (MQ=255)
        tCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCt                                                                                                                                              <  2:7500484/151‑1 (MQ=255)
          cGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTgg                                                                                                                                            >  1:10738148/1‑151 (MQ=255)
           gTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGACGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCCGCGCCTgga                                                                                                                                           >  2:7700298/1‑151 (MQ=255)
           gTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTgga                                                                                                                                           >  2:1008736/1‑151 (MQ=255)
           gTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTgga                                                                                                                                           >  1:3102384/1‑151 (MQ=255)
            tATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTggag                                                                                                                                          <  1:4837845/151‑1 (MQ=255)
             aTCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTggagg                                                                                                                                         >  2:5386420/1‑151 (MQ=255)
             aTCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTggagg                                                                                                                                         <  2:2676680/151‑1 (MQ=255)
             aTCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTggagg                                                                                                                                         <  1:8682693/151‑1 (MQ=255)
              tCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTggagga                                                                                                                                        >  1:8942734/1‑151 (MQ=255)
                  cgACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATc                                                                                                                                    >  2:10376114/1‑151 (MQ=255)
                  cgACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATc                                                                                                                                    <  1:8617522/151‑1 (MQ=255)
                   gACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCa                                                                                                                                   >  2:1508584/1‑151 (MQ=255)
                      gTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAAcc                                                                                                                                <  1:9474462/151‑1 (MQ=255)
                      gTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAAcc                                                                                                                                <  2:7667828/151‑1 (MQ=255)
                      gTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAAcc                                                                                                                                <  2:311790/151‑1 (MQ=255)
                        cGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCgg                                                                                                                              >  2:1822869/1‑151 (MQ=255)
                          ccGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGAc                                                                                                                            <  2:5712997/151‑1 (MQ=255)
                          ccGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGAc                                                                                                                            <  1:5626045/151‑1 (MQ=255)
                            gAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGAccc                                                                                                                          >  1:10002783/1‑151 (MQ=255)
                             aGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCt                                                                                                                         <  2:3352941/151‑1 (MQ=255)
                             aGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCt                                                                                                                         <  2:6342746/151‑1 (MQ=255)
                                   gCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCaa                                                                                                                 >  1:7584958/1‑150 (MQ=255)
                                    ccTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCaag                                                                                                                >  2:9167333/1‑149 (MQ=255)
                                     cTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAagg                                                                                                               <  2:1125299/151‑1 (MQ=255)
                                      tGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGACGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAaggg                                                                                                              >  2:3377611/1‑151 (MQ=255)
                                        ggCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAaggggg                                                                                                            >  1:6834609/1‑150 (MQ=255)
                                         gCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAagggggt                                                                                                           <  2:7191572/151‑3 (MQ=255)
                                          cAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTc                                                                                                           <  1:9834968/151‑1 (MQ=255)
                                           aaGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCa                                                                                                          <  2:2139814/151‑1 (MQ=255)
                                              gcgCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGCAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaagg                                                                                                       <  2:9907881/149‑4 (MQ=255)
                                              gAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaagg                                                                                                       <  1:8374314/151‑4 (MQ=255)
                                               aGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGTACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaaggc                                                                                                      >  1:4159920/1‑147 (MQ=255)
                                               aGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaaggc                                                                                                      <  1:10597818/151‑5 (MQ=255)
                                               aGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaaggc                                                                                                      <  2:608337/151‑5 (MQ=255)
                                               aGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCaaggc                                                                                                      <  1:5563936/151‑5 (MQ=255)
                                                    cgcACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCa                                                                                                 >  2:9463892/1‑151 (MQ=255)
                                                    cgcACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCa                                                                                                 <  1:9891618/151‑1 (MQ=255)
                                                         cGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTc                                                                                            <  2:3514005/151‑1 (MQ=255)
                                                           gcgcGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTa                                                                                          <  1:10547777/151‑1 (MQ=255)
                                                            cgcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                         >  2:8768518/1‑151 (MQ=255)
                                                              cgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGa                                                                                       >  2:10268082/1‑151 (MQ=255)
                                                              cgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGa                                                                                       >  1:3634066/1‑151 (MQ=255)
                                                               ggCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAc                                                                                      <  1:8030902/151‑1 (MQ=255)
                                                                gCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGAcc                                                                                     <  2:7624430/151‑1 (MQ=255)
                                                                 ccACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGACGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCg                                                                                    >  2:8457447/1‑151 (MQ=255)
                                                                     cGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacc                                                                                >  2:5263247/1‑151 (MQ=255)
                                                                      gagaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcacca                                                                               >  2:8641904/1‑151 (MQ=255)
                                                                       agaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGcaccac                                                                              >  1:6261832/1‑151 (MQ=255)
                                                                        gaACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCaccacc                                                                             >  1:8191577/1‑151 (MQ=255)
                                                                          aCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCt                                                                           <  2:9688858/151‑1 (MQ=255)
                                                                           cGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTg                                                                          <  2:6813747/151‑1 (MQ=255)
                                                                            ggCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTgg                                                                         <  1:5317448/151‑2 (MQ=255)
                                                                             gCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggt                                                                        <  1:4384406/151‑3 (MQ=255)
                                                                                  gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggaca                                                                   >  1:9704978/1‑144 (MQ=255)
                                                                                  gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggaca                                                                   >  1:7333613/1‑144 (MQ=255)
                                                                                  gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggaca                                                                   <  2:4275176/151‑8 (MQ=255)
                                                                                                        tCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCTAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaa                                            <  1:4665094/151‑1 (MQ=255)
                                                                                                         cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                           <  1:3407902/151‑1 (MQ=255)
                                                                                                         cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                           <  1:10368633/151‑1 (MQ=255)
                                                                                                          aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaaga                                          >  2:238084/1‑151 (MQ=255)
                                                                                                          aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaaga                                          >  1:9146177/1‑151 (MQ=255)
                                                                                                              tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGTAGAAGAAcct                                      >  1:4019027/1‑150 (MQ=255)
                                                                                                              tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcct                                      <  1:6070686/151‑2 (MQ=255)
                                                                                                               gATCCGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGAGCCTGGAGGAGATCAAACGGACACTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCTCACCACCCTGGTGTACAAGGCCGTGGCCACGGTgaagaagaaacag                                     >  2:5227427/1‑147 (MQ=255)
                                                                                                                aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                                    >  2:9584025/1‑148 (MQ=255)
                                                                                                                 tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctggt                                   <  2:368029/151‑5 (MQ=255)
                                                                                                                 tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctggt                                   <  1:5077205/151‑5 (MQ=255)
                                                                                                                  cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtg                                 <  1:1243617/151‑1 (MQ=255)
                                                                                                                    gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGa                               <  2:8569974/151‑1 (MQ=255)
                                                                                                                      gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                             <  2:251421/151‑1 (MQ=255)
                                                                                                                       agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAgg                            <  1:1508584/151‑1 (MQ=255)
                                                                                                                        gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggc                           >  1:212386/1‑151 (MQ=255)
                                                                                                                          aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgc                         >  2:4991022/1‑151 (MQ=255)
                                                                                                                           cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcg                        >  2:8170613/1‑151 (MQ=255)
                                                                                                                            aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcgg                       <  1:8309847/151‑1 (MQ=255)
                                                                                                                              ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcg                     <  2:5981272/151‑1 (MQ=255)
                                                                                                                                    aGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTg               <  2:4513031/151‑1 (MQ=255)
                                                                                                                                     gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGa              <  2:3971505/151‑1 (MQ=255)
                                                                                                                                     gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGa              >  2:6048357/1‑151 (MQ=255)
                                                                                                                                      tGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGag             >  2:7115173/1‑150 (MQ=255)
                                                                                                                                       ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGAta            <  2:2783809/151‑2 (MQ=255)
                                                                                                                                       ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc            <  1:573609/151‑1 (MQ=255)
                                                                                                                                       ggCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATc            <  1:3235180/151‑1 (MQ=255)
                                                                                                                                        gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcg           >  2:4449423/1‑151 (MQ=255)
                                                                                                                                          ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgcg         <  2:6685060/151‑1 (MQ=255)
                                                                                                                                          ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgcg         <  2:170252/151‑1 (MQ=255)
                                                                                                                                           caggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCCAGGGGGTCAAGGCGAGCACCGTCTCCGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg        >  2:5504383/4‑151 (MQ=255)
                                                                                                                                           caggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCgg        <  1:9323116/148‑1 (MQ=255)
                                                                                                                                            aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt       <  1:2802459/149‑1 (MQ=255)
                                                                                                                                            aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt       <  1:2509870/149‑1 (MQ=255)
                                                                                                                                             ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGACGGCGCGGCGCTGGTGATCGCGGTg      <  2:3283347/150‑1 (MQ=255)
                                                                                                                                             ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg      >  1:3489236/2‑151 (MQ=255)
                                                                                                                                              gcgGTCGGCCAGCGGCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCTCCCAGGTGGACAAGGCCGTGGCCTCGGTGAAGAAGAACCTGGTGGAAGGCGCGTCGCTGGTGATCGCGGTGc     <  2:7800173/151‑1 (MQ=255)
                                                                                                                                              gcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGc     <  1:1881813/151‑1 (MQ=255)
                                                                                                                                               cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGACGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct    <  2:9029893/151‑1 (MQ=255)
                                                                                                                                                 gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctct  <  1:6357953/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGTGCCGATCCGTATCCGCGACGTCGCCGAGGTCGGCCTGGGCAAGGAGCTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCT  >  NC_008463/2787636‑2787931

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: