Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,254,096 G→A 10.2% L196L (CTG→TTG)  PP_4630 ← MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,254,0960GA10.2% 103.8 / 3.6 39L196L (CTG→TTG) PP_4630MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base A (2/2);  total (20/19)
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.64e-01

TGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCATTCGCTGGCCTGCCCATAGCCGCCT  >  NC_002947/5253953‑5254244
                                                                                                                                               |                                                                                                                                                    
tgctgTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTc                                                                                                                                                >  2:403085/1‑150 (MQ=255)
     tgGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCt                                                                                                                                           >  1:210124/1‑150 (MQ=255)
        gCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCAt                                                                                                                                        >  1:511741/1‑150 (MQ=255)
          gACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATgg                                                                                                                                      <  2:541494/150‑1 (MQ=255)
             tGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCAc                                                                                                                                   <  1:166561/150‑1 (MQ=255)
              gATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACt                                                                                                                                  <  2:511741/150‑1 (MQ=255)
                 gggAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGcc                                                                                                                               >  2:175485/1‑150 (MQ=255)
                  ggAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGccc                                                                                                                              >  2:82528/1‑150 (MQ=255)
                          gCACCAGGTCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCAGGCaaaa                                                                                                                      <  2:519856/150‑3 (MQ=255)
                                gggCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCa                                                                                                                                         >  1:846572/1‑125 (MQ=255)
                                gggCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCa                                                                                                                                         <  2:846572/125‑1 (MQ=255)
                                       cGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCAc                                                                                                         <  1:403085/150‑1 (MQ=255)
                                        ggCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATg                                                                                                                                       <  1:54520/119‑1 (MQ=255)
                                        ggCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATg                                                                                                                                       >  2:54520/1‑119 (MQ=255)
                                         gCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCaa                                                                                                                                              >  2:782646/1‑111 (MQ=255)
                                         gCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCaa                                                                                                                                              <  1:782646/111‑1 (MQ=255)
                                         gCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACgg                                                                                                       >  2:387807/1‑150 (MQ=255)
                                                 cGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTa                                                                                               <  1:82528/150‑1 (MQ=255)
                                                          aCCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGcaca                                                                                      <  1:387807/150‑1 (MQ=255)
                                                                 cAAACTACTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCa                                                                               >  2:11108/1‑150 (MQ=255)
                                                                            aGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCgg                                                                    <  2:210124/150‑1 (MQ=255)
                                                                             gCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAACAACTCAAGCTCa                                                                                                                                         >  1:242097/1‑80 (MQ=255)
                                                                             gCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAACAACTCAAGCTCa                                                                                                                                         <  2:242097/80‑1 (MQ=255)
                                                                             gCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAACAACTCAAGCTCa                                                                                                                                         <  2:242194/80‑1 (MQ=255)
                                                                             gCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAACAACTCAAGCTCa                                                                                                                                         >  1:242194/1‑80 (MQ=255)
                                                                                  tGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAg                                                                                                                                                      >  1:625194/1‑62 (MQ=255)
                                                                                  tGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAg                                                                                                                                                      <  2:625194/62‑1 (MQ=255)
                                                                                        ttCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGa                                                        <  1:294664/150‑1 (MQ=255)
                                                                                            ccAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACcgc                                                    >  1:772040/1‑150 (MQ=255)
                                                                                             cAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCg                                                                                                                            >  2:449214/1‑77 (MQ=255)
                                                                                             cAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCg                                                                                                                            <  1:449214/77‑1 (MQ=255)
                                                                                                    aTCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAAc                                            >  2:731164/1‑150 (MQ=255)
                                                                                                                  gtcgCTGCCCAGGATCAGGCTGGTTACCAGCAACTCCAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAg                              <  2:522960/150‑1 (MQ=255)
                                                                                                                     gCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCa                           >  1:697702/1‑150 (MQ=255)
                                                                                                                         cccAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCATTCg                       <  1:532418/150‑1 (MQ=255)
                                                                                                                         cccAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCATTCg                       >  2:532418/1‑150 (MQ=255)
                                                                                                                                    gCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCaa                                                                        >  2:177541/1‑90 (MQ=255)
                                                                                                                                    gCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCaa                                                                        <  1:177541/90‑1 (MQ=255)
                                                                                                                                       ggTTACCAGCAACTCAAGCTCAGGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGAGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCATTCGCTGGCCTGCCCATa         <  2:317121/150‑1 (MQ=255)
                                                                                                                                              aGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGAGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCATTCGCTGGCCTGCCCATAGCCGCCt  >  1:278379/1‑150 (MQ=255)
                                                                                                                                              aGCAACTCAAGCTCATGACACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCATTCGCTGGCCTGCCCATAGCCGCCt  >  1:793150/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAGGATCAGGCTGGTTACCAGCAACTCAAGCTCATGGCACTGCCCCGGCAATGGCGCCAATAGCACGGTGCGGCTACACGGCACACGGGCCAGCTGCAAGCGGGCGAAAACCCGACCGCGCAGGAACTGGTCATAGAACAGCCATTCGCTGGCCTGCCCATAGCCGCCT  >  NC_002947/5253953‑5254244

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