Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,733,426 C→T 32.4% R7R (CGC→CGT yheL → predicted intracellular sulfur oxidation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,733,4260CT32.4% 114.6 / 66.4 129R7R (CGC→CGTyheLpredicted intracellular sulfur oxidation protein
Reads supporting (aligned to +/- strand):  ref base C (53/34);  new base T (13/29);  total (66/63)
Fisher's exact test for biased strand distribution p-value = 2.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

ACTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCTAT  >  minE/2733365‑2733490
                                                             |                                                                
aCTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCt                                                           <  1:1957137/69‑1 (MQ=255)
aCTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCt                                                           <  1:2122311/69‑1 (MQ=255)
aCTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCt                                                           <  1:2564181/69‑1 (MQ=255)
 cTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCTg                                                          >  1:1032150/1‑69 (MQ=255)
 cTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCTg                                                          >  1:1863228/1‑69 (MQ=255)
 cTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCTg                                                          >  1:3235872/1‑69 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCt                                                                 >  1:851013/1‑36 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                >  1:326153/1‑69 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                >  1:2223399/1‑69 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                >  1:224470/1‑69 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                >  1:2519492/1‑69 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                <  1:160373/69‑1 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                >  1:3353324/1‑69 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                >  1:10345/1‑69 (MQ=255)
                           gttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCtgctgc                                >  1:765698/1‑69 (MQ=255)
                            ttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCg                               >  1:1936286/1‑69 (MQ=255)
                            ttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCg                               >  1:127969/1‑69 (MQ=255)
                            ttttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCg                               >  1:594664/1‑69 (MQ=255)
                              ttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGt                              >  1:1576856/1‑68 (MQ=255)
                              ttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:1516890/1‑69 (MQ=255)
                              ttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:3403547/1‑69 (MQ=255)
                              ttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:2031195/1‑69 (MQ=255)
                              ttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:61814/1‑69 (MQ=255)
                              ttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:1194723/1‑69 (MQ=255)
                              ttgaggCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGAGTc                             >  1:897982/1‑69 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCg                               >  1:1388604/1‑67 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGt                              >  1:3090086/1‑68 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:575583/1‑69 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:2752169/1‑69 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:2825165/1‑69 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:1030487/1‑69 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:1640426/1‑69 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:1558352/1‑69 (MQ=255)
                              ttgaggCATTTATGCTGCACACATTACATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTc                             >  1:1069025/1‑69 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACgg                                                  >  1:1541546/1‑41 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 <  1:2277538/42‑1 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 <  1:25709/42‑1 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 <  1:1576175/42‑1 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 >  1:3179773/1‑42 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 <  1:2397883/42‑1 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 <  1:2282037/42‑1 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 >  1:612547/1‑42 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 >  1:1097931/1‑42 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 >  1:765773/1‑42 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 >  1:775647/1‑42 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 <  1:1022477/42‑1 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 <  1:1874921/42‑1 (MQ=255)
                                     gTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGa                                                 >  1:2014586/1‑42 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACgg                                                  >  1:2081476/1‑38 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:601250/69‑1 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   >  1:928163/1‑69 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:1545665/69‑1 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   >  1:3027440/1‑69 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   >  1:3002479/1‑69 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:902650/69‑1 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   >  1:2672917/1‑69 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2217876/69‑1 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2347050/69‑1 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2408634/69‑1 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   >  1:1694204/1‑69 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:1848708/69‑1 (MQ=255)
                                        tATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   >  1:1683689/1‑69 (MQ=255)
                                                 cacaTTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGa          <  1:918909/69‑1 (MQ=255)
                                                     ttACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGa          >  1:1945938/1‑65 (MQ=255)
                                                     ttACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGa          >  1:3322393/1‑65 (MQ=255)
                                                     ttACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGa          >  1:763296/1‑65 (MQ=255)
                                                       aCATCGTTCACCCTGGCTTACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2850373/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:740272/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:629094/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:108255/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2257152/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:1246221/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:221791/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:1264219/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:925672/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:366052/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:1302701/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:3374362/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:3325552/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2345082/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:1485016/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:3146053/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:3144855/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:3123897/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:310859/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2753626/54‑1 (MQ=255)
                                                       aCATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2843561/54‑1 (MQ=255)
                                                       aCATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCt    >  1:2495211/1‑69 (MQ=255)
                                                       aCATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCt    >  1:2017410/1‑69 (MQ=255)
                                                       aCATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCt    >  1:1852837/1‑69 (MQ=255)
                                                       aCATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCt    >  1:1661579/1‑69 (MQ=255)
                                                        cATCGTTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGa                   <  1:2883567/53‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:2652893/68‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:2783455/68‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:2159631/68‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:295423/68‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:2532280/68‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:693724/68‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:670640/68‑1 (MQ=255)
                                                          tCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCtat  <  1:575657/68‑1 (MQ=255)
                                                             |                                                                
ACTCGCCAACTACGATGTTATTTTGAGGTTTTGAGGCGTTTATGCTGCACACATTACATCGCTCACCCTGGCTGACGGATTTTGCTGCGCTGCTGCGTCTGCTCAGTGAAGGAGACGAACTGCTAT  >  minE/2733365‑2733490

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