Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,718,692 T→G 8.7% intergenic (+36/+95) fabB → / ← PP_4176 3‑oxoacyl‑ACP synthase/5‑oxo‑L‑prolinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,718,6920TG 8.7% 124.5 / 3.8 46intergenic (+36/+95)fabB/PP_41763‑oxoacyl‑ACP synthase/5‑oxo‑L‑prolinase
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base G (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

ATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTTGT  >  NC_002947/4718546‑4718839
                                                                                                                                                  |                                                                                                                                                   
aTCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTc                                                                                                                                                   <  4:187326/149‑1 (MQ=255)
                   aaaaCGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGt                                                                                                                                <  7:67755/149‑1 (MQ=255)
                       cGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgttt                                                                                                                                    <  6:304629/141‑1 (MQ=255)
                       cGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgttt                                                                                                                                    >  5:304629/1‑141 (MQ=255)
                        gCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCt                                                                                                                           <  4:23835/149‑1 (MQ=255)
                                  acacGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGG‑AGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCctc                                                                                                                >  4:306274/1‑149 (MQ=255)
                                      ggTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTc                                                                                                             >  6:305289/1‑149 (MQ=255)
                                          aTGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTcgc                                                                                                           >  8:65858/1‑147 (MQ=255)
                                          aTGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTcgc                                                                                                           <  7:65858/147‑1 (MQ=255)
                                          aTGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGACCTCTTCGCgg                                                                                                         <  3:18348/149‑1 (MQ=255)
                                             aGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGt                                                                                                       >  5:206494/1‑148 (MQ=255)
                                               cAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTaaa                                                                                                    >  1:308555/1‑149 (MQ=255)
                                                aaCAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAc                                                                                                   >  6:23709/1‑149 (MQ=255)
                                                              cggcggCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                       <  2:308555/147‑1 (MQ=255)
                                                               ggcggcACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                       <  5:32500/146‑1 (MQ=255)
                                                               ggcggcACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                       >  6:32500/1‑146 (MQ=255)
                                                                      ccAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAccc                                                                                                 >  8:62208/1‑129 (MQ=255)
                                                                      ccAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAccc                                                                                                 <  7:62208/129‑1 (MQ=255)
                                                                        aaCGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTtttgtttg                                                                                                                                   >  8:265254/1‑93 (MQ=255)
                                                                        aaCGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCa                                                                           <  6:206494/149‑1 (MQ=255)
                                                                         aCGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAg                                                                          <  2:196374/149‑1 (MQ=255)
                                                                          cGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAtt                                                                               <  2:48809/143‑1 (MQ=255)
                                                                          cGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAtt                                                                               >  1:48809/1‑143 (MQ=255)
                                                                                     tGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCa                                                              >  7:301655/1‑149 (MQ=255)
                                                                                            gcgcTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAg                                                                          <  5:217530/130‑1 (MQ=255)
                                                                                            gcgcTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAg                                                                          >  6:217530/1‑130 (MQ=255)
                                                                                            gcgcTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACCCCCCGACTGGTTCGGGGCCTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTGATTTGCAGTGCTTGTGGCCATTGCAgg                                                       >  4:328470/1‑149 (MQ=255)
                                                                                              gcTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                        >  2:106580/1‑114 (MQ=255)
                                                                                              gcTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTaca                                                                                        <  1:106580/114‑1 (MQ=255)
                                                                                                    aaGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTTGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCa                                                                                 <  6:202603/115‑1 (MQ=255)
                                                                                                    aaGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCa                                                                                 >  5:202603/1‑115 (MQ=255)
                                                                                                    aaGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGgc                                               <  5:305289/149‑1 (MQ=255)
                                                                                                            tAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCa                                                                           <  2:9428/113‑1 (MQ=255)
                                                                                                            tAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCa                                                                           >  1:9428/1‑113 (MQ=255)
                                                                                                               gACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGc                                                                                                                       <  7:113657/66‑1 (MQ=255)
                                                                                                               gACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGc                                                                                                                       >  8:113657/1‑66 (MQ=255)
                                                                                                                   ctgctgctgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAaca                                <  5:23709/149‑1 (MQ=255)
                                                                                                                      ctgctgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCg                             >  4:130085/1‑149 (MQ=255)
                                                                                                                          tgTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAg                                                                                     >  4:107097/1‑89 (MQ=255)
                                                                                                                                 aaaaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGgcgcgc                                           >  1:66174/2‑124 (MQ=37)
                                                                                                                                 aaaaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGgcgcgc                                           <  2:66174/123‑1 (MQ=37)
                                                                                                                                   aaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAttt                                                                              >  5:239505/1‑87 (MQ=21)
                                                                                                                                   aaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAttt                                                                              <  6:239505/87‑1 (MQ=21)
                                                                                                                                              cTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGc                                                                >  6:360545/1‑90 (MQ=255)
                                                                                                                                                 gTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGaaaa                     <  2:144351/130‑1 (MQ=255)
                                                                                                                                                 gTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTtgt  >  5:351039/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCGCCAGAAAATTCTCCCGGGCCACCTTGT  >  NC_002947/4718546‑4718839

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: ‑