Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,616,492 G→A 65.3% T38T (ACC→ACT eutJ ← predicted chaperonin, ethanolamine utilization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,616,4920GA65.3% 10.8 / 9.2 23T38T (ACC→ACTeutJpredicted chaperonin, ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base G (6/2);  new base A (2/13);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 6.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCGAGGCATACCGCCACTGGCTGACCGTCGCGGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCGGGCGTCTGGTTACACAG  >  minE/1616431‑1616554
                                                             |                                                              
aGTCGAGGCATACCGCCACCGGCTGACCGTCGCGGTCGACAACCATCGACACCACATCGCAAGTg                                                             >  1:3184739/1‑65 (MQ=255)
aGTCGAGGCATACCGCCACCGGCTGACCGTCGCGGTCGACAACCATCGACACCACATCGCAAGTg                                                             >  1:2480433/1‑65 (MQ=255)
            ccGCCACTGGCTGACCGTCGCGGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCa                                             >  1:692588/1‑69 (MQ=255)
            ccGCCACTGGCTGACCGTCGCGGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCa                                             >  1:1612329/1‑69 (MQ=255)
                       tGACCGTCGCGGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGga                                  <  1:1901026/69‑1 (MQ=255)
                               gcgGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGt                          >  1:1128233/1‑69 (MQ=255)
                               gcgGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGt                          >  1:372783/1‑69 (MQ=255)
                               gcgGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGt                          >  1:1724607/1‑69 (MQ=255)
                               gcgGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGt                          >  1:1119644/1‑69 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:2071188/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:2111584/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:2480323/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:156069/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:3010285/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:3037729/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:3157800/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:1275880/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:3203010/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:3205924/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:1179188/69‑1 (MQ=255)
                                                      caTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:846060/69‑1 (MQ=255)
                                                       aTCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCGGGCGTCTGGTTACACAg  <  1:741766/69‑1 (MQ=255)
                                                       aTCGCAAGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCAGGCGTCTGGTTAcaca   <  1:3382719/68‑1 (MQ=255)
                                                             |                                                              
AGTCGAGGCATACCGCCACTGGCTGACCGTCGCGGTCGACAACCATCGACACCACATCGCAGGTGCCCAGATCCACCCCCAGCCACAGCGGAGATTCCGTCGCGGCGGGCGTCTGGTTACACAG  >  minE/1616431‑1616554

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: