Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 171,888 A→T 32.0% T64T (ACA→ACT frr → ribosome recycling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE171,8880AT32.0% 106.3 / 72.5 100T64T (ACA→ACTfrrribosome recycling factor
Reads supporting (aligned to +/- strand):  ref base A (21/47);  new base T (20/12);  total (41/59)
Fisher's exact test for biased strand distribution p-value = 4.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TTACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGATT  >  minE/171828‑171942
                                                            |                                                      
ttACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGaaaa                                                 >  1:1093497/1‑68 (MQ=255)
ttACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGaaaa                                                 >  1:3272931/1‑68 (MQ=255)
ttACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGaaaa                                                 >  1:243314/1‑68 (MQ=255)
ttACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGaaaa                                                 >  1:1382587/1‑68 (MQ=255)
ttACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGaaaa                                                 >  1:3140011/1‑68 (MQ=255)
ttACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGaaaa                                                 >  1:762519/1‑68 (MQ=255)
     gCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAAc                                            >  1:3410631/1‑68 (MQ=255)
          ccgaggccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCa                                              >  1:158550/1‑61 (MQ=255)
          ccgacgccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGttt                                      >  1:2952908/1‑69 (MQ=255)
          ccgacgccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGttt                                      >  1:1969909/1‑69 (MQ=255)
          ccgacgccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGttt                                      >  1:1700307/1‑69 (MQ=255)
          ccgacgccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGttt                                      >  1:1682760/1‑69 (MQ=255)
          ccgacgccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGttt                                      >  1:1678727/1‑69 (MQ=255)
          ccgacgccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGttt                                      >  1:1553278/1‑69 (MQ=255)
          ccgacgccgCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGttt                                      >  1:1201523/1‑69 (MQ=255)
                  gcTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTc                              <  1:612114/69‑1 (MQ=255)
                  gcTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTc                              <  1:769489/69‑1 (MQ=255)
                        tCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGtc                        <  1:68391/69‑1 (MQ=255)
                        tCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGtc                        <  1:1761290/69‑1 (MQ=255)
                        tCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGtc                        <  1:1770197/69‑1 (MQ=255)
                        tCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGtc                        <  1:3320011/69‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2153490/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2653435/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:3230568/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2137041/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:3227686/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2298099/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:236037/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2570213/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2920486/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2438953/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2441894/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:1188569/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:674603/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:392360/39‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCg                <  1:2623778/68‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCg                <  1:3380795/68‑1 (MQ=255)
                                 aaGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCg                <  1:1544346/68‑1 (MQ=255)
                                   gCGTAACGGTAGAAGATTCCCGTACACTGAAAATCaa                                             <  1:2348712/37‑1 (MQ=255)
                                   gCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTc                              <  1:9682/52‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:2075344/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:1269415/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:3116866/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:2483857/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:3165537/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:450298/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:3358045/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:3169482/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:1952265/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:1974482/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:3262090/69‑1 (MQ=255)
                                             agaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:210265/69‑1 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTg             >  1:1097795/1‑58 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:415152/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:2062015/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:2171226/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:3228404/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:3277526/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:3187082/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:1951497/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:3387473/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:3411668/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:1540033/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:302623/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:576172/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:1434331/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:2548950/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:2940480/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:265993/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:828447/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGa            >  1:984762/1‑59 (MQ=255)
                                              gaagaTTCCCGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTATCCGGCCGTTg             >  1:2721586/1‑58 (MQ=39)
                                              gaagaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAt   <  1:2738073/68‑1 (MQ=255)
                                                agaTTCCCTTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCg     >  1:2519894/1‑64 (MQ=255)
                                                agaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCg                >  1:649397/1‑53 (MQ=255)
                                                agaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCg     >  1:1922638/1‑64 (MQ=255)
                                                agaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCg     >  1:2469511/1‑64 (MQ=255)
                                                agaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCg     >  1:2498470/1‑64 (MQ=255)
                                                agaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCg     >  1:63889/1‑64 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:2283617/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:2863827/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:1172251/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:1222603/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:1324153/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:598969/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:485452/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:1501044/41‑1 (MQ=255)
                                                 gaTTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAAt                           <  1:2292530/41‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:525355/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:567669/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:2964684/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:2963358/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:648579/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:3078966/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:1731021/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:1296495/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:1289204/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:191127/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:2891540/61‑1 (MQ=255)
                                                      ccGTACTCTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGAtt  <  1:1978938/61‑1 (MQ=255)
                                                            |                                                      
TTACGGCACGCCGACGCCGCTGCGTCAGCTGGCAAGCGTAACGGTAGAAGATTCCCGTACACTGAAAATCAACGTGTTTGATCGTTCAATGTCTCCGGCCGTTGAAAAAGCGATT  >  minE/171828‑171942

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: