Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,733,673 C→T 31.5% P90S (CCC→TCC)  yheL → predicted intracellular sulfur oxidation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,733,6730CT31.5% 105.4 / 56.1 114P90S (CCC→TCC) yheLpredicted intracellular sulfur oxidation protein
Reads supporting (aligned to +/- strand):  ref base C (44/34);  new base T (20/16);  total (64/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

Only 100 of 118 total aligned reads displayed.

TGGTCAAATTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCGC  >  minE/2733606‑2733736
                                                                    |                                                               
tGGTCAAATTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCTGACTTACGGTTAAGCAcc                                                                 <  1:2661421/68‑1 (MQ=255)
tGGTCAAATTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCAcc                                                                 <  1:2186732/68‑1 (MQ=255)
tGGTCAAATTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCAcc                                                                 <  1:2287360/68‑1 (MQ=255)
tGGTCAAATTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCAcc                                                                 <  1:1843147/68‑1 (MQ=255)
  gTCAAATTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCcca                                                              >  1:3337698/1‑69 (MQ=255)
       atatCGAACGACATCATTCTCAATGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:743673/65‑1 (MQ=255)
       aTTTCTAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:572801/68‑1 (MQ=255)
       aTTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:1195793/68‑1 (MQ=255)
       aTTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:3089593/68‑1 (MQ=255)
        tttCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:1750740/67‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:677104/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:1421216/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:641976/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:1328933/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:412983/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:335883/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:1072249/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:740934/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCccagcca                                                          <  1:838110/53‑1 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATgg                                          >  1:931649/1‑69 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATgg                                          >  1:2309196/1‑69 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATgg                                          >  1:2916439/1‑69 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATgg                                          >  1:2910773/1‑69 (MQ=255)
                      atTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATg                                           >  1:1970172/1‑68 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:2894487/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:3199689/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:906781/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:2563160/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:2512831/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:2509060/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:2819004/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:3348391/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:967662/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:1081599/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:1575001/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:1900832/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGTGATGGc                                         <  1:149598/65‑1 (MQ=255)
                        tcccATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACTCCAGCCAGATGGCCTGGCGATGGc                                         <  1:2180207/65‑1 (MQ=255)
                         ctcATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCgg                                       >  1:1814551/1‑69 (MQ=255)
                               gACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATcat                                 >  1:500557/1‑67 (MQ=255)
                               gACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGt                                 >  1:2509559/1‑69 (MQ=255)
                               gACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGt                                 >  1:514435/1‑69 (MQ=255)
                               gACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGt                                 >  1:818758/1‑69 (MQ=255)
                                   ataCTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGta                             <  1:1301977/69‑1 (MQ=255)
                                    taCTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGtat                            >  1:504640/1‑69 (MQ=255)
                                    taCTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGtat                            >  1:1853993/1‑69 (MQ=255)
                                     aCTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGta                             <  1:3065899/67‑1 (MQ=255)
                                      cTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGtata                           >  1:225185/1‑68 (MQ=255)
                                      cTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGtata                           >  1:3386104/1‑68 (MQ=255)
                                      cTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGtata                           >  1:693766/1‑68 (MQ=255)
                                      cTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGtata                           >  1:2087296/1‑68 (MQ=255)
                                      cTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTatat                          >  1:3049854/1‑69 (MQ=255)
                                      cTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTatat                          >  1:1671473/1‑69 (MQ=255)
                                      cTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTatat                          >  1:165728/1‑69 (MQ=255)
                                         aTTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTc                       >  1:367709/1‑69 (MQ=255)
                                           ttCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCtt                     <  1:1732909/69‑1 (MQ=255)
                                           ttCGTCAGACTTACGGTTAAGAACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCtt                     <  1:590375/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCTTTGTATATTTCTTGacac                <  1:878939/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:349133/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:3375325/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:3418632/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:2142694/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:1274035/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:1134881/69‑1 (MQ=255)
                                                cAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:2863671/69‑1 (MQ=255)
                                                 aGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:3156651/68‑1 (MQ=255)
                                                 aGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGacac                <  1:2753449/68‑1 (MQ=255)
                                                   aCTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCtt             >  1:3184484/1‑69 (MQ=255)
                                                   aCTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCtt             >  1:24277/1‑69 (MQ=255)
                                                   aCTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCtt             >  1:1472448/1‑69 (MQ=255)
                                                   aCTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCt              >  1:2945250/1‑68 (MQ=255)
                                                    cTTACGGTTAAGCACCCCCAGCCAGATGGCCTGGTGATGGCGGGATCTTTGTATATTTCTTGACACCtt             <  1:3273924/69‑1 (MQ=255)
                                                    cTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGa                   >  1:598991/1‑62 (MQ=255)
                                                    cTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGa                   >  1:2206832/1‑62 (MQ=255)
                                                    cTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGa                   >  1:1221081/1‑62 (MQ=255)
                                                    cTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGa                   >  1:2518216/1‑62 (MQ=255)
                                                    cTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGa                   >  1:425526/1‑62 (MQ=255)
                                                    cTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGa                   >  1:1947033/1‑62 (MQ=255)
                                                            tAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCAt     >  1:2703901/1‑68 (MQ=255)
                                                            tAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATc    >  1:376687/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTatat                          >  1:2100476/1‑46 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCtt                     >  1:2727684/1‑51 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCt                      >  1:771512/1‑50 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCtt             >  1:2688178/1‑59 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:1513959/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:1528803/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:1116579/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:1802713/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:885645/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:516410/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:3099633/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:2304440/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:3142810/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:3149577/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:2337506/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:1407635/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:3233214/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:1602315/1‑69 (MQ=255)
                                                             aaGCACTCCAGCCAGATGGCCTGGTGATGGAGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCg   >  1:727864/1‑69 (MQ=255)
                                                               gCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCGc  <  1:754088/68‑1 (MQ=255)
                                                                    |                                                               
TGGTCAAATTTCGAACGACATCATTCTCATTGACTATACTGATTTCGTCAGACTTACGGTTAAGCACCCCAGCCAGATGGCCTGGTGATGGCGGGATCGTTGTATATTTCTTGACACCTTTTCGGCATCGC  >  minE/2733606‑2733736

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