Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 749,866 G→A 30.7% R1118R (CGC→CGT mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE749,8660GA30.7% 76.7 / 37.0 81R1118R (CGC→CGTmfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (4/52);  new base A (15/10);  total (19/62)
Fisher's exact test for biased strand distribution p-value = 8.31e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGC  >  minE/749798‑749928
                                                                    |                                                              
tCACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAgg                                                                >  1:1051051/1‑69 (MQ=255)
 cACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGgc                                                               >  1:2407816/1‑69 (MQ=255)
 cACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGgc                                                               >  1:2515792/1‑69 (MQ=255)
 cACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGgc                                                               >  1:1197205/1‑69 (MQ=255)
 cACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGgc                                                               >  1:571424/1‑69 (MQ=255)
 cACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGgc                                                               >  1:1473935/1‑69 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTTAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:194710/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:3304977/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:3435220/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:3183896/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:3176887/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:3090157/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:3034743/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2966624/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2910566/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2908805/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2861061/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2435807/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:3347698/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2757822/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:420039/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:420428/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:509396/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:527751/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:541642/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:547761/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:604591/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:781130/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:864999/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:945834/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:965540/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:985911/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2339520/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2412659/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:1637919/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:1096422/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:1271416/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:151972/69‑1 (MQ=255)
                      aTTCGATACGCGTTTTCCGCTCACTCAAATCATGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGt                                          <  1:2990381/69‑1 (MQ=255)
                         cGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAt                                       <  1:1506793/69‑1 (MQ=255)
                         cGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAt                                       <  1:364191/69‑1 (MQ=255)
                         cGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAt                                       <  1:360225/69‑1 (MQ=255)
                         cGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAt                                       <  1:355267/69‑1 (MQ=255)
                         cGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAt                                       >  1:2329359/1‑69 (MQ=255)
                         cGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAt                                       <  1:2476876/69‑1 (MQ=255)
                         cGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAt                                       >  1:2155669/1‑69 (MQ=255)
                            tACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAtgc                                     <  1:2595058/68‑1 (MQ=255)
                            tACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAAtgc                                     <  1:2006992/68‑1 (MQ=255)
                                          tcactcaAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTgc                      <  1:1335019/69‑1 (MQ=255)
                                          tcactcaAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTgc                      <  1:3418802/69‑1 (MQ=255)
                                          tcactcaAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTgc                      <  1:1089827/69‑1 (MQ=255)
                                          tcactcaAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTgc                      <  1:2643118/69‑1 (MQ=255)
                                           cactcaAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTgc                      <  1:1960012/68‑1 (MQ=255)
                                               caAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCaa                 >  1:1277162/1‑69 (MQ=255)
                                               caAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCaa                 >  1:1178388/1‑69 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:2812642/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:250938/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:1313893/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:1360579/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:367602/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:2819273/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:1100423/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:2854538/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:1644487/68‑1 (MQ=255)
                                                  aTCCTGAATAAACTTCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAAc               <  1:2491791/68‑1 (MQ=255)
                                                   tCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCa             <  1:899334/69‑1 (MQ=255)
                                                   tCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCa             <  1:2133520/69‑1 (MQ=255)
                                                   tCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCa             <  1:697099/69‑1 (MQ=255)
                                                   tCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCa             <  1:2636353/69‑1 (MQ=255)
                                                   tCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCa             <  1:1199265/69‑1 (MQ=255)
                                                   tCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCa             <  1:2288607/69‑1 (MQ=255)
                                                   tCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCa             <  1:2818197/69‑1 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAgga  >  1:2046548/1‑67 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:723731/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:2083519/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:811353/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:2448287/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:1265731/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:262859/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:1438952/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:406983/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:2778616/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:3379349/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:1533786/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:3293738/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACAATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:1732701/1‑68 (MQ=255)
                                                               ttCAGACGCGTCGGACAATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGc  >  1:306643/1‑68 (MQ=255)
                                                                    |                                                              
TCACGCATAAACTGGCGTACCCATTCGATACGCGTTTTCCGCTCACTCAAATCCTGAATAAATTTCAGGCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGC  >  minE/749798‑749928

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