Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 434,040 C→T 20.4% intergenic (‑93/+196) fur ← / ← fldA DNA‑binding transcriptional dual regulator/flavodoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE434,0400CT20.4% 111.8 / 22.5 54intergenic (‑93/+196)fur/fldADNA‑binding transcriptional dual regulator/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base C (27/16);  new base T (4/7);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

ACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAA  >  minE/433977‑434109
                                                               |                                                                     
aCTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCATGCCTaa                                                                 <  1:1842672/70‑1 (MQ=255)
 cTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAAt                                                               >  1:826344/1‑71 (MQ=255)
  ttCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCAc                                                                       <  1:744299/62‑1 (MQ=255)
  ttCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAAt                                                               >  1:3193989/1‑70 (MQ=255)
    cATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGc                                                                     >  1:3670245/1‑62 (MQ=255)
    cATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGa                                                             >  1:2127251/1‑70 (MQ=255)
     aTTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGaa                                                            <  1:1803114/70‑1 (MQ=255)
         gaagaaGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAAcc                                                        <  1:1930052/70‑1 (MQ=255)
            gaagTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACaa                                                    <  1:3056422/71‑1 (MQ=255)
             aagTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAg                                                   >  1:1645647/1‑71 (MQ=255)
               gTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCAc                                                                       >  1:3155798/1‑49 (MQ=255)
                tGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTcc                                                <  1:171630/71‑1 (MQ=255)
                  aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTccc                                               <  1:2589069/70‑1 (MQ=255)
                  aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTccc                                               <  1:2589974/70‑1 (MQ=255)
                  aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTCCCt                                              >  1:3874686/1‑71 (MQ=255)
                  aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTCACt                                              >  1:2988077/1‑71 (MQ=255)
                  aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTaa                                                                 <  1:2801211/52‑1 (MQ=255)
                  aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTccc                                               <  1:2497498/70‑1 (MQ=255)
                  aCAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCt                                              >  1:3276945/1‑71 (MQ=255)
                    aGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTgg                                            >  1:1452451/1‑71 (MQ=255)
                    aGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTgg                                            >  1:3668141/1‑71 (MQ=255)
                     gAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTcc                                                >  1:2449317/1‑66 (MQ=255)
                      aaTTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGc                                           >  1:2284710/1‑70 (MQ=255)
                      aaTTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGc                                           >  1:1973637/1‑70 (MQ=255)
                      aaTTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCa                                          <  1:21083/71‑1 (MQ=255)
                      aaTTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCa                                          <  1:3839987/71‑1 (MQ=255)
                               ttGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGt                                  >  1:2567255/1‑70 (MQ=255)
                                tGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGt                                  >  1:3183985/1‑69 (MQ=255)
                                 gagaTAATGCGTATCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTCCCTGGCAAATAtt                                    >  1:2546663/1‑66 (MQ=255)
                                 gagaTAATGCGTATCATTATAGAATTGCCACGCCTAAATGaaaa                                                          >  1:3809828/1‑44 (MQ=255)
                                   gaTAATGAGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAatat                                     >  1:1235505/1‑63 (MQ=255)
                                    aTAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaataa                            <  1:1758173/71‑1 (MQ=255)
                                       aTGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCa                                          >  1:2683460/1‑54 (MQ=255)
                                       aTGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATaaaa                          <  1:3292223/70‑1 (MQ=255)
                                        tGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAAc                        >  1:170418/1‑71 (MQ=255)
                                         gCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATAtt                                    >  1:282383/1‑58 (MQ=255)
                                          cGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACg                       <  1:2969514/70‑1 (MQ=255)
                                            tATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGt                                                  <  1:1844917/41‑1 (MQ=255)
                                             aTCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCACt                                              >  1:1365728/1‑44 (MQ=255)
                                              tCATTATAGAATTGCCATGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTgaat                  >  1:184715/1‑68 (MQ=255)
                                               cATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTtaataa                            <  1:2387722/60‑1 (MQ=255)
                                                aTTATAGAATTGCCATGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGcac                 <  1:1912007/70‑1 (MQ=255)
                                                aTTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGca                  >  1:1744752/1‑69 (MQ=255)
                                                aTTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGca                  >  1:3583159/1‑69 (MQ=255)
                                                 ttATAGAATTGCCATGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAg               <  1:1493181/71‑1 (MQ=255)
                                                 ttATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAg               >  1:2431517/1‑71 (MQ=255)
                                                 ttATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAg               >  1:2322199/1‑71 (MQ=255)
                                                 ttATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAg               >  1:3471661/1‑71 (MQ=255)
                                                  tataGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGGACAg               >  1:3713505/1‑70 (MQ=255)
                                                    taGAATTGCCATGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCt             <  1:947789/70‑1 (MQ=255)
                                                    taGAATTGCCATGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCt             <  1:952323/70‑1 (MQ=255)
                                                    taGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTta                                >  1:3396593/1‑51 (MQ=255)
                                                      gAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTa           >  1:2664680/1‑70 (MQ=255)
                                                      gAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTAc          <  1:3350299/71‑1 (MQ=255)
                                                         ttGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAatat                                     <  1:2390854/41‑1 (MQ=255)
                                                              aCGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACaa  <  1:1503887/71‑1 (MQ=255)
                                                               |                                                                     
ACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAATTGCCACGCCTAAATGAAAACCACAAGTCCCTGGCAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAA  >  minE/433977‑434109

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: