Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,866,358 C→G 9.1% G174R (GGG→CGG)  PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,866,3580CG 9.1% 135.6 / 3.4 44G174R (GGG→CGG) PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base G (2/2);  total (21/23)
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.82e-01

CCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTCGAAGCCGGTCAGCAAAATGCCGCAGGGCT  >  NC_002947/3866211‑3866504
                                                                                                                                                   |                                                                                                                                                  
ccTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCg                                                                                                                                                   >  2:200174/1‑149 (MQ=255)
  tCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGc                                                                                                                                                 >  5:136158/1‑149 (MQ=255)
    cAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAACTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGcc                                                                                                                                                <  6:80451/148‑1 (MQ=255)
           gCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGAtgctg                                                                                                                                        >  4:135893/1‑149 (MQ=255)
               ccTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCgag                                                                                                                                    <  2:20497/149‑1 (MQ=255)
                   cAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGaa                                                                                                                                  <  1:200174/147‑1 (MQ=255)
                   cAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAg                                                                                                                                 <  3:214963/148‑1 (MQ=255)
                        cAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGc                                                                                                                                <  4:182974/144‑1 (MQ=255)
                        cAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGc                                                                                                                                >  3:182974/1‑144 (MQ=255)
                                  cAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGccc                                                                                                                 >  3:15136/1‑149 (MQ=255)
                                  cAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGccc                                                                                                                 >  3:19279/1‑149 (MQ=255)
                                   aaaGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGcccc                                                                                                                <  2:713/149‑1 (MQ=255)
                                    aaGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCt                                                                                                               <  6:136158/149‑1 (MQ=255)
                                       ccTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGc                                                                                                            >  1:96000/1‑149 (MQ=255)
                                               aaaCCCTCGGCTCGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAg                                                                                                                                 <  6:271791/120‑1 (MQ=255)
                                               aaaCCCTCGGCTCGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAg                                                                                                                                 >  5:271791/1‑120 (MQ=255)
                                                 aCCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGa                                                                                                  <  4:19279/149‑1 (MQ=255)
                                                 aCCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGa                                                                                                  <  4:15136/149‑1 (MQ=255)
                                                   ccTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAACACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTg                                                                                                              >  2:297185/1‑135 (MQ=255)
                                                   ccTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAACACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTg                                                                                                              <  1:297185/135‑1 (MQ=255)
                                                         cACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCcatc                                                                                          <  3:135893/149‑1 (MQ=255)
                                                                      ccTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGcacac                                                                                >  3:105432/1‑146 (MQ=255)
                                                                      ccTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGcacac                                                                                <  4:105432/146‑1 (MQ=255)
                                                                          gggTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCa                                                                          <  1:213610/148‑1 (MQ=255)
                                                                                aGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCt                                                                   >  7:2488/1‑149 (MQ=255)
                                                                                 gCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTg                                                                  <  6:289907/149‑1 (MQ=255)
                                                                                        ggCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTg                                                                  <  8:51099/142‑1 (MQ=255)
                                                                                        ggCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTg                                                                  >  7:51099/1‑142 (MQ=255)
                                                                                                         gcTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGgcagca                                          >  7:298874/1‑149 (MQ=255)
                                                                                                           ttGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGc                                                                                                                                >  1:265059/1‑61 (MQ=255)
                                                                                                           ttGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGc                                                                                                                                <  2:265059/61‑1 (MQ=255)
                                                                                                                 ccGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTgcg                                  >  3:27901/1‑149 (MQ=255)
                                                                                                                         gATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTgcg                                  <  2:37133/141‑1 (MQ=255)
                                                                                                                         gATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTgcg                                  >  1:37133/1‑141 (MQ=255)
                                                                                                                         gATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCAGCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTCGAAGc                          >  3:223849/1‑149 (MQ=255)
                                                                                                                             aGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTCGAAGCCgg                       <  8:2488/148‑1 (MQ=255)
                                                                                                                                      gcagcatcaGGCCGGGCGTGCTGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTc                               >  7:64757/1‑131 (MQ=255)
                                                                                                                                      gcagcatcaGGCCGGGCGTGCTGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTc                               >  7:69098/1‑131 (MQ=255)
                                                                                                                                      gcagcatcaGGCCGGGCGTGCTGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTc                               <  8:64757/131‑1 (MQ=255)
                                                                                                                                      gcagcatcaGGCCGGGCGTGCTGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTc                               <  8:69098/131‑1 (MQ=255)
                                                                                                                                                gCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGgca                                             <  3:205624/107‑1 (MQ=255)
                                                                                                                                                gCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGgca                                             >  4:205624/1‑107 (MQ=255)
                                                                                                                                                gCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTgcg                                  >  6:125754/1‑118 (MQ=255)
                                                                                                                                                gCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTgcg                                  <  5:125754/118‑1 (MQ=255)
                                                                                                                                                  ccGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTCGAAGCCGGTCAGCAAAATGCCGCAGGGCt  >  4:78751/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTGCAGCAGCACGCCCGGCCTGATGCTGCGAGAAGCCGACCGACAAAGCCCCTGGCTCGCCATTGAGCTCCATCATGTGCACACACGGCACGCCGCTGGCAACCAGCATCTGCCGGGCAGCATCGCTGCGCTCGAAGCCGGTCAGCAAAATGCCGCAGGGCT  >  NC_002947/3866211‑3866504

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