Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,882 +G pseudogene (960/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,8811.G99.4% 1336.2 / ‑1.1 346pseudogene (959/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (158/186);  total (159/187)
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.72e-01

Only 100 of 346 total aligned reads displayed.

CATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAAC  >  NC_008463/2787741‑2788023
                                                                                                                                                |                                                                                                                                              
 aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                                                                                       <  2:7362746/150‑1 (MQ=255)
      gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcct                                                                                                                                  <  2:5468100/150‑2 (MQ=255)
      gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcct                                                                                                                                  <  1:3518454/150‑2 (MQ=255)
        tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctgg                                                                                                                                >  1:5834707/1‑147 (MQ=255)
             gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGaa                                                                                                                          >  2:2520642/1‑150 (MQ=255)
             gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGaa                                                                                                                          <  2:238083/150‑1 (MQ=255)
             gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGaa                                                                                                                          <  1:4226987/150‑1 (MQ=255)
              agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                                                                         >  1:5928284/1‑150 (MQ=255)
              agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                                                                         >  2:6066050/1‑150 (MQ=255)
              agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                                                                         <  1:9007050/150‑1 (MQ=255)
                aaCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggc                                                                                                                       >  1:1525988/1‑150 (MQ=255)
                 aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcg                                                                                                                      <  2:10544842/150‑1 (MQ=255)
                  cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgc                                                                                                                     >  1:8465127/1‑150 (MQ=255)
                    gCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcgg                                                                                                                   <  1:8784703/150‑1 (MQ=255)
                     ccGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGGCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggc                                                                                                                  >  2:2719125/1‑150 (MQ=255)
                        ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                                                                               >  2:9146031/1‑150 (MQ=255)
                        ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                                                                               <  1:462662/150‑1 (MQ=255)
                         gAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTg                                                                                                              >  1:4189471/1‑150 (MQ=255)
                           aGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGt                                                                                                            >  1:1332103/1‑150 (MQ=255)
                           aGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGt                                                                                                            >  1:5534068/1‑150 (MQ=255)
                            gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTg                                                                                                           >  1:486634/1‑150 (MQ=255)
                               gCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCGGGTGGAAGGCGCGGCGCTGGTGATc                                                                                                        <  2:7558973/150‑1 (MQ=255)
                                 ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgc                                                                                                      <  1:7131391/150‑1 (MQ=255)
                                 ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgc                                                                                                      >  1:6526658/1‑150 (MQ=255)
                                  caggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                                                                                     >  1:10655092/4‑150 (MQ=255)
                                     gcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                                                                  <  1:3013901/150‑1 (MQ=255)
                                     gcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                                                                  >  2:43186/1‑150 (MQ=255)
                                        gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctc                                                                                               >  2:341248/1‑150 (MQ=255)
                                          cggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCtt                                                                                             >  2:2693145/3‑150 (MQ=255)
                                            gCCAGCGCCTGGAGGATATAAACCTGAACATGCCAAAGGGGTTAAATGCTATCACCGTCTACGACAGAACCACCCTGGTGGACAAGGCCGTGGCCCCGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGTTGATCGCGGTGCTCTTcc                                                                                           >  2:8986126/1‑150 (MQ=255)
                                            gCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                                           <  1:4017739/150‑1 (MQ=255)
                                            gCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTATTcc                                                                                           <  2:7472561/150‑1 (MQ=255)
                                            gCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCACGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                                           <  1:2702332/150‑1 (MQ=255)
                                             ccAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                                                                          <  1:884592/150‑1 (MQ=255)
                                              cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGGGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                                                                         >  2:7178780/1‑150 (MQ=255)
                                               aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTCCGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                                                                        <  2:1578251/150‑1 (MQ=255)
                                                gcgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                                       <  2:10401194/150‑1 (MQ=255)
                                                gcgcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                                       >  1:5138047/1‑150 (MQ=255)
                                                  gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                                                                     >  1:7989712/1‑150 (MQ=255)
                                                  gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTcc                                                                                     <  1:6764262/150‑1 (MQ=255)
                                                    cTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTc                                                                                   >  1:2790472/1‑150 (MQ=255)
                                                     tGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                                                                  >  1:588301/1‑150 (MQ=255)
                                                      ggaggaGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                                                                 <  1:3466797/150‑1 (MQ=255)
                                                          gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                                                                             >  1:5252028/1‑150 (MQ=255)
                                                          gagaTCAACCGGACCCTGCCCACGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                                                                             <  1:757411/150‑1 (MQ=255)
                                                           agaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                                                                            <  2:7096278/150‑1 (MQ=255)
                                                             aTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATc                                                                          >  2:10135768/1‑150 (MQ=255)
                                                               cAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                                                                        <  1:4446315/150‑1 (MQ=255)
                                                                aaCCGGACCCTGCCCACGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgc                                                                       <  2:924953/150‑1 (MQ=255)
                                                                 aCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                                                                      <  2:4600119/150‑1 (MQ=255)
                                                                     gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                                                  <  2:6467795/150‑1 (MQ=255)
                                                                         cTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg                                                              <  1:9152094/150‑1 (MQ=255)
                                                                         cTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg                                                              >  2:3078752/1‑150 (MQ=255)
                                                                         cTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTg                                                              <  1:9587439/150‑1 (MQ=255)
                                                                             cccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATc                                                           <  2:5603957/146‑1 (MQ=255)
                                                                              ccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACCAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa                                                          >  1:2200040/4‑150 (MQ=255)
                                                                              ccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa                                                          >  2:8605641/4‑150 (MQ=255)
                                                                               caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc                                                         <  1:7434992/148‑1 (MQ=255)
                                                                               caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc                                                         <  1:10048566/148‑1 (MQ=255)
                                                                               caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc                                                         <  2:3684793/148‑1 (MQ=255)
                                                                                aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc                                                        >  2:1897310/2‑150 (MQ=255)
                                                                                 aGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                                       >  2:2865611/1‑150 (MQ=255)
                                                                                       gtcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacc                                                  <  2:8852100/144‑1 (MQ=255)
                                                                                        tcaaggCGATCACCGACTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccacca                                                 <  2:8470136/145‑1 (MQ=255)
                                                                                          aaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatc                                               >  2:6800238/4‑150 (MQ=255)
                                                                                           aggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatca                                              <  1:10278439/148‑1 (MQ=255)
                                                                                           aggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatca                                              <  2:2150900/148‑1 (MQ=255)
                                                                                            ggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatcat                                             >  1:7077215/2‑150 (MQ=255)
                                                                                             gCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCatcatc                                            <  2:9657067/150‑1 (MQ=255)
                                                                                             gCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCatcatc                                            <  2:10778305/150‑1 (MQ=255)
                                                                                              cGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATcc                                           <  2:5346791/150‑1 (MQ=255)
                                                                                                aTCACCGTCTACGACCGCACCACCCTGGTGGACCAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCg                                         >  2:9513347/1‑150 (MQ=255)
                                                                                                aTCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCg                                         >  1:10593202/1‑150 (MQ=255)
                                                                                                aTCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCg                                         >  2:4043851/1‑150 (MQ=255)
                                                                                                   aCCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTg                                      <  1:7781386/150‑1 (MQ=255)
                                                                                                    ccGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGt                                     <  2:1117084/150‑1 (MQ=255)
                                                                                                    ccGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGt                                     <  2:5826106/150‑1 (MQ=255)
                                                                                                     cGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                                    >  1:7608656/1‑150 (MQ=255)
                                                                                                      gTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACTTCCGCGCGGCGCTGATCACCGCCACCTTCATCCCGCTGTCg                                   >  2:8028157/1‑150 (MQ=255)
                                                                                                      gTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCg                                   >  2:9367035/1‑150 (MQ=255)
                                                                                                         tACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATg                                >  1:3306279/1‑150 (MQ=255)
                                                                                                         tACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATg                                >  2:8231780/1‑150 (MQ=255)
                                                                                                            gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg                             >  1:9855832/1‑150 (MQ=255)
                                                                                                            gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg                             <  1:6194962/150‑1 (MQ=255)
                                                                                                            gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg                             <  2:278746/150‑1 (MQ=255)
                                                                                                            gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg                             <  1:6862451/150‑1 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGCAGAAGAACCTGGTGGACGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                            <  2:7801057/150‑1 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                            >  1:6255412/1‑150 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                            >  2:7984054/1‑150 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                            <  1:8107247/150‑1 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt                            >  1:8588437/1‑150 (MQ=255)
                                                                                                                gCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttca                         <  1:381514/150‑1 (MQ=255)
                                                                                                                  accaccCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttcacc                       >  1:1674886/1‑150 (MQ=255)
                                                                                                                                 caaggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGt         >  1:769902/8‑150 (MQ=255)
                                                                                                                                  aaggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTc        <  2:7706985/144‑1 (MQ=255)
                                                                                                                                  aaggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTc        <  2:6876877/144‑1 (MQ=255)
                                                                                                                                   aggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCg       <  2:3232958/145‑1 (MQ=255)
                                                                                                                                    ggccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCgg      <  2:4829072/146‑1 (MQ=255)
                                                                                                                                     gccgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGc     >  2:5103243/4‑150 (MQ=255)
                                                                                                                                       cgTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGCTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCCTCCCGCTGTCGATGATGTTCACCTTCACCGGGATGGTCGGCaa   >  1:4900877/2‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCTTCACCGGGATGGTCGGCAAC  >  NC_008463/2787741‑2788023

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: