Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,871 Δ1 bp pseudogene (949/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,8710A.98.4% 1734.7 / 23.3 366pseudogene (949/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (151/207);  minor base G (3/3);  total (154/212)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 366 total aligned reads displayed.

TCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCT  >  NC_008463/2787740‑2788003
                                                                                                                                     |                                                                                                                                     
tCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaa                                                                                                                    >  1:6004205/1‑151 (MQ=255)
 cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                                                                   <  1:6558935/151‑1 (MQ=255)
 cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaag                                                                                                                   >  2:10495071/1‑151 (MQ=255)
  aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACCAGGCCGTGGCCACGGTgaagaaga                                                                                                                  >  1:2269764/1‑151 (MQ=255)
   tGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaagaa                                                                                                                 >  1:11152008/1‑151 (MQ=255)
   tGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTgaagaagaa                                                                                                                 >  2:7181846/1‑151 (MQ=255)
     cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcc                                                                                                               <  2:9105516/151‑1 (MQ=255)
      tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcct                                                                                                              <  1:2510238/151‑2 (MQ=255)
       gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAAcctg                                                                                                             >  1:7078596/1‑149 (MQ=255)
         tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGACGAAGAAcctggt                                                                                                           >  1:8669180/1‑147 (MQ=255)
          cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtg                                                                                                         >  1:9051667/1‑151 (MQ=255)
          cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCtggtg                                                                                                         >  1:1877478/1‑151 (MQ=255)
             cGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGaa                                                                                                      >  1:7926529/1‑151 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                                                     <  2:11057280/151‑1 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                                                     <  1:8324630/151‑1 (MQ=255)
              gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAg                                                                                                     >  2:9358356/1‑151 (MQ=255)
                 aaCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcg                                                                                                  <  2:10770916/151‑1 (MQ=255)
                  aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgc                                                                                                 <  2:5904088/151‑1 (MQ=255)
                       cGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACCAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcgc                                                                                            >  2:3315976/1‑151 (MQ=255)
                       cGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAggcgcggcgc                                                                                            <  1:3074890/151‑1 (MQ=255)
                        gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                                                           <  2:9169198/151‑1 (MQ=255)
                        gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                                                           >  2:5137113/1‑151 (MQ=255)
                        gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCt                                                                                           <  1:1690200/151‑1 (MQ=255)
                         ggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTg                                                                                          <  2:8571524/151‑1 (MQ=255)
                           aaGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGt                                                                                        >  2:2985228/1‑151 (MQ=255)
                             gTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGa                                                                                      <  1:6779389/151‑1 (MQ=255)
                              tGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGAt                                                                                     <  2:4019266/151‑1 (MQ=255)
                                 cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATcgc                                                                                  >  2:1779938/1‑151 (MQ=255)
                                    aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCTAGGGGGTCTAGGCGATCTCCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt                                                                               <  2:6985228/149‑1 (MQ=255)
                                    aggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGt                                                                               >  2:8412424/3‑151 (MQ=255)
                                     ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                                              <  1:9327438/150‑1 (MQ=255)
                                       ctgtcgtccAGCGACTGGAGGAGATCAACCGGACCCTGCCCTAGGGGGTCAAGCCGTACACCTTCTACGTCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                                                            <  2:2442854/144‑1 (MQ=255)
                                       cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                                                            <  1:9374979/151‑1 (MQ=255)
                                       cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                                                            >  1:4816851/1‑151 (MQ=255)
                                           cggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTc                                                                        <  1:3640308/149‑1 (MQ=255)
                                            ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                       <  1:3446989/150‑1 (MQ=255)
                                            ggCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                                                       >  1:7334245/2‑151 (MQ=255)
                                             gCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                                                      <  1:4073845/151‑1 (MQ=255)
                                               cAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                                                    <  1:492530/151‑1 (MQ=255)
                                                aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                   <  2:6348994/151‑1 (MQ=255)
                                                aGCGCCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                                                   >  2:234575/1‑151 (MQ=255)
                                                   gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                                                                <  1:551028/151‑1 (MQ=255)
                                                   gcCTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACACGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                                                                <  1:8003120/151‑1 (MQ=255)
                                                    ccTGGAGGAGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTc                                                               <  1:10248989/151‑1 (MQ=255)
                                                       ggaggaGATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                                                            <  2:3268069/151‑1 (MQ=255)
                                                       ggaggaGATCAACCGGACCCTGCCCATGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGc                                                            <  1:5127546/151‑1 (MQ=255)
                                                         aggagATCACCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCaa                                                          <  1:7512182/151‑1 (MQ=255)
                                                          ggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                                                         <  2:10022940/151‑1 (MQ=255)
                                                          ggagATCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAAc                                                         <  1:10063984/151‑1 (MQ=255)
                                                           gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                                                        >  1:865002/1‑151 (MQ=255)
                                                           gagaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                                                        >  1:1555727/1‑151 (MQ=255)
                                                            agaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACAt                                                       >  1:5277425/1‑151 (MQ=255)
                                                            agaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACAt                                                       <  2:1526057/151‑1 (MQ=255)
                                                             gaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATc                                                      <  1:5893721/151‑1 (MQ=255)
                                                             gaTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATc                                                      >  2:160863/1‑151 (MQ=255)
                                                              aTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                                                     <  1:9300082/151‑1 (MQ=255)
                                                              aTCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                                                     <  1:3998928/151‑1 (MQ=255)
                                                               tCAACCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                                                    <  2:2658225/151‑1 (MQ=255)
                                                                 aaCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                                                  <  2:3883511/151‑1 (MQ=255)
                                                                  aCCGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                                                 >  2:9549529/1‑151 (MQ=255)
                                                                   ccGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgcg                                                >  1:4166946/1‑151 (MQ=255)
                                                                    cGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg                                               <  1:10888805/151‑1 (MQ=255)
                                                                    cGGACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg                                               >  2:2111437/1‑151 (MQ=255)
                                                                     ggACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggc                                              <  2:4094805/151‑1 (MQ=255)
                                                                      gACCCTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcggcg                                             >  2:5668243/1‑151 (MQ=255)
                                                                        cccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCTGCGCTGGTGGACGCGGTGCTCTACGCGTTCCTCGGCAACAGCCGCGCGGCGCt                                           >  2:11000811/1‑151 (MQ=255)
                                                                        cccTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCt                                           <  1:10750104/151‑1 (MQ=255)
                                                                          cTGCCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGa                                         <  2:4605722/151‑1 (MQ=255)
                                                                             cccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCa                                      >  1:1144846/5‑151 (MQ=255)
                                                                               caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc                                    >  2:5156982/3‑151 (MQ=255)
                                                                                aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg                                   >  2:4249390/2‑151 (MQ=255)
                                                                                 aGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGc                                  <  2:11226875/151‑1 (MQ=255)
                                                                                    ggggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGcca                                >  2:5479729/1‑151 (MQ=255)
                                                                                        tcaaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcat                            <  2:7893488/146‑1 (MQ=255)
                                                                                          aaggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatca                          <  1:9476135/148‑1 (MQ=255)
                                                                                           aggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcatcat                         >  1:9498516/3‑151 (MQ=255)
                                                                                           aggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACcagcat                         >  2:6475901/3‑151 (MQ=255)
                                                                                            ggCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCCCCatcatc                        >  2:5882717/2‑151 (MQ=255)
                                                                                               gATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCg                     >  1:4379011/1‑151 (MQ=255)
                                                                                               gATCACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCg                     <  2:2798219/151‑1 (MQ=255)
                                                                                                  cACCGTCTACGACCGCACCCCCCTGGTGGACAAGTCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGGTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTg                  <  2:10258631/151‑1 (MQ=255)
                                                                                                  cACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTg                  >  2:7599688/1‑151 (MQ=255)
                                                                                                  cACCGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTg                  <  1:1267066/151‑1 (MQ=255)
                                                                                                    ccGTCTACGACCGCACCACCCTGGTGGACTAGGCCGTGGCCACGGTGAAGAAGAACCGGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                <  2:4403304/151‑1 (MQ=255)
                                                                                                    ccGTCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTc                >  1:10574000/1‑151 (MQ=255)
                                                                                                       tCTACGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGAt             <  1:2765172/151‑1 (MQ=255)
                                                                                                         tACGACCGCACCTCCCTGGTGGACATGGCCGTGGCCACGGTGAAGACGAACCGGGTGGAAGGCGCGGCGCGGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGc           <  2:7998376/151‑1 (MQ=255)
                                                                                                          aCGACCGCACCCCCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTTGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGGTCCTCGGCAACTTGCGCGCGGCGCTGATCACCGCCACCTTCTTCCCGCTGTCGATGCt          <  2:7811715/151‑1 (MQ=255)
                                                                                                          aCGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCt          <  1:9660581/151‑1 (MQ=255)
                                                                                                          aCGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCt          >  1:2689024/1‑151 (MQ=255)
                                                                                                           cGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCTCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg         <  1:5647561/151‑1 (MQ=255)
                                                                                                           cGACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTg         >  1:7432246/1‑151 (MQ=255)
                                                                                                            gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt        >  1:1350264/1‑151 (MQ=255)
                                                                                                            gACCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGt        >  2:10629199/1‑151 (MQ=255)
                                                                                                             aCCGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAACGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGtt       <  2:7253107/151‑1 (MQ=255)
                                                                                                              ccGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttc      <  1:2064160/151‑1 (MQ=255)
                                                                                                               cGCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttca     <  1:1241904/151‑1 (MQ=255)
                                                                                                                gCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttcac    <  1:6795927/151‑1 (MQ=255)
                                                                                                                gCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttcac    <  1:3821567/151‑1 (MQ=255)
                                                                                                                gCACCACCCTGGTGGACAAGGCCGTGGCCACGGTGAAGAAGAACCGGGTGGAATGCGCGGCGCTGGTGATCGCGGTGCTCTTCCGGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGttcac    >  2:3610751/1‑151 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCAGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCACGTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACCATCATCCCGCTGTCGATGCTGTTCACCT  >  NC_008463/2787740‑2788003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: