Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 684,846:1 +A 7.4% coding (1401/2400 nt) cadA‑II ← cadmium translocating P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947684,8461.A 7.4% 371.8 / 23.3 95coding (1401/2400 nt)cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base . (41/47);  new base A (3/4);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

GCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCA  >  NC_002947/684699‑684989
                                                                                                                                                     |                                                                                                                                               
gcCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCa                                                                                                                                                <  5:68706/149‑1 (MQ=255)
    gCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCaaggat                                                                                                                                           >  7:122508/1‑144 (MQ=255)
    gCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCGGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGAt                                                                                                                                            <  8:300806/149‑1 (MQ=255)
        ggCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCg                                                                                                                                        >  3:218897/1‑149 (MQ=255)
          cGGTAATGGTCTGCGCTGTTGCCGTCGAGCGCCTGGCTGTGAACCACCTGCGGGCTGCAGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGAGCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGCCgaa                                                                                                                                      <  4:217754/149‑3 (MQ=255)
          cGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTGATCAAGGATGCCgaa                                                                                                                                      <  8:122508/149‑3 (MQ=255)
             tAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTggcg                                                                                                                                   >  5:343381/1‑149 (MQ=255)
                aGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTggcgggc                                                                                                                                <  3:348663/149‑1 (MQ=255)
                 ggTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTggctggcg                                                                                                                               <  4:228255/149‑1 (MQ=255)
                  gTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATACTTGATCAGGATGCCGTGgcgggcgg                                                                                                                              <  1:211727/149‑1 (MQ=255)
                   tCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGggcggc                                                                                                                             <  5:132045/149‑1 (MQ=255)
                   tCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGggcggc                                                                                                                             <  8:268490/149‑1 (MQ=255)
                   tCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGggcggc                                                                                                                             >  3:255501/1‑149 (MQ=255)
                       cgcTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGAcc                                                                                                                         >  1:316809/1‑149 (MQ=255)
                       cgcTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGAcc                                                                                                                         >  2:205409/1‑149 (MQ=255)
                       cgcTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGAcc                                                                                                                         >  8:346425/1‑149 (MQ=255)
                          ttgggCCGTCCCGCGCCTGGCGTTGAACCTCCTGCGGGCTGCCCGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGGTGACCGCATGGGCACGCTCCTGGGCTTCGTCAGCCTTGATCTGGATGCCGTGGCGGGCGGCGACCCCg                                                                                                                      <  4:219239/145‑1 (MQ=255)
                            ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACCACTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCg                                                                                                                      <  1:102723/148‑1 (MQ=255)
                            ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCg                                                                                                                      <  5:376501/147‑1 (MQ=255)
                            ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCg                                                                                                                      <  6:196598/147‑1 (MQ=255)
                            ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGCATTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGAACAGGATCCCGTGGCGGGCGGCGAcccccg                                                                                                                     >  3:255605/1‑146 (MQ=255)
                              gCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTAc                                                                                                                  >  2:78745/1‑149 (MQ=255)
                               ccGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTAcc                                                                                                                 <  2:367190/149‑1 (MQ=255)
                                cGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCg                                                                                                                <  4:255501/149‑1 (MQ=255)
                                 gTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                               >  1:340761/1‑149 (MQ=255)
                                  tCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCg                                                                                                                                        <  6:75855/123‑1 (MQ=255)
                                  tCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCg                                                                                                                                        >  5:75855/1‑123 (MQ=255)
                                  tCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                              <  2:316809/149‑1 (MQ=255)
                                    cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCa                                                                                                            <  6:321038/149‑1 (MQ=255)
                                     aGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATc                                                                                                                                                  >  6:368092/1‑111 (MQ=255)
                                     aGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATc                                                                                                                                                  <  5:368092/111‑1 (MQ=255)
                                     aGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCAt                                                                                                           >  3:240603/1‑149 (MQ=255)
                                        gcCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGAt                                                                                                        >  6:192483/1‑149 (MQ=255)
                                         ccTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                       <  1:301576/149‑1 (MQ=255)
                                         ccTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                       <  5:371569/149‑1 (MQ=255)
                                             gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTg                                                                                                                                      <  7:144493/114‑1 (MQ=255)
                                             gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATcgcg                                                                                                    <  7:229443/148‑1 (MQ=255)
                                             gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCCGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTg                                                                                                                                      >  8:144493/1‑114 (MQ=255)
                                                 tgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGcggcgg                                                                                                <  2:319345/148‑1 (MQ=255)
                                                 tgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCACGGTACCGGCCATGATCGCggcggg                                                                                               >  2:125238/1‑149 (MQ=255)
                                                  gAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCggcgggc                                                                                              >  4:5602/1‑149 (MQ=255)
                                                   aaCCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCggcgggcg                                                                                             >  6:113272/1‑149 (MQ=255)
                                                   aaCCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCggcgggcg                                                                                             <  2:216517/149‑1 (MQ=255)
                                                   aaCCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTGCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCggcgggc                                                                                              <  6:23589/148‑1 (MQ=255)
                                                    accaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGt                                                                                            <  3:56330/149‑1 (MQ=255)
                                                     ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                               >  4:292637/1‑129 (MQ=255)
                                                     ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                               <  3:292637/129‑1 (MQ=255)
                                                     ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                           >  4:273740/1‑149 (MQ=255)
                                                     ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATACTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                           >  4:385341/1‑149 (MQ=255)
                                                         cTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCa                                                                                       <  1:205409/149‑1 (MQ=255)
                                                                 tGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAgg                                                                               <  6:252733/149‑1 (MQ=255)
                                                                         gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAggg                                                                              >  6:21268/1‑142 (MQ=255)
                                                                         gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAggg                                                                              <  5:21268/142‑1 (MQ=255)
                                                                         gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAggg                                                                       <  8:289075/149‑1 (MQ=255)
                                                                            aGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGc                                                                     <  8:18996/148‑1 (MQ=255)
                                                                              ggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATTATCGCGGCGGGCGTGTCCAGGCCGATGGCGCAGGGGCAgg                                                                  >  6:67183/1‑149 (MQ=255)
                                                                              ggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAgg                                                                  <  6:375087/149‑1 (MQ=255)
                                                                               gTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGc                                                                 <  3:5602/149‑1 (MQ=255)
                                                                                  ccGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                              <  4:240603/149‑1 (MQ=255)
                                                                                  ccGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                              <  6:343381/149‑1 (MQ=255)
                                                                                          gTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGc                                                       <  4:218897/148‑1 (MQ=255)
                                                                                           tCGAACACTACGCGGTTGACCGGATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCac                                                     >  8:253689/1‑149 (MQ=255)
                                                                                           tCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCac                                                     >  7:270065/1‑149 (MQ=255)
                                                                                              aaCACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggc                                                  >  3:290917/1‑149 (MQ=255)
                                                                                               acacTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCATTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggag                                                 >  7:30597/1‑147 (MQ=255)
                                                                                                 acTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgac                                               <  5:126585/149‑1 (MQ=255)
                                                                                                     cgcgGTTGACCGCATGGGCACGCTGCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCGCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacggcg                                           <  1:125238/149‑1 (MQ=255)
                                                                                                      gcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGt                                          >  5:244226/1‑149 (MQ=255)
                                                                                                      gcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGt                                          <  8:270065/149‑1 (MQ=255)
                                                                                                       cgGTTGACCGCATGGGCACGCTCCAGGGCTTCTGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGtt                                         <  7:346425/149‑1 (MQ=255)
                                                                                                        ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTg                                        <  5:192483/149‑1 (MQ=255)
                                                                                                          ttGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGAt                                      >  5:122204/1‑149 (MQ=255)
                                                                                                             aCCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAg                                   <  2:67937/149‑1 (MQ=255)
                                                                                                              ccGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgac                                               <  3:282395/136‑1 (MQ=255)
                                                                                                              ccGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgac                                               >  4:282395/1‑136 (MQ=255)
                                                                                                                 cATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGcgcg                               >  5:237962/1‑149 (MQ=255)
                                                                                                                 cATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGcgcg                               <  4:290917/149‑1 (MQ=255)
                                                                                                                    gggCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGtt                            >  2:273221/1‑149 (MQ=255)
                                                                                                                       cACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTc                          >  5:45072/1‑148 (MQ=255)
                                                                                                                        aCGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTc                          >  4:205745/1‑147 (MQ=255)
                                                                                                                             tCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCgg                    >  6:357033/1‑149 (MQ=255)
                                                                                                                               cAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                  <  3:385341/149‑1 (MQ=255)
                                                                                                                               cAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                  >  7:43804/1‑149 (MQ=255)
                                                                                                                                 gggCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAcc                >  6:139558/1‑149 (MQ=255)
                                                                                                                                 gggCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAcc                <  5:113272/149‑1 (MQ=255)
                                                                                                                                     ttCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccag            >  2:253292/1‑149 (MQ=255)
                                                                                                                                     cacGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGtt                            >  6:279241/3‑133 (MQ=255)
                                                                                                                                     cacGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGtt                            <  5:279241/131‑1 (MQ=255)
                                                                                                                                     cacGGCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccag            >  6:359482/3‑149 (MQ=255)
                                                                                                                                      acGGCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc           <  5:359482/148‑1 (MQ=255)
                                                                                                                                       cGGCATCCTGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagcc          >  1:100031/1‑149 (MQ=255)
                                                                                                                                            tCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAgccagcca     >  1:286846/1‑149 (MQ=255)
                                                                                                                                              cctGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAgg                                                                  >  6:56013/3‑87 (MQ=255)
                                                                                                                                              cctGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAgg                                                                  <  5:56013/85‑1 (MQ=255)
                                                                                                                                              cTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccacc   <  7:253689/149‑1 (MQ=255)
                                                                                                                                              cTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGGTTCCAGCGGTACCCCAGccagccacc   >  3:205288/1‑149 (MQ=255)
                                                                                                                                               ttGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccac    >  7:53381/1‑147 (MQ=255)
                                                                                                                                               ttGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccac    >  4:362762/1‑147 (MQ=255)
                                                                                                                                                 gATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccacca  <  5:357033/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATCAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCACCA  >  NC_002947/684699‑684989

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: