Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 685,193 T→C 7.9% S352G (AGC→GGC) ‡ 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,1930TC 7.9% 177.4 / 4.7 63S352G (AGC→GGC) ‡cadA‑IIcadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base C (2/3);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ACGGTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGG  >  NC_002947/685047‑685336
                                                                                                                                                  |                                                                                                                                                 
acGGTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGt                                                                                                                                                 >  6:169750/1‑149 (MQ=255)
  ggTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCa                                                                                                                                               >  2:71060/1‑149 (MQ=255)
     cGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑g                                                                                                                                           <  7:243599/149‑1 (MQ=255)
      gACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑gg                                                                                                                                          <  6:81742/149‑1 (MQ=255)
      gACCAGCTTCTGGATCGGCGCGTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑gg                                                                                                                                          <  5:168179/149‑1 (MQ=255)
        ccAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑ggtg                                                                                                                                        >  1:189031/1‑147 (MQ=255)
        ccAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑ggtg                                                                                                                                        <  2:189031/149‑3 (MQ=255)
            cTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCt                                                                                                                                    >  8:98492/1‑149 (MQ=255)
            cTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCt                                                                                                                                    >  8:1651/1‑149 (MQ=255)
               cTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTg                                                                                                                                  <  8:141415/148‑1 (MQ=255)
               cTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTg                                                                                                                                  <  5:169750/148‑1 (MQ=255)
                     cGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTg                                                                                                                           >  5:95701/1‑149 (MQ=255)
                      ggCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTgg                                                                                                                          >  3:189010/1‑149 (MQ=255)
                       gcgcTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTgg                                                                                                                          >  6:2370/1‑148 (MQ=255)
                       gcgcTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTgg                                                                                                                          >  1:270896/1‑148 (MQ=255)
                        cgcTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCa                                                                                                                        <  1:236365/149‑1 (MQ=255)
                            ttGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGct                                                                                                                    >  6:225993/1‑149 (MQ=255)
                                          gTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCTCCGGTTACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCa                                                                                                      <  7:167131/149‑1 (MQ=255)
                                          gTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCt                                                                                                                                                   >  6:81256/1‑105 (MQ=255)
                                          gTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCt                                                                                                                                                   <  5:81256/105‑1 (MQ=255)
                                             ttCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑ggg                                                                                                  >  6:75175/1‑149 (MQ=255)
                                                 accaGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTg                                                                                                                           <  1:127741/121‑1 (MQ=255)
                                                 accaGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTg                                                                                                                           >  2:127741/1‑121 (MQ=255)
                                                         gatgatGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCAt                                                                                      <  8:142736/149‑1 (MQ=255)
                                                             atgCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGgct                                                                                  <  5:131879/149‑1 (MQ=255)
                                                               gcgcgcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGc                                                                                                 <  4:55735/132‑1 (MQ=255)
                                                               gcgcgcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGc                                                                                                 >  3:55735/1‑132 (MQ=255)
                                                               gcgcgcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGgctgc                                                                                >  8:248972/1‑149 (MQ=255)
                                                                   gcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGcc                                                                            <  3:107158/149‑1 (MQ=255)
                                                                   gcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGcc                                                                            <  6:259126/149‑1 (MQ=255)
                                                                       gCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCt                                                                        >  7:78369/1‑149 (MQ=255)
                                                                        cACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCtt                                                                       <  4:189010/149‑1 (MQ=255)
                                                                          cGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCg                                                                     <  3:171664/149‑1 (MQ=255)
                                                                            gTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGAc                                                                   <  5:2370/149‑1 (MQ=255)
                                                                             tctcGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGAcc                                                                  >  8:245750/1‑149 (MQ=255)
                                                                             tctcGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGAc                                                                   >  3:131342/1‑148 (MQ=255)
                                                                                      ccAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTTTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCATAGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTc                                                          <  1:188398/149‑1 (MQ=255)
                                                                                          cgcCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCg                                                     <  6:95701/149‑1 (MQ=255)
                                                                                           gcCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑gg                                                                                                   >  3:33973/1‑102 (MQ=255)
                                                                                           gcCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑gg                                                                                                   <  4:33973/102‑1 (MQ=255)
                                                                                                  tACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGt                                             >  3:73259/1‑149 (MQ=255)
                                                                                                   aCGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGtt                                            >  7:44073/1‑149 (MQ=255)
                                                                                                    cGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGtt                                            <  8:78369/148‑1 (MQ=255)
                                                                                                     gCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCg                                          >  6:67744/1‑149 (MQ=255)
                                                                                                          agcagcCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCggg                                     <  5:67744/149‑1 (MQ=255)
                                                                                                          agcagcCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCggg                                     <  4:73259/149‑1 (MQ=255)
                                                                                                               ccGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACCGGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGATTg                                 >  1:51950/1‑149 (MQ=255)
                                                                                                               ccGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGa                                >  4:141322/1‑149 (MQ=255)
                                                                                                                  gCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCtt                                  <  7:163275/144‑1 (MQ=255)
                                                                                                                  gCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCtt                                  >  8:163275/1‑144 (MQ=255)
                                                                                                                        gCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTAc                         >  6:239676/1‑147 (MQ=255)
                                                                                                                          cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        >  4:2913/1‑82 (MQ=255)
                                                                                                                          cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        >  4:6767/1‑82 (MQ=255)
                                                                                                                          cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        <  3:6767/82‑1 (MQ=255)
                                                                                                                          cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        <  3:2913/82‑1 (MQ=255)
                                                                                                                           aTTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCGTGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGt                    >  2:206712/1‑149 (MQ=255)
                                                                                                                           aTTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCg                                                                     >  2:231849/1‑100 (MQ=255)
                                                                                                                           aTTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCg                                                                     <  1:231849/100‑1 (MQ=255)
                                                                                                                               aTGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGcc                <  2:51950/149‑1 (MQ=255)
                                                                                                                                  gcgcCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCa              <  7:1651/148‑1 (MQ=255)
                                                                                                                                     ccACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCAcgcg          <  7:245750/149‑1 (MQ=255)
                                                                                                                                     aacccGGTGACAGCGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCAcgcg          >  6:189571/4‑149 (MQ=255)
                                                                                                                                           gTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCCCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGCGAAGCGTTCGCCCGCCTTGACCAGTACCAGGTCGCCCACGCCCAAGtc    >  8:125467/1‑148 (MQ=255)
                                                                                                                                            tGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTgg   <  4:131342/149‑1 (MQ=255)
                                                                                                                                              aCGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGGCGACTGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTggg  <  1:71060/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACGGTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGG  >  NC_002947/685047‑685336

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: ‑