Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,252,246 G→C 28.5% intergenic (‑170/+70) 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,2460GC28.5% 35.2 / 10.5 21intergenic (‑170/+70)PP_3726/rocEenoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE
Reads supporting (aligned to +/- strand):  ref base G (10/5);  new base C (3/3);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGG  >  NC_002947/4252104‑4252392
                                                                                                                                              |                                                                                                                                                  
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)
                                                                      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)
                                                                                                                                           cgcgTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAggg  >  2:233495/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGG  >  NC_002947/4252104‑4252392

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