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.1% 1275.7 / 1.5 329pseudogene (939/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (150/176);  minor base . (1/2);  total (151/178)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
*NC_0084632,787,8620AT98.8% 1373.9 / ‑3.4 329pseudogene (940/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base T (149/175);  minor base G (2/2);  total (151/178)
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,8630TG98.5% 1266.7 / ‑1.1 329pseudogene (941/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base G (150/174);  total (151/178)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 329 total aligned reads displayed.

CATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑C‑A‑GGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACC  >  NC_008463/2787741‑2788002
                                                                                                                           |||                                                                                                                                            
   gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaagaa                                                                                                                >  2:5219776/1‑150 (MQ=255)
       aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctg                                                                                                            >  2:197343/1‑148 (MQ=255)
        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctgg                                                                                                           <  1:7583091/150‑4 (MQ=255)
        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctgg                                                                                                           <  2:6971240/150‑4 (MQ=255)
           gCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                                                       >  2:10114309/1‑150 (MQ=255)
           gCGAGAACAGCCGGGAAGAGGCCCAGGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                                                       >  2:7944417/1‑150 (MQ=255)
              agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                                                    >  2:6066050/1‑150 (MQ=255)
              agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                                                    >  1:4233228/1‑150 (MQ=255)
               gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAgg                                                                                                   <  2:3770722/150‑1 (MQ=255)
               gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAgg                                                                                                   >  2:2271167/1‑150 (MQ=255)
                 aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcg                                                                                                 <  2:10544842/150‑1 (MQ=255)
                  cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgc                                                                                                >  1:8465127/1‑150 (MQ=255)
                  cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgc                                                                                                >  1:4869478/1‑150 (MQ=255)
                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                                              >  1:6899719/1‑150 (MQ=255)
                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                                              <  1:3009157/150‑1 (MQ=255)
                     ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGTACCCTTCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGTCAAGG‑CCGTGGCCACGGTGAATAAGAA‑CCTGGTGGAAggcgcggc                                                                                             <  2:3393360/150‑1 (MQ=255)
                     ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcggc                                                                                             >  1:4759921/1‑150 (MQ=255)
                     ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGGCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcggc                                                                                             >  2:2719125/1‑150 (MQ=255)
                        ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCt                                                                                          >  2:9146031/1‑150 (MQ=255)
                          aaGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGTAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTgg                                                                                        <  1:1254768/150‑1 (MQ=255)
                               gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCGGGTGGAAGGCGCGGCGCTGGTGATc                                                                                   <  2:7558973/150‑1 (MQ=255)
                                 ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcgc                                                                                 <  1:7131391/150‑1 (MQ=255)
                                  caggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACAACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                                                                <  1:371883/147‑1 (MQ=255)
                                     gcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                                             >  2:1606220/1‑150 (MQ=255)
                                        gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                                                          >  1:7812334/1‑150 (MQ=255)
                                        gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                                                          >  2:4425992/1‑150 (MQ=255)
                                          cggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCtt                                                                        >  2:2702839/3‑150 (MQ=255)
                                          cggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCtt                                                                        <  1:7706722/148‑1 (MQ=255)
                                            gCCAGCGCCTGGAGGATATAAACCTGAACATGCCAA‑A‑GGGGTTAAATGCTATCACCGTCTACGACAGAACCACCCTGGTGGACAAGGCCGTGGCCCCG‑GTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGTTGATCGCGGTGCTCTTcc                                                                      >  2:8986126/1‑150 (MQ=255)
                                            gCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                      >  1:4118419/1‑150 (MQ=255)
                                             ccAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                                                     <  1:884592/150‑1 (MQ=255)
                                             ccAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                                                     >  2:8297132/1‑150 (MQ=255)
                                             ccAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                                                     <  1:347378/150‑1 (MQ=255)
                                              cAGCGCCTGGAGGTGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTG‑GTGGTGAAGGCCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                                                    <  2:8511011/150‑1 (MQ=255)
                                              cAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                                                    <  1:1569412/150‑1 (MQ=255)
                                              cAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                                                    <  2:1114451/150‑1 (MQ=255)
                                               aGCGCCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGACGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                                                   <  2:5417708/150‑1 (MQ=255)
                                                gcgcCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                  <  1:5230108/150‑1 (MQ=255)
                                                gcgcCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                  <  2:10399304/150‑1 (MQ=255)
                                                 cgcCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTc                                                                 <  2:5532006/150‑1 (MQ=255)
                                                  gcCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                                                <  1:6764262/150‑1 (MQ=255)
                                                  gcCTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                                                >  1:7989712/1‑150 (MQ=255)
                                                   ccTGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                                                               >  2:6694503/1‑150 (MQ=255)
                                                     tGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                             >  2:7483071/1‑150 (MQ=255)
                                                     tGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                             >  2:4195041/1‑150 (MQ=255)
                                                     tGGAGGAGATCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                             >  1:5887480/1‑150 (MQ=255)
                                                          gagaTCAACCGGACCCTGCC‑C‑ACGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                                                        <  1:757411/150‑1 (MQ=255)
                                                           agaTCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                                                       <  2:7096278/150‑1 (MQ=255)
                                                            gaTCAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACAt                                                      <  1:8943002/150‑1 (MQ=255)
                                                               cAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                                                   <  1:622869/150‑1 (MQ=255)
                                                               cAACCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                                                   <  1:5651346/150‑1 (MQ=255)
                                                                aaCCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgc                                                  >  1:3912290/1‑150 (MQ=255)
                                                                 aCCGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                                                 <  2:4600119/150‑1 (MQ=255)
                                                                  ccGGACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                                                >  2:7719828/1‑150 (MQ=255)
                                                                    ggACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg                                              >  2:2936618/1‑150 (MQ=255)
                                                                     gACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                             <  2:6467795/150‑1 (MQ=255)
                                                                     gACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                             >  1:4511815/1‑150 (MQ=255)
                                                                     gACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                             <  1:2736397/150‑1 (MQ=255)
                                                                     gACCCTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                             <  2:2565365/150‑1 (MQ=255)
                                                                       cccTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGgcgc                                           >  2:1290441/1‑150 (MQ=255)
                                                                         cTGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg                                         >  2:3078752/1‑150 (MQ=255)
                                                                          tGCC‑CAA‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGa                                        <  1:10641803/150‑1 (MQ=255)
                                                                             cccaa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATc                                      <  1:4973712/146‑1 (MQ=255)
                                                                             cccaa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATc                                      <  1:8750131/146‑1 (MQ=255)
                                                                              ccaa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa                                     >  2:8605641/4‑150 (MQ=255)
                                                                               caa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc                                    <  1:10048566/148‑1 (MQ=255)
                                                                               caa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc                                    <  1:7434992/148‑1 (MQ=255)
                                                                                aa‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc                                   <  2:7405881/149‑1 (MQ=255)
                                                                                 a‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                  >  1:3452749/1‑150 (MQ=255)
                                                                                 a‑GGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                  <  1:5289186/150‑1 (MQ=255)
                                                                                      gggcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccac                              >  1:4217607/8‑150 (MQ=255)
                                                                                      ggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccac                              <  2:3575726/150‑1 (MQ=255)
                                                                                       gtcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacc                             <  2:8852100/144‑1 (MQ=255)
                                                                                       gtcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacc                             <  1:10503033/144‑1 (MQ=255)
                                                                                        tcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacca                            >  1:1771575/6‑150 (MQ=255)
                                                                                          aaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatc                          <  2:2485559/147‑1 (MQ=255)
                                                                                          aaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatc                          <  1:9574641/147‑1 (MQ=255)
                                                                                            ggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatcat                        <  2:8021973/149‑1 (MQ=255)
                                                                                            ggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatcat                        >  1:7077215/2‑150 (MQ=255)
                                                                                             gCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCatcatc                       <  2:9657067/150‑1 (MQ=255)
                                                                                             gCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCatcatc                       <  1:1461950/150‑1 (MQ=255)
                                                                                              cGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATcc                      <  1:6537335/150‑1 (MQ=255)
                                                                                              cGATCACCGACTACGACCTCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATcc                      <  2:8031171/150‑1 (MQ=255)
                                                                                                aTCACCGTCTACGACCGCACCACCCTGGTGGACCAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCg                    >  2:9513347/1‑150 (MQ=255)
                                                                                                aTCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCg                    <  1:587362/150‑1 (MQ=255)
                                                                                                 tCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGc                   <  2:10753852/150‑1 (MQ=255)
                                                                                                    ccGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGt                <  2:2802838/150‑1 (MQ=255)
                                                                                                     cGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc               <  2:5966375/150‑1 (MQ=255)
                                                                                                     cGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc               <  2:4899442/150‑1 (MQ=255)
                                                                                                      gTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCg              >  2:9367035/1‑150 (MQ=255)
                                                                                                       tCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGa             <  1:3077553/150‑1 (MQ=255)
                                                                                                          aCGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGc          <  2:7210714/150‑1 (MQ=255)
                                                                                                          aCGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGc          <  1:9084648/150‑1 (MQ=255)
                                                                                                            gACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg        >  1:9855832/1‑150 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt       >  1:6255656/1‑150 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt       <  1:8107247/150‑1 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt       >  1:8588437/1‑150 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt       >  1:6038645/1‑150 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt       <  1:2578965/150‑1 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt       >  2:9216694/1‑150 (MQ=255)
                                                                                                                           |||                                                                                                                                            
CATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCC‑C‑A‑GGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACC  >  NC_008463/2787741‑2788002

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