Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 685,327 C→T 6.8% R307H (CGC→CAC) ‡ cadA‑II ← cadmium translocating P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947685,3270CT 6.8% 215.9 / 3.9 73R307H (CGC→CAC) ‡cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base T (3/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

GCCACCGGTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAGG  >  NC_002947/685179‑685470
                                                                                                                                                    |                                                                                                                                               
gcCACCGGTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCAcgc                                                                                                                                                 <  7:99319/149‑1 (MQ=255)
  acccGGTGACAGCGTCACCAGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCAcgcg                                                                                                                                                >  2:264738/3‑148 (MQ=255)
    ccGGTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACgcgc                                                                                                                                               >  4:81425/1‑147 (MQ=255)
       gTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTg                                                                                                                                          >  5:177223/1‑149 (MQ=255)
        tGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTgg                                                                                                                                         <  3:163603/149‑1 (MQ=255)
         gACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTggg                                                                                                                                        >  5:1898/1‑149 (MQ=255)
                tCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGGCGCCCACGCGCAAGTGGGCGATg                                                                                                                                   <  6:166454/147‑1 (MQ=255)
                   ccGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGCGCGATGGCCAc                                                                                                                              >  8:33754/1‑149 (MQ=255)
                      ggTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTc                                                                                                                           <  5:184288/149‑1 (MQ=255)
                        ttGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCct                                                                                                                         <  4:120531/149‑1 (MQ=255)
                        ttGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACAACGCCGTCGACCGGGAAGCGTTAGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCct                                                                                                                         >  6:242051/1‑149 (MQ=255)
                          gCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCctct                                                                                                                       >  1:223745/1‑149 (MQ=255)
                            ttGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTc                                                                                                                     >  8:94696/1‑149 (MQ=255)
                            ttGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTc                                                                                                                     >  6:50151/1‑149 (MQ=255)
                            ttGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTc                                                                                                                     >  6:184994/1‑149 (MQ=255)
                              gggCACTGGCGGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCa                                                                                                                   >  2:191254/1‑149 (MQ=255)
                                gCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACg                                                                                                                 <  4:91764/149‑1 (MQ=255)
                                    tGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTAGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGt                                                                                                             >  6:285481/1‑149 (MQ=255)
                                       cAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGa                                                                                                          <  1:232412/149‑1 (MQ=255)
                                                gccgcTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCCTTCGCCGGCCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACgcgc                                                                                                  >  8:133926/1‑148 (MQ=255)
                                                    cTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGgc                                                                                             >  3:70084/1‑149 (MQ=255)
                                                       aTCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGgcgc                                                                                           <  3:81425/148‑1 (MQ=255)
                                                         cAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGtt                                                                                        >  3:152240/1‑149 (MQ=255)
                                                            ggCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGGTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCg                                                                                     >  8:120435/1‑149 (MQ=255)
                                                            ggCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCg                                                                                     >  4:173390/1‑149 (MQ=255)
                                                                  aTCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGgcg                                                                                 >  2:99676/1‑147 (MQ=255)
                                                                        aTGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCtt                                                                         <  4:152240/149‑1 (MQ=255)
                                                                        aTGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCtt                                                                         >  4:299634/1‑149 (MQ=255)
                                                                           gctgctGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGa                                                                      <  6:1898/149‑1 (MQ=255)
                                                                           gctgctGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGa                                                                      <  5:184994/149‑1 (MQ=255)
                                                                           gctgctGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGa                                                                      <  5:50151/149‑1 (MQ=255)
                                                                             tgctgcCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGcgcg                                                                                          <  3:10366/127‑1 (MQ=255)
                                                                             tgctgcCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGcgcg                                                                                          >  4:10366/1‑127 (MQ=255)
                                                                                tgcCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCt                                                                 >  1:172155/1‑149 (MQ=255)
                                                                                  ccGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCg                                                               <  7:94696/149‑1 (MQ=255)
                                                                                             cACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTgg                                                    >  6:62819/1‑149 (MQ=255)
                                                                                                gCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCgg                                                 <  5:242051/149‑1 (MQ=255)
                                                                                                  cGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTc                                               >  1:55545/1‑149 (MQ=255)
                                                                                                   gTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCt                                              <  5:62819/149‑1 (MQ=255)
                                                                                                    tCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTg                                             <  2:223745/149‑1 (MQ=255)
                                                                                                      gACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCGACGCGCAAGTGGGCGATGGCGACGTCCTGTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGc                                           <  7:33754/149‑1 (MQ=255)
                                                                                                      gACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTgg                                            <  3:299634/148‑1 (MQ=255)
                                                                                                      gACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGc                                           >  8:3513/1‑149 (MQ=255)
                                                                                                       aCCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCg                                          <  1:191254/149‑1 (MQ=255)
                                                                                                          gggAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGttt                                       <  1:12091/149‑1 (MQ=255)
                                                                                                          gggAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGttt                                       >  7:119742/1‑149 (MQ=255)
                                                                                                            gAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTgg                                     >  8:131471/1‑149 (MQ=255)
                                                                                                            gAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTgg                                     <  7:131471/149‑1 (MQ=255)
                                                                                                               gCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGgcgc                                  <  7:133926/149‑1 (MQ=255)
                                                                                                                 gTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCgg                                >  5:293377/1‑149 (MQ=255)
                                                                                                                 gTTCGCCGGGCTTGACCAGTACCACGTCGCCCACCCGCAAGTGGGCGATCGCCACGTCCTCTTCCACGCCGTCGACCACCCGCAACCCGCGATCCCCGCGCAGGCCTTCGAGCGCTCCGATCGCACTGGCGGTCTCGCGTTTGGcgcgc                                >  5:198776/1‑148 (MQ=255)
                                                                                                                    cGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGt                                                <  3:170989/130‑1 (MQ=255)
                                                                                                                    cGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGt                                                >  4:170989/1‑130 (MQ=255)
                                                                                                                        gggCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCa                         <  6:177223/149‑1 (MQ=255)
                                                                                                                         ggCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAg                        <  3:173390/149‑1 (MQ=255)
                                                                                                                         ggCTTGACCAGTACCAGGGCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAg                        <  1:264738/149‑1 (MQ=255)
                                                                                                                          gCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTcc                          <  7:127749/146‑1 (MQ=255)
                                                                                                                          gCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTcc                          >  8:127749/1‑146 (MQ=255)
                                                                                                                            ttGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTGTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTa                     >  5:8976/1‑149 (MQ=255)
                                                                                                                               aCCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACtt                  >  4:55899/1‑149 (MQ=255)
                                                                                                                                 cAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGTCCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGc                >  8:142941/1‑149 (MQ=255)
                                                                                                                                   gTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGccc              <  4:70084/149‑1 (MQ=255)
                                                                                                                                      ccAGGTCGCCCAGGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGGTCCGGGCGGAGGGCTTCGCGGGCTCGGATGGGACGGGCGGTCTGGCGGGTGGCGCGGGGTTGCAGGTACtagagcg             >  4:245204/1‑141 (MQ=255)
                                                                                                                                      ccAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTAc                    <  3:245204/140‑1 (MQ=255)
                                                                                                                                       cAGGTCGCCCACGCGCAAGTGGTCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGc                >  7:47606/1‑143 (MQ=255)
                                                                                                                                         ggTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACtt                  <  5:64415/139‑1 (MQ=255)
                                                                                                                                         ggTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACtt                  >  6:64415/1‑139 (MQ=255)
                                                                                                                                          gTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATa       >  5:66350/1‑149 (MQ=255)
                                                                                                                                          atCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAg                        >  6:247024/2‑132 (MQ=255)
                                                                                                                                          atCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAg                        <  5:247024/131‑1 (MQ=255)
                                                                                                                                          atCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACtt                  >  8:106966/2‑138 (MQ=255)
                                                                                                                                          atCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACtt                  <  7:106966/137‑1 (MQ=255)
                                                                                                                                          atCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATa       >  3:300038/2‑149 (MQ=255)
                                                                                                                                               ccACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCTGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAgg  <  1:99676/149‑1 (MQ=255)
                                                                                                                                               ccACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAgg  <  7:3513/149‑1 (MQ=255)
                                                                                                                                               ccACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAgg  <  2:172155/149‑1 (MQ=255)
                                                                                                                                                   gcgcAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAgg                       <  4:116333/124‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCCACCGGTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAGG  >  NC_002947/685179‑685470

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