Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,252,241 A→G 26.1% intergenic (‑165/+75) PP_3726 ← / ← rocE enoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,252,2410AG26.1% 29.8 / 9.3 23intergenic (‑165/+75)PP_3726/rocEenoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE
Reads supporting (aligned to +/- strand):  ref base A (11/6);  new base G (3/3);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

GTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCAT  >  NC_002947/4252104‑4252383
                                                                                                                                         |                                                                                                                                              
gatgtaggcCTTGAGGTGGGAGATTCCAAGCGCGTGGTCAGCAACGCGGGCTTGGCTGTTACAACTTCGCTCAACTCGGTGTCTGGCGTTACACAAGTGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAggcgg                                                                                                                                    <  2:143245/145‑1 (MQ=255)
 ttGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTAcgcggt                                                                                                                                   >  2:176893/1‑149 (MQ=255)
   ggaggaCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGgcgc                                                                                                                                 <  2:257675/150‑1 (MQ=255)
     aggaCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTAcgcggtgcgc                                                                                                                               >  1:606528/1‑145 (MQ=255)
               gaGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTCACAg                                                                                                                     <  1:176893/150‑1 (MQ=37)
                 gggAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGgc                                                                                                                   <  2:606528/150‑1 (MQ=37)
                                        gCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGCCCTCACAGGCGCTGAAAATGCAGCGACCAacac                                                                                            >  1:350971/1‑150 (MQ=255)
                                        gCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAacac                                                                                            >  2:735268/1‑150 (MQ=255)
                                                      gCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAg                                                                              >  2:355256/1‑150 (MQ=255)
                                                              aaCTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTACGCGgtgc                                                                                                                                 >  1:355034/1‑91 (MQ=25)
                                                              aaCTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTACGCGgtgc                                                                                                                                 <  2:355034/91‑1 (MQ=25)
                                                                cTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACcgc                                                                                                                                            <  1:675572/78‑1 (MQ=255)
                                                                cTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACcgc                                                                                                                                            >  2:675572/1‑78 (MQ=255)
                                                                      acaACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAccc                                                              >  1:90168/1‑150 (MQ=255)
                                                                                 tCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTg                                                   >  2:32176/1‑150 (MQ=255)
                                                                                  cTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTgg                                                  <  1:735268/150‑1 (MQ=255)
                                                                                      cGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGc                                              <  2:90168/150‑1 (MQ=255)
                                                                                           cacaAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAgg                                         >  2:8849/1‑150 (MQ=255)
                                                                                                              aCGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCCAACCCCTCGTGTGCTGGGTGCGAAGGCCCCCTGCCGAcggcggct                      >  2:129795/1‑148 (MQ=255)
                                                                                                              aCGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCCACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                      <  2:350971/150‑1 (MQ=255)
                                                                                                              aCGGCTTTATGGCCTGTGAGATCGCGCACCGCGGAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAGCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                      >  2:216809/1‑150 (MQ=255)
                                                                                                                   tttATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACga                 >  2:47938/1‑150 (MQ=255)
                                                                                                                                  aTCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCAt  >  1:229108/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                              
GTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCAT  >  NC_002947/4252104‑4252383

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