Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,699,098 C→T 36.4% F326F (TTC→TTT csiE → stationary phase inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,699,0980CT36.4% 39.3 / 36.3 63F326F (TTC→TTTcsiEstationary phase inducible protein
Reads supporting (aligned to +/- strand):  ref base C (26/14);  new base T (11/12);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

GGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGG  >  minE/1699033‑1699158
                                                                 |                                                            
ggATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGt                                                           <  1:516141/69‑1 (MQ=255)
     tGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                                      <  1:930618/69‑1 (MQ=255)
     tGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                                      <  1:3224477/69‑1 (MQ=255)
     tGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                                      <  1:573378/69‑1 (MQ=255)
      gAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                                      <  1:3086634/68‑1 (MQ=255)
      gAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                                      <  1:3106256/68‑1 (MQ=255)
         aGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCt                                                   <  1:1686863/68‑1 (MQ=255)
                aCGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAg                                            >  1:1338400/1‑68 (MQ=255)
                aCGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAg                                            >  1:3007300/1‑68 (MQ=255)
                 cGGCGTCCATCTTTCTGATGAGGAAGGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAgg                                           >  1:1474800/1‑68 (MQ=255)
                 cGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATTCAgg                                           >  1:122002/1‑68 (MQ=255)
                 cGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAgg                                           >  1:1236708/1‑68 (MQ=255)
                 cGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAgg                                           >  1:1995739/1‑68 (MQ=255)
                   gCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGaa                                         <  1:2726659/68‑1 (MQ=255)
                      tCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAAc                                      >  1:840798/1‑68 (MQ=255)
                      tCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACg                                     >  1:1343292/1‑69 (MQ=255)
                        cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                                   >  1:445025/1‑69 (MQ=255)
                        cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                                   >  1:3333019/1‑69 (MQ=255)
                        cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                                   >  1:288099/1‑69 (MQ=255)
                        cATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                                   >  1:2185535/1‑69 (MQ=255)
                          tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGTTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                                  >  1:2409344/1‑68 (MQ=255)
                          tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGAt                                   >  1:802418/1‑67 (MQ=255)
                          tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                                  >  1:1658932/1‑68 (MQ=255)
                          tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                                  >  1:307662/1‑68 (MQ=255)
                          tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                                  >  1:1485925/1‑68 (MQ=255)
                          tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                                  >  1:2772686/1‑68 (MQ=255)
                          tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATc                                  >  1:1888167/1‑68 (MQ=255)
                           cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                                >  1:2100424/1‑69 (MQ=255)
                           cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                                >  1:1522807/1‑69 (MQ=255)
                           cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                                >  1:1893693/1‑69 (MQ=255)
                           cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                                >  1:2236664/1‑69 (MQ=255)
                           cTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTg                                >  1:2365187/1‑69 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:2849114/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:2347926/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:2903541/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:3040980/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:316880/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:774549/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:2547270/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:2519207/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:1879509/1‑59 (MQ=255)
                            tttCTGATGAGGAAAGCGGTCTGGTCGCGGTGAATTTTGGt                                                           >  1:1559450/1‑41 (MQ=37)
                            tttCTGATGAGCAAAGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGaa                                         >  1:2086869/1‑59 (MQ=255)
                                   tgatgaAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTg                                                      <  1:1637531/39‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:948341/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:2631683/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:1523948/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:1581362/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:1700660/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:1718305/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:1870099/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:197341/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:2140985/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:2158215/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:2966101/69‑1 (MQ=255)
                                         aaGCGGTCTGGTCGCGGTGATTTTTGGTGCCTGGCTAATGCAGGAAAACGACCTGCATGAGAAACAGat                  <  1:2472251/69‑1 (MQ=255)
                                                         gTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCggg  <  1:3154899/69‑1 (MQ=255)
                                                         gTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCggg  <  1:178383/69‑1 (MQ=255)
                                                         gTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCggg  <  1:2205071/69‑1 (MQ=255)
                                                         gTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCggg  <  1:2223506/69‑1 (MQ=255)
                                                         gTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCggg  <  1:1170872/69‑1 (MQ=255)
                                                           gATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGaaa                         >  1:692972/1‑44 (MQ=255)
                                                           gATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGaaa                         >  1:2475609/1‑44 (MQ=255)
                                                                 |                                                            
GGATTTGAAAGTGAATACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGG  >  minE/1699033‑1699158

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: