Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,662,872 T→C 6.0% intergenic (‑44/‑145) rtcR ← / → rtcB sigma 54‑dependent transcriptional regulator of rtcBA expression/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,662,8720TC 6.0% 189.6 / 2.2 66intergenic (‑44/‑145)rtcR/rtcBsigma 54‑dependent transcriptional regulator of rtcBA expression/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base C (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

AAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAA  >  minE/2662803‑2662938
                                                                     |                                                                  
aagccaaaagccaCTGTTTTACTCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                                   <  1:3111058/71‑1 (MQ=255)
aagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                                   >  1:4572051/1‑71 (MQ=255)
aagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                                   <  1:3192046/71‑1 (MQ=255)
 agccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTgg                                                                  >  1:1486444/1‑71 (MQ=255)
    caaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCtt                                                               >  1:3703833/1‑71 (MQ=255)
    caaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCtt                                                               >  1:1860532/1‑71 (MQ=255)
      aaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTg                                                             <  1:5627366/71‑1 (MQ=255)
       aagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGc                                                            <  1:4933504/71‑1 (MQ=255)
        agccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCg                                                           <  1:3614063/71‑1 (MQ=255)
        agccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCg                                                           <  1:2357371/71‑1 (MQ=255)
        agccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCg                                                           <  1:4324270/71‑1 (MQ=255)
           caCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCggg                                                         <  1:4702398/70‑1 (MQ=255)
              tGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAt                                                     >  1:1734513/1‑71 (MQ=255)
              tGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAt                                                     >  1:1478383/1‑71 (MQ=255)
               gTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGc                                                                 <  1:3583512/58‑1 (MQ=255)
               gTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGc                                                                 <  1:4968941/58‑1 (MQ=255)
                ttttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtc                                                    <  1:3434125/70‑1 (MQ=255)
                 tttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtct                                                   >  1:3723354/1‑70 (MQ=255)
                  ttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGaaa                                                      >  1:5109758/1‑66 (MQ=255)
                  ttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtctc                                                  <  1:4391898/70‑1 (MQ=255)
                   tACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTcc                                                 <  1:1698806/70‑1 (MQ=255)
                   tACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTcc                                                 <  1:2869614/70‑1 (MQ=255)
                   tACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTcc                                                 <  1:5048090/70‑1 (MQ=255)
                   tACGCATCTTAGATATCCTTATAAAAGACGAGATTTCATTATCGTTAAATCGGCTTTGCGGGAAATCTcc                                                 <  1:689451/70‑1 (MQ=255)
                   tACGCATCTTAGATATCCTTATAAAAGACGAGATTTCATTATCGGTAAATCGGCTTTGCGGGAAATCTcc                                                 <  1:4562811/70‑1 (MQ=255)
                    aCGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCg                                                >  1:1312398/1‑70 (MQ=255)
                      gCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGtt                                             <  1:4430214/71‑1 (MQ=255)
                         tCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTgg                                                                  <  1:1641511/47‑1 (MQ=255)
                         tCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCt                                                                >  1:2250792/1‑49 (MQ=255)
                         tCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCt                                                                >  1:4344058/1‑49 (MQ=255)
                           ttAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCg                                                           >  1:1521712/1‑52 (MQ=255)
                           ttAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTAc                                           >  1:4710040/1‑68 (MQ=255)
                            tAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGAc                                       >  1:142300/1‑71 (MQ=255)
                             aGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACgg                                         <  1:634806/68‑1 (MQ=255)
                              gATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACg                                      <  1:1063710/70‑1 (MQ=255)
                               atatCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTa                                    <  1:1375955/71‑1 (MQ=255)
                                   ccTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGaaa                                                      >  1:729746/1‑49 (MQ=255)
                                    cTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGa                                >  1:1647707/1‑70 (MQ=255)
                                        taAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCggg                                                         <  1:4779068/41‑1 (MQ=255)
                                           aaGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACg                                          >  1:3197803/1‑53 (MQ=255)
                                           aaGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACg                                          >  1:2608953/1‑53 (MQ=255)
                                           aaGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACg                                          >  1:5165286/1‑53 (MQ=255)
                                           aaGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCt                        >  1:4464027/1‑71 (MQ=255)
                                            aGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCtt                       <  1:137904/71‑1 (MQ=255)
                                              aCGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtctc                                                  <  1:2744316/42‑1 (MQ=255)
                                              aCGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtctc                                                  <  1:3631012/42‑1 (MQ=255)
                                              aCGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTaa                                   >  1:1277313/1‑57 (MQ=255)
                                               cGATATTTCATTATCTGTAAATTGGCTTTGCGGGa                                                        <  1:3941655/35‑1 (MQ=255)
                                               cGAGATTTCATTATCGGTAAATCGGCTTTGCGGGAAATCTCCGTTTACGGACGTAAAAAtata                            >  1:5322624/1‑63 (MQ=255)
                                               cGAGATTTCATTATCGGTAAATCGGCTTTGCGGGAAATCTCCGTTTACGGACGTAAAAAtata                            >  1:2565239/1‑63 (MQ=255)
                                                 atatTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCtttttt                   <  1:600098/70‑1 (MQ=255)
                                                   atTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCttt                      <  1:2161843/65‑1 (MQ=255)
                                                       cATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGa                                        >  1:1591627/1‑43 (MQ=255)
                                                           aTCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGaa        >  1:1521795/1‑71 (MQ=255)
                                                           aTCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGaa        >  1:2231509/1‑71 (MQ=255)
                                                           aTCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGa         <  1:729549/70‑1 (MQ=255)
                                                              ggTAAATTGGCTTTGCGGGAAGTCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTa     <  1:2292976/71‑1 (MQ=255)
                                                              ggTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTa     <  1:2616705/71‑1 (MQ=255)
                                                              ggTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTa     <  1:3118744/71‑1 (MQ=255)
                                                                 aaaTTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTaaa   >  1:3281735/1‑70 (MQ=255)
                                                                  aaTTGGCTTTGCGGGAAATCTCCGGTTACGGACGTaaa                                  <  1:356332/38‑1 (MQ=255)
                                                                  aaTTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCtt                       >  1:2612550/1‑49 (MQ=255)
                                                                  aaTTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTaaaa  <  1:2538305/70‑1 (MQ=255)
                                                                  aaTTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTaaaa  <  1:4439980/70‑1 (MQ=255)
                                                                  aaTTGGCTTTGCCGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAAc              <  1:3526366/58‑1 (MQ=255)
                                                                   aTTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTaaaa               >  1:787774/1‑56 (MQ=255)
                                                                     tGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCt                        >  1:4225044/1‑45 (MQ=255)
                                                                     |                                                                  
AAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAA  >  minE/2662803‑2662938

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: