Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,293,839 C→A 44.8% G129G (GGG→GGT cbl ← DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,293,8390CA44.8% 24.0 / 53.7 58G129G (GGG→GGTcblDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (13/19);  new base A (16/10);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

ACGCTCGCTGGCGATACCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTC  >  minE/1293773‑1293898
                                                                  |                                                           
aCGCTCGCTGGCGATACCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTc                                                             <  1:1259942/67‑1 (MQ=255)
aCGCTCGCTGGCGATACCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTc                                                             <  1:3339827/67‑1 (MQ=255)
               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:2216290/69‑1 (MQ=255)
                   atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                        <  1:1956732/69‑1 (MQ=255)
                   atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                        <  1:1881122/69‑1 (MQ=255)
                   atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                        <  1:3382019/69‑1 (MQ=255)
                   atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                        <  1:1417802/69‑1 (MQ=255)
                   atatCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGa                                        <  1:2461205/69‑1 (MQ=255)
                        aGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTgg                                                       <  1:2517724/49‑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:320809/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:3243634/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:3299173/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:331489/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:3034404/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:2716171/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:467437/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:712178/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:738839/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:868437/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:1675227/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:1762984/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:1823508/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:2411853/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:2190662/1‑63 (MQ=255)
                                   tttGCAACAACGTCGCAATTTCCTGTGGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGa                              >  1:2226343/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:2302878/69‑1 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    >  1:1977045/1‑69 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    <  1:1042751/69‑1 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    >  1:2432338/1‑69 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    >  1:1261594/1‑69 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    <  1:2677240/69‑1 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    >  1:1734634/1‑69 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    >  1:1493766/1‑69 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    >  1:14402/1‑69 (MQ=255)
                                                       tCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGAcc    <  1:353439/69‑1 (MQ=255)
                                                          tgtgGTGTACCCTTGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTATTGACctc  <  1:1035155/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:478415/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:1336575/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:1328644/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:1387458/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:3270729/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:785902/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:2436497/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCAGCCGAAACTCCGGGAAAAGTTCGCTAAAAGCTTTAATGACctc  <  1:1795364/68‑1 (MQ=255)
                                                          tgtgGTGTACCCTGGATTAGCTCCACCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACctc  <  1:2842121/68‑1 (MQ=255)
                                                                  |                                                           
ACGCTCGCTGGCGATACCAATATCAGCTTCGCCATTTTGCAACAATGTCGCAATTTCCTGTGGCGTCCCCTGGATTAGCTCGAGCCGAACCTCCGGGAAAAGTTCGCGAAAAGCTTTAATGACCTC  >  minE/1293773‑1293898

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