Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,718,690 G→C 9.1% intergenic (+34/+97) 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,6900GC 9.1% 136.9 / 4.6 44intergenic (+34/+97)fabB/PP_41763‑oxoacyl‑ACP synthase/5‑oxo‑L‑prolinase
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base C (2/2);  total (23/21)
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.48e-01

GCCAATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCG  >  NC_002947/4718542‑4718812
                                                                                                                                                    |                                                                                                                          
gCCAATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTg                                                                                                                            >  3:187326/1‑149 (MQ=255)
    aTCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTc                                                                                                                        <  4:187326/149‑1 (MQ=255)
             aaaaCCGAAGAAAACGCCAAGATCGACACGGACATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTg                                                                                                                            >  6:160132/1‑136 (MQ=255)
             aaaaCCGAAGAAAACGCCAAGATCGACACGGACATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTg                                                                                                                            <  5:160132/136‑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                                                            >  6:32500/1‑146 (MQ=255)
                                                                   ggcggcACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                            <  5:32500/146‑1 (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                                                    >  1:48809/1‑143 (MQ=255)
                                                                              cGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAtt                                                    <  2:48809/143‑1 (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                                                >  1:9428/1‑113 (MQ=255)
                                                                                                                tAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCa                                                <  2:9428/113‑1 (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                                                   <  6:239505/87‑1 (MQ=21)
                                                                                                                                       aaaCGCCCCGAACCAGTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCAttt                                                   >  5:239505/1‑87 (MQ=21)
                                                                                                                                                  cTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGc                                     >  6:360545/1‑90 (MQ=255)
                                                                                                                                                    |                                                                                                                          
GCCAATCCTGCGTAAAACCGAAGAAAACGCCAAGATCGACACGGTCATGAGCAACAGCTTCGGCTTCGGCGGCACCAACGCCACCCTGGTGCTCAAGCGCTGGGAAGGCAAGTAAGACGCTGCTGCTGTAAATGAAAACGCCCCGACTGGTTCGGGGCGTTTTTGTTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGTCATTTGCAGTGCTTGTGGCCATTGCAGGCTAATGGCGCGCCTTGGGCAACACCG  >  NC_002947/4718542‑4718812

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