Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,252,257 T→G 11.2% intergenic (‑181/+59) 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,2570TG11.2% 86.7 / 2.7 35intergenic (‑181/+59)PP_3726/rocEenoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE
Reads supporting (aligned to +/- strand):  ref base T (17/14);  new base G (2/2);  total (19/16)
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.93e-01

GAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAA  >  NC_002947/4252119‑4252395
                                                                                                                                          |                                                                                                                                          
gaGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAg                                                                                                                                 <  1:18274/150‑1 (MQ=255)
    gagaTACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                             >  1:547198/1‑150 (MQ=255)
       aTACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGatc                                                                                                                                       >  2:130157/1‑135 (MQ=37)
              aCGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCa                                                                                                                   <  2:547198/150‑1 (MQ=255)
                    gtTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTTTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAcc                                                                                                             >  2:481307/1‑150 (MQ=255)
                           aaGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCt                                                                                                      >  2:166583/1‑150 (MQ=255)
                            aGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTg                                                                                                     >  1:387881/1‑150 (MQ=255)
                                       gCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAg                                                                                          >  2:521165/1‑150 (MQ=255)
                                                 cTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTc                                                                                <  2:106094/150‑1 (MQ=255)
                                                       acaACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAAccc                                                                          >  1:275642/1‑150 (MQ=255)
                                                                 gtCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCt                                                                <  2:275642/150‑1 (MQ=255)
                                                                  tCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                             <  1:463371/88‑1 (MQ=255)
                                                                  tCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                             >  2:463371/1‑88 (MQ=255)
                                                                     ggCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGt                                                            <  1:521165/150‑1 (MQ=255)
                                                                             acaAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGc                                                    >  1:131020/1‑150 (MQ=255)
                                                                             acaAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGc                                                    >  1:487680/1‑150 (MQ=255)
                                                                               aaGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGccc                                                  >  1:339033/1‑150 (MQ=255)
                                                                               aaGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGccc                                                  <  1:14745/150‑1 (MQ=255)
                                                                                        cGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGAc                                         <  2:487680/150‑1 (MQ=255)
                                                                                          tCATGACGGCTTTATGTCCGGTGAGATGGCGCGCTGCGTAGGGGGCGCTCGACCGCAGAGGGGCTGAAAATGGAGCGAGCAACACCTGACGACGACGAGCTGCGGGGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCGTGCCGCcgg                                       >  2:309009/1‑150 (MQ=255)
                                                                                            aTGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGAcggcg                                     <  2:339033/150‑1 (MQ=255)
                                                                                               aCGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                                  <  1:166583/150‑1 (MQ=255)
                                                                                                  gCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGAc                               >  1:474906/1‑150 (MQ=255)
                                                                                                  gCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGAc                               >  1:474854/1‑150 (MQ=255)
                                                                                                       aTGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTgg                                                                      <  1:104419/106‑1 (MQ=255)
                                                                                                       aTGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTgg                                                                      >  2:104419/1‑106 (MQ=255)
                                                                                                            cTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAc                     >  1:380843/1‑150 (MQ=37)
                                                                                                            cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcggc                                                                                   <  2:201439/88‑1 (MQ=255)
                                                                                                            cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcggc                                                                                   >  1:201439/1‑88 (MQ=255)
                                                                                                                gaggTCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcgccg                 <  1:580568/150‑1 (MQ=38)
                                                                                                                          aCCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                                  >  2:281782/1‑123 (MQ=255)
                                                                                                                          aCCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                                  >  2:281776/1‑123 (MQ=255)
                                                                                                                          aCCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                                  <  1:281782/123‑1 (MQ=255)
                                                                                                                          aCCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                                  <  1:281776/123‑1 (MQ=255)
                                                                                                                               ctacgcggtgCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCaa  <  2:380843/146‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
GAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAAGGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCATAGCACAGGGCAA  >  NC_002947/4252119‑4252395

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