Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,293,854 G→C 34.6% L124L (CTC→CTG cbl ← DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,293,8540GC34.6% 71.5 / 58.3 75L124L (CTC→CTGcblDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (15/34);  new base C (16/10);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 1.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAG  >  minE/1293788‑1293921
                                                                  |                                                                   
acgaaTATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAg                                                                   >  1:720523/4‑69 (MQ=255)
aCCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTTTGGCGTCCCCTGGATTAGCTCGAg                                                                   >  1:2608954/1‑69 (MQ=255)
aCCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCgg                                                                    >  1:3281865/1‑67 (MQ=255)
aCCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAg                                                                   >  1:3037754/1‑69 (MQ=255)
aCCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAg                                                                   >  1:3108470/1‑69 (MQ=255)
aCCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAg                                                                   >  1:3434523/1‑69 (MQ=255)
aCCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAg                                                                   >  1:2258749/1‑69 (MQ=255)
    atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCgt                                                               <  1:3346705/69‑2 (MQ=255)
    atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                                               <  1:3382019/69‑1 (MQ=255)
    atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                                               <  1:2216290/69‑1 (MQ=255)
    atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                                               <  1:1956732/69‑1 (MQ=255)
    atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                                               <  1:1881122/69‑1 (MQ=255)
    atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                                               <  1:1417802/69‑1 (MQ=255)
    atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                                               <  1:2461205/69‑1 (MQ=255)
                 cATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGaaa                                                   <  1:577810/68‑1 (MQ=255)
                  aTTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAg                                                 <  1:2475816/69‑1 (MQ=255)
                  aTTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAg                                                 <  1:435476/69‑1 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:331489/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:2411853/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:868437/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:3299173/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:2226343/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:2716171/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:3243634/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:320809/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:3034404/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:2190662/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:1823508/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:1762984/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:738839/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:1675227/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:712178/1‑63 (MQ=255)
                    tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                                                     >  1:467437/1‑63 (MQ=255)
                      tGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAg                                                 <  1:124799/65‑1 (MQ=255)
                                  gCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCttt                                  <  1:537777/68‑1 (MQ=255)
                                    aaTTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAAt                               <  1:2463215/69‑1 (MQ=255)
                                    aaTTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAAt                               <  1:620338/69‑1 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           >  1:1977045/1‑69 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           <  1:1042751/69‑1 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           >  1:1261594/1‑69 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           <  1:2677240/69‑1 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           >  1:14402/1‑69 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           >  1:1493766/1‑69 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           >  1:2432338/1‑69 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           >  1:1734634/1‑69 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           <  1:2302878/69‑1 (MQ=255)
                                        tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc                           <  1:353439/69‑1 (MQ=255)
                                           tgtgGTGTACCCTTGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTATTGACctc                         <  1:1035155/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:785902/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:3270729/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:1328644/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:2436497/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:1336575/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:1387458/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:478415/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCAGCCGAAACTCCGGGAAAAGTTCGCTAAAAGCTTTAATGACctc                         <  1:1795364/68‑1 (MQ=255)
                                           tgtgGTGTACCCTGGATTAGCTCCACCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc                         <  1:2842121/68‑1 (MQ=255)
                                                         gATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACg          >  1:2888766/1‑69 (MQ=255)
                                                         gATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACg          >  1:663151/1‑69 (MQ=255)
                                                         gATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACg          >  1:1657988/1‑69 (MQ=255)
                                                               cTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCt     <  1:1031476/68‑1 (MQ=255)
                                                               cTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCt     <  1:3361445/68‑1 (MQ=255)
                                                               cTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCt     <  1:2659370/68‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:1818005/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:14584/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:2006351/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:2268265/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:2519329/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:629334/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:638216/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:1277168/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:1272992/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:735574/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:1133479/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:2841694/69‑1 (MQ=255)
                                                                 cGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:857490/69‑1 (MQ=255)
                                                                  gAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAg  <  1:3024130/68‑1 (MQ=255)
                                                                  |                                                                   
ACCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTCTGGCAAGCTATAACGTGCCTGAG  >  minE/1293788‑1293921

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: