Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 685,330 A→G 6.3% L306S (TTG→TCG) ‡ 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,3300AG 6.3% 227.3 / 4.5 80L306S (TTG→TCG) ‡cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base A (41/33);  new base G (3/2);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAGGGCAAT  >  NC_002947/685186‑685475
                                                                                                                                                |                                                                                                                                                 
gTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTg                                                                                                                                               >  7:100681/1‑149 (MQ=255)
   aCGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACTTGCGCGTGGGc                                                                                                                                            <  5:48708/149‑1 (MQ=255)
   aCGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTggg                                                                                                                                             <  1:70754/148‑1 (MQ=255)
         tCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATg                                                                                                                                        <  8:170972/147‑1 (MQ=255)
             cGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCACGTCGCCAACGCGCAACTGGGCGATGGcccc                                                                                                                                   >  7:80271/1‑146 (MQ=255)
                        ggCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCa                                                                                                                        <  8:209950/148‑1 (MQ=255)
                        ggCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCa                                                                                                                        >  7:209950/1‑148 (MQ=255)
                         gCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGc                                                                                                                                            <  7:202008/127‑1 (MQ=255)
                         gCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGc                                                                                                                                            >  8:202008/1‑127 (MQ=255)
                               gCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTACACGCCGTCg                                                                                                                >  7:136643/1‑149 (MQ=255)
                                    ctctCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCctc                                                                                                                             >  6:295570/1‑131 (MQ=255)
                                    ctctCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGGGGGCGATGGCCACGTCctc                                                                                                                             <  5:295570/131‑1 (MQ=255)
                                     tctcGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTTCCAGGTCGCCCACGCGCAAGTGGTCGATGGCCACGTCCTCTTCCACGCCTTCGACAAcc                                                                                                          >  4:212455/1‑148 (MQ=255)
                                       tcGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCAcgcg                                                                                                        >  6:84294/1‑149 (MQ=255)
                                       tcGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCAcgcg                                                                                                        >  7:143831/1‑149 (MQ=255)
                                         gccgcTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACgcgc                                                                                                       >  7:77804/1‑148 (MQ=255)
                                               gATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTc                                                                                                                                <  6:224583/117‑1 (MQ=255)
                                                aTCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGCGCAGCGTTCGCCGGGCTTGACCAGTACAAGGTCGCCGCCGCGCAAGTGGGCGCTGGCCACGTCGGCTTCCACGCCGTCGACCACGCGCAAGGcgcg                                                                                               >  6:250648/1‑149 (MQ=255)
                                                 tCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGt                                                                                              >  1:46458/1‑149 (MQ=255)
                                                 tCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGt                                                                                              >  7:293175/1‑149 (MQ=255)
                                                    gggCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTcc                                                                                           <  7:74408/149‑1 (MQ=255)
                                                     ggCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCAc                                                                                                                                   >  3:260049/1‑108 (MQ=255)
                                                     ggCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCAc                                                                                                                                   <  4:260049/108‑1 (MQ=255)
                                                       ccTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCggg                                                                                        <  4:71123/149‑1 (MQ=255)
                                                       ccTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCggg                                                                                        >  6:320742/1‑149 (MQ=255)
                                                        cTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGgc                                                                                       >  6:217209/1‑149 (MQ=255)
                                                         tCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCAc                                                                                                                       <  5:78242/116‑1 (MQ=255)
                                                         tCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCAc                                                                                                                       >  6:78242/1‑116 (MQ=255)
                                                               gCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGc                                                                                >  1:186552/1‑149 (MQ=255)
                                                                   ggCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTc                                                                             >  2:249444/1‑148 (MQ=255)
                                                                      tgctgcCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCAATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGCgg                                                                         >  4:319275/1‑149 (MQ=255)
                                                                           ccGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTc                                                                     <  8:143831/148‑1 (MQ=255)
                                                                           ccGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTc                                                                     <  3:241092/148‑1 (MQ=255)
                                                                            cGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAg                                                                          <  4:297196/142‑1 (MQ=255)
                                                                            cGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAg                                                                          >  3:297196/1‑142 (MQ=255)
                                                                                  cGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCgg                                                                                         <  6:218945/121‑1 (MQ=255)
                                                                                  cGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATg                                                                <  2:186552/146‑1 (MQ=255)
                                                                                     ccACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTg                                                          >  2:195639/1‑149 (MQ=255)
                                                                                      cACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGa                                                                           <  4:127101/131‑1 (MQ=255)
                                                                                      cACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGa                                                                           >  3:127101/1‑131 (MQ=255)
                                                                                      cACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTgg                                                         <  1:222145/149‑1 (MQ=255)
                                                                                      cACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTgg                                                         >  4:314123/1‑149 (MQ=255)
                                                                                         gCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAg                                                                          <  8:262220/129‑1 (MQ=255)
                                                                                         gCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAg                                                                          >  7:262220/1‑129 (MQ=255)
                                                                                            gTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCt                                                   <  5:202067/149‑1 (MQ=255)
                                                                                               gACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTg                                                  <  5:101560/147‑1 (MQ=255)
                                                                                               gACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTg                                                  >  6:101560/1‑147 (MQ=255)
                                                                                               gACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGc                                                <  5:250648/149‑1 (MQ=255)
                                                                                                aCCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCg                                               >  8:33594/1‑149 (MQ=255)
                                                                                                aCCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCg                                               <  5:84294/149‑1 (MQ=255)
                                                                                                 ccGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAggg                                                                        <  1:202470/123‑1 (MQ=255)
                                                                                                 ccGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAggg                                                                        >  2:202470/1‑123 (MQ=255)
                                                                                                   gggAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGttt                                            <  5:320742/149‑1 (MQ=255)
                                                                                                     gAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTgg                                          >  8:32287/1‑149 (MQ=255)
                                                                                                        gCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGgcgc                                       <  8:136643/149‑1 (MQ=255)
                                                                                                                cGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTcc                               <  1:249444/149‑1 (MQ=255)
                                                                                                                cGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTc                                >  2:277364/1‑148 (MQ=255)
                                                                                                                  ggCTTGACCATTTCCAGGTCGCCCACGCGCATGTGGGCGATGGGCACGTGCTCTTGCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAg                             <  3:212455/149‑1 (MQ=255)
                                                                                                                       gACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTa                          <  7:89275/147‑1 (MQ=255)
                                                                                                                        aCCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACtt                       >  3:74753/1‑149 (MQ=255)
                                                                                                                          cAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGc                     >  5:82042/1‑149 (MQ=255)
                                                                                                                            gTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGccc                   <  3:314123/149‑1 (MQ=255)
                                                                                                                             tccgcggTCGCCCACGCGCAAGTGGTCGATGGCGACGTCCTCTTCACCGCCGGCGCCCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTGGAGGGGTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCGAGGTACTTGCCCa                  <  8:80271/144‑1 (MQ=255)
                                                                                                                             tACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCa                  <  5:217209/149‑1 (MQ=255)
                                                                                                                                cAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGc                >  2:312830/1‑148 (MQ=255)
                                                                                                                                 ccatCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTc                                >  4:253725/4‑131 (MQ=255)
                                                                                                                                 ccatCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTc                                <  3:253725/128‑1 (MQ=255)
                                                                                                                                  catCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCaa              >  5:116485/3‑148 (MQ=255)
                                                                                                                                   gTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCa                              <  4:57219/131‑1 (MQ=255)
                                                                                                                                   gTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCa                              >  3:57219/1‑131 (MQ=255)
                                                                                                                                   atCGCCCACTTGCGCGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATa            >  5:167249/2‑149 (MQ=255)
                                                                                                                                    tCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGgcgc                                       >  3:252469/1‑121 (MQ=255)
                                                                                                                                    tCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGgcgc                                       <  4:252469/121‑1 (MQ=255)
                                                                                                                                      gCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATAc           >  6:44218/1‑147 (MQ=255)
                                                                                                                                       cccACGCGCTAGTGGGCGATGGCCACGTAGTCTTACTCTCTGTCGAGCACGCGCAAGGCGCGTTCCGGGCGCTGGGCTTCGAGGGGGCTGATGGCACTGTCGTTCTGGCGTTTGGCTAGGCTTTCCATTTACTTGCCCAGCAATAGCAg        <  8:100681/149‑1 (MQ=255)
                                                                                                                                       cccACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAg        >  8:196580/1‑149 (MQ=255)
                                                                                                                                       cccACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGAAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCAAGCAATACCAg        >  6:44420/1‑149 (MQ=255)
                                                                                                                                        ccACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAgg       >  3:75346/1‑149 (MQ=255)
                                                                                                                                        ccACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAgg       >  1:215463/1‑149 (MQ=255)
                                                                                                                                              gcAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAGGGCAAt  >  4:96139/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTGACGCTGTCACCGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCAGGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGGCGATGGCCACGTCCTCTTCCACGCCGTCGACCACGCGCAAGGCGCGTTCCGGGCGCAGGGCTTCGAGGGCTCGGATGGCACTGGCGGTCTGGCGTTTGGCGCGGCTTTCCAGGTACTTGCCCAGCAATACCAGGGCAAT  >  NC_002947/685186‑685475

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: