Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,787,829 2 bp→AG pseudogene (907‑908/3090 nt) PA14_RS13085 → CusA/CzcA family heavy metal efflux RND transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,787,8290CA99.7% 1361.9 / NA 328pseudogene (907/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (136/191);  total (137/191)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084632,787,8300TG99.7% 1270.6 / NA 328pseudogene (908/3090 nt)PA14_RS13085CusA/CzcA family heavy metal efflux RND transporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (136/191);  total (137/191)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

CTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCC‑AGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACC  >  NC_008463/2787685‑2787975
                                                                                                                                                  ||                                                                                                                                                  
cTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGa                                                                                                                                                <  2:2444854/151‑1 (MQ=255)
   cgcACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCa                                                                                                                                             >  1:10581485/1‑151 (MQ=255)
     cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCACCCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCAcc                                                                                                                                           <  2:3004888/151‑1 (MQ=255)
     cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCAcc                                                                                                                                           >  1:9960293/1‑151 (MQ=255)
     cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCAcc                                                                                                                                           <  2:6786817/151‑1 (MQ=255)
     cACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCAcc                                                                                                                                           <  2:3910632/151‑1 (MQ=255)
      aCCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCTTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCg                                                                                                                                          >  1:2323362/1‑151 (MQ=255)
      aCCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTGCCCC‑AGGGGGTCAAGGCGATCACCg                                                                                                                                          >  2:8937652/1‑151 (MQ=255)
       ccGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGAGTTCATGCTGAACGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGt                                                                                                                                         <  2:2300117/151‑1 (MQ=255)
         ggCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCt                                                                                                                                       <  2:4128864/151‑1 (MQ=255)
          gcgcGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTa                                                                                                                                      <  2:2948582/151‑1 (MQ=255)
           cgcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTAc                                                                                                                                     <  1:5544886/151‑1 (MQ=255)
            gcgGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACg                                                                                                                                    <  2:5331224/151‑1 (MQ=255)
             cgGCCACCGAGAACCGCCGCGAAGTGGTTCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCATGGAGGAGATCAACCGGACCCAG‑CCCAAGGGGGTCAATGCGAACACCGTCAACGa                                                                                                                                   >  2:1930871/1‑151 (MQ=255)
              ggCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACCGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGAc                                                                                                                                  <  2:8633302/151‑1 (MQ=255)
               gCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGAcc                                                                                                                                 <  2:4995970/151‑1 (MQ=255)
                ccACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCg                                                                                                                                <  2:941086/151‑1 (MQ=255)
                ccACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCg                                                                                                                                >  2:7675824/1‑151 (MQ=255)
                 cACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGc                                                                                                                               <  2:1674024/151‑1 (MQ=255)
                 cACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGc                                                                                                                               >  2:4263447/1‑151 (MQ=255)
                        aaCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCAccc                                                                                                                        <  1:3675703/151‑1 (MQ=255)
                         aCGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCt                                                                                                                       >  1:3366625/1‑151 (MQ=255)
                            gCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggt                                                                                                                    >  1:3950285/1‑149 (MQ=255)
                            gCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggt                                                                                                                    <  2:2212257/151‑3 (MQ=255)
                            gCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggt                                                                                                                    >  2:5156056/1‑149 (MQ=255)
                              cgcgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtgg                                                                                                                  >  1:9609198/1‑147 (MQ=255)
                                cgAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggac                                                                                                                >  2:574978/1‑145 (MQ=255)
                                 gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggaca                                                                                                               >  2:7484310/1‑144 (MQ=255)
                                 gAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTggtggaca                                                                                                               <  1:6031213/151‑8 (MQ=255)
                                                        cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaag                                                                                       <  1:6880827/151‑1 (MQ=255)
                                                        cATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaag                                                                                       <  1:3118703/151‑1 (MQ=255)
                                                         aTGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTgaagaaga                                                                                      >  2:1256231/1‑151 (MQ=255)
                                                           gCTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAc                                                                                    <  1:9469592/151‑1 (MQ=255)
                                                            cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCACGGCGATCACCGTCTACGACCGCACCACCCAGGAGGACCAGGCCGGGGCCACG‑GTGACGAAGAAcc                                                                                   >  2:3604697/1‑151 (MQ=255)
                                                            cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCTCCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcc                                                                                   <  2:5957017/151‑1 (MQ=255)
                                                            cTGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcc                                                                                   <  1:7538448/151‑1 (MQ=255)
                                                             tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcct                                                                                  >  2:955136/1‑150 (MQ=255)
                                                             tGATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcct                                                                                  <  1:10041038/151‑2 (MQ=255)
                                                              gATCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctg                                                                                 <  1:10003861/151‑3 (MQ=255)
                                                               aTCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGGCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctgg                                                                                >  1:1272528/1‑148 (MQ=255)
                                                                tCGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAAcctggt                                                                               >  1:4152046/1‑147 (MQ=255)
                                                                 cGGCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtg                                                                             <  2:9025394/151‑1 (MQ=255)
                                                                  ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCTGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                            >  2:3334578/1‑151 (MQ=255)
                                                                  ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                            <  1:114591/151‑1 (MQ=255)
                                                                  ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                            >  1:5842987/1‑151 (MQ=255)
                                                                  ggCGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCtggtgg                                                                            >  1:5302885/1‑151 (MQ=255)
                                                                    cGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGaa                                                                          <  2:731611/151‑1 (MQ=255)
                                                                    cGAGAACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGaa                                                                          <  1:488011/151‑1 (MQ=255)
                                                                     gagaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAg                                                                         >  2:8113192/1‑151 (MQ=255)
                                                                      agaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAgg                                                                        >  2:6307830/1‑151 (MQ=255)
                                                                       gaACAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggc                                                                       >  2:2514209/1‑151 (MQ=255)
                                                                         aCAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgc                                                                     <  1:10302604/151‑1 (MQ=255)
                                                                          cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcg                                                                    <  1:6344766/151‑1 (MQ=255)
                                                                          cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcg                                                                    <  1:8620559/151‑1 (MQ=255)
                                                                          cAGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcg                                                                    <  1:10720594/151‑1 (MQ=255)
                                                                           aGCCGGGAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAggcgcgg                                                                   >  1:8100775/1‑151 (MQ=255)
                                                                               gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCt                                                               >  2:1123761/1‑151 (MQ=255)
                                                                               gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCt                                                               <  1:10382255/151‑1 (MQ=255)
                                                                               gggAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCt                                                               >  1:1357754/1‑151 (MQ=255)
                                                                                 gAAGTGGCCCAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTgg                                                             <  1:2784008/151‑1 (MQ=255)
                                                                                        cccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcgc                                                      <  1:2965159/151‑1 (MQ=255)
                                                                                         ccAGGCGGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATcgcg                                                     >  1:501488/1‑151 (MQ=255)
                                                                                            ggcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTg                                                  >  2:6913548/2‑151 (MQ=255)
                                                                                             gcgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGc                                                 <  1:4798226/151‑1 (MQ=255)
                                                                                              cgGTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGct                                                <  2:3471417/151‑1 (MQ=255)
                                                                                                gTCGGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGAACGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGctct                                              >  2:1371082/1‑151 (MQ=255)
                                                                                                  cggCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTc                                            >  2:10372323/3‑151 (MQ=255)
                                                                                                   ggCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTcc                                           <  2:9011278/150‑1 (MQ=255)
                                                                                                    gCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGACGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                          >  2:8460820/1‑151 (MQ=255)
                                                                                                    gCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCt                                          <  1:5426187/151‑1 (MQ=255)
                                                                                                     ccAGCGCCTGGAGGAGATCAACCGGACCCCG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTg                                         <  1:198305/151‑1 (MQ=255)
                                                                                                      cAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGt                                        <  1:8531054/151‑1 (MQ=255)
                                                                                                       aGCGCCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGtt                                       <  2:7967673/151‑1 (MQ=255)
                                                                                                          gcCTGGATGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGACAACGACCGCACCCCCCTGGTGGACCAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTTGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCACTTCCTGATCCt                                    >  2:2727998/1‑151 (MQ=255)
                                                                                                          gcCTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCt                                    >  2:3736966/1‑151 (MQ=255)
                                                                                                            cTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                  >  1:253060/1‑151 (MQ=255)
                                                                                                            cTGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCg                                  >  1:446473/1‑151 (MQ=255)
                                                                                                             tGGAGGAGATCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCgg                                 <  1:7389737/151‑1 (MQ=255)
                                                                                                                  gagaTCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACa                            >  1:4628139/1‑151 (MQ=255)
                                                                                                                    gaTCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATc                          <  2:4936398/151‑1 (MQ=255)
                                                                                                                     aTCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAG‑ACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATcc                         <  2:6698425/151‑1 (MQ=255)
                                                                                                                      tCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                        <  1:8034908/151‑1 (MQ=255)
                                                                                                                      tCAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcg                        <  1:3286432/151‑1 (MQ=255)
                                                                                                                       cAACCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgc                       <  1:8832390/151‑1 (MQ=255)
                                                                                                                        aaCCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                      >  2:9312897/1‑151 (MQ=255)
                                                                                                                        aaCCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcg                      >  2:1829282/1‑151 (MQ=255)
                                                                                                                         aCCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                     <  1:7247423/151‑1 (MQ=255)
                                                                                                                         aCCGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCcgcgc                     <  1:10599991/151‑1 (MQ=255)
                                                                                                                           cGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg                   <  2:9123325/151‑1 (MQ=255)
                                                                                                                           cGGACCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCgcgg                   <  1:8305942/151‑1 (MQ=255)
                                                                                                                              aCCCTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGgcgc                >  1:6653232/1‑151 (MQ=255)
                                                                                                                               cccTG‑CCCAAGGGGGTCAAGGCGCTCACCGTCTACGACCGCACCCCCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCt               <  1:4282628/151‑1 (MQ=255)
                                                                                                                               cccTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCt               <  1:6234169/151‑1 (MQ=255)
                                                                                                                                 cTG‑CCCAAGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGa             >  1:172559/1‑151 (MQ=255)
                                                                                                                                      ccaaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAc         >  1:9003486/4‑151 (MQ=255)
                                                                                                                                       caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc        <  1:9817336/149‑1 (MQ=255)
                                                                                                                                       caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc        <  1:7029955/149‑1 (MQ=255)
                                                                                                                                       caaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCAcc        <  1:1200030/149‑1 (MQ=255)
                                                                                                                                        aaGGGGGTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCg       >  1:6693095/2‑151 (MQ=255)
                                                                                                                                            gggTCAAGGCGATCACCGTCTACGACCGCACCACCCTGGTGGACAAGG‑CCGTGGCCACGGTGAAGAAGAA‑CCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGGTCCTCGGCAACATCCGCGCGGCGCTGATCACCGccac   <  2:6196489/151‑1 (MQ=255)
                                                                                                                                                  ||                                                                                                                                                  
CTGCGCACCGGCGCGGCCACCGAGAACGGCCGCGAAGTGGTGCTCGGCACGGTGTTCATGCTGATCGGCGAGAACAGCCGGGAAGTGGCCCACGCGGTCTGCCAGCGCCTGGAGGAGATCAACCGGACCCTG‑CCC‑AGGGGGTCACTGCGATCACCGTCTACGACCGCACCACCCTGCATGACTAGGACAGTGCCCAC‑GTGAAGAAGAACCCTGGTGGAAGGCGCGGCGCTGGTGATCGCGGTGCTCTTCCTGTTCCTCGGCAACATCCGCGCGGCGCTGATCACCGCCACC  >  NC_008463/2787685‑2787975

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