Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,872,269 T→C 77.1% S174P (TCA→CCA)  ulaF → L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,872,2690TC77.1% 102.7 / 22.4 78S174P (TCA→CCA) ulaFL‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (11/7);  new base C (30/30);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

GCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAA  >  minE/2872201‑2872315
                                                                    |                                              
gCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCt                                                >  1:2783859/1‑69 (MQ=255)
gCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCt                                                >  1:1600677/1‑69 (MQ=255)
   gAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:944996/69‑1 (MQ=255)
   gAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:2435873/69‑1 (MQ=255)
   gAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:2620115/69‑1 (MQ=255)
   gAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:2898504/69‑1 (MQ=255)
   gAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:549244/69‑1 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGATAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:74650/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:506326/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:340682/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:3059951/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:3034215/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:520097/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:2885186/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:681360/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:2533167/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:984315/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:2370091/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:2156777/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:1881764/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:1614997/1‑68 (MQ=255)
                         ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg                        >  1:1708479/1‑68 (MQ=255)
                          tgtgGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt                      <  1:3303790/69‑1 (MQ=255)
                           gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                       >  1:912962/1‑67 (MQ=255)
                           gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGgtggtg                     >  1:2115442/1‑69 (MQ=255)
                                   tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga             >  1:2077280/1‑69 (MQ=255)
                                   tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga             >  1:3089497/1‑69 (MQ=255)
                                   tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga             >  1:2184702/1‑69 (MQ=255)
                                   tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga             >  1:2710180/1‑69 (MQ=255)
                                   tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga             >  1:1289280/1‑69 (MQ=255)
                                     aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGga                >  1:1180736/1‑64 (MQ=255)
                                     aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag           >  1:2089168/1‑69 (MQ=255)
                                     aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag           >  1:540956/1‑69 (MQ=255)
                                     aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag           >  1:2721466/1‑69 (MQ=255)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:554411/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:357796/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:3020342/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:1790361/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:2956298/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:1838338/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:3186756/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:2177427/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:1326801/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:1647915/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:3197161/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:2061916/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:1401308/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:745873/1‑37 (MQ=25)
                                        aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc                                        >  1:2468780/1‑37 (MQ=25)
                                            gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa    <  1:902955/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGGCGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:3265477/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1327710/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:2049078/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:383387/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1379368/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:747051/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:830459/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:627055/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:77709/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1258623/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:328907/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:2217703/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:2252857/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1782504/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:2375466/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1758086/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:244404/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:172296/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:2601116/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:26156/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:2735877/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:279014/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1633490/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:2976774/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1476252/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:1404523/69‑1 (MQ=255)
                                              cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:3397796/69‑1 (MQ=255)
                                               gTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:284224/68‑1 (MQ=255)
                                                      tGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:114505/61‑1 (MQ=255)
                                                       gggtgAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGTAAGAAGTGGCGaaaa  <  1:2198662/60‑1 (MQ=255)
                                                                    |                                              
GCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAA  >  minE/2872201‑2872315

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