Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,662,853 T→G 6.3% intergenic (‑25/‑164) 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,8530TG 6.3% 204.4 / 5.0 78intergenic (‑25/‑164)rtcR/rtcBsigma 54‑dependent transcriptional regulator of rtcBA expression/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/40);  new base G (2/3);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTT  >  minE/2662784‑2662921
                                                                     |                                                                    
aTCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGata                                                                     >  1:59063/1‑71 (MQ=255)
aTCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGata                                                                     >  1:2069659/1‑71 (MQ=255)
aTCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGata                                                                     >  1:2703249/1‑71 (MQ=255)
   cAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATAttt                                                                  >  1:1903551/1‑71 (MQ=255)
   cAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATAttt                                                                  <  1:2779685/71‑1 (MQ=255)
       aCGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAt                                                               >  1:942755/1‑70 (MQ=255)
        cGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAtt                                                              <  1:1676593/70‑1 (MQ=255)
        cGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAtt                                                              <  1:2782211/70‑1 (MQ=255)
        cGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTa                                                             <  1:2124237/71‑1 (MQ=255)
        cGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTa                                                             >  1:3656181/1‑71 (MQ=255)
              cGACAAAGCCAAAAGCGACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTa                                                       <  1:1705867/71‑1 (MQ=255)
               gACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTc                                                                 >  1:2226927/1‑60 (MQ=255)
                aCAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGt                                                        <  1:1146851/68‑1 (MQ=255)
                  aaagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATAttt                                                                  <  1:3429813/56‑1 (MQ=255)
                  aaagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATAttt                                                                  <  1:5025205/56‑1 (MQ=255)
                  aaagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTc                                                                 >  1:1368732/1‑57 (MQ=255)
                   aagccaaaagccaCTGTTTTACTCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                  <  1:3111058/71‑1 (MQ=255)
                   aagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                  <  1:3192046/71‑1 (MQ=255)
                   aagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                  >  1:4572051/1‑71 (MQ=255)
                    agccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTgg                                                 >  1:1486444/1‑71 (MQ=255)
                       caaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCtt                                              >  1:3703833/1‑71 (MQ=255)
                       caaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCtt                                              >  1:1860532/1‑71 (MQ=255)
                         aaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAtt                                                              <  1:4661556/53‑1 (MQ=255)
                         aaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTAt                                                            >  1:3494873/1‑55 (MQ=255)
                         aaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTg                                            <  1:5627366/71‑1 (MQ=255)
                          aagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGc                                           <  1:4933504/71‑1 (MQ=255)
                           agccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCg                                          <  1:2357371/71‑1 (MQ=255)
                           agccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCg                                          <  1:4324270/71‑1 (MQ=255)
                           agccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCg                                          <  1:3614063/71‑1 (MQ=255)
                            gccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGAGATTTCATTATCGGTAAAtcg                                                  <  1:4545447/62‑3 (MQ=255)
                              caCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCggg                                        <  1:4702398/70‑1 (MQ=255)
                                 tGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAt                                    >  1:1734513/1‑71 (MQ=255)
                                 tGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAt                                    >  1:1478383/1‑71 (MQ=255)
                                  gTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAt                                                               >  1:1111529/1‑43 (MQ=255)
                                  gTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAt                                                               >  1:1230870/1‑43 (MQ=255)
                                  gTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGc                                                <  1:3583512/58‑1 (MQ=255)
                                  gTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGc                                                <  1:4968941/58‑1 (MQ=255)
                                   ttttACGCATCTTAGATATCCTTATAAAAGACGATAttt                                                                  <  1:3705698/39‑1 (MQ=255)
                                   ttttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtc                                   <  1:3434125/70‑1 (MQ=255)
                                    tttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtct                                  >  1:3723354/1‑70 (MQ=255)
                                     ttACGCATCTTAGATATCCTTATAAAAGACGATATTTCa                                                                >  1:351370/1‑39 (MQ=255)
                                     ttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAAt                                                    >  1:3513199/1‑51 (MQ=255)
                                     ttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGaaa                                     >  1:5109758/1‑66 (MQ=255)
                                     ttACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAAtctc                                 <  1:4391898/70‑1 (MQ=255)
                                      tACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTcc                                <  1:5048090/70‑1 (MQ=255)
                                      tACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTcc                                <  1:2869614/70‑1 (MQ=255)
                                      tACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTcc                                <  1:1698806/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)
                                            tCTTAGATATCCTTATAAAAGACGATATTTCATTAt                                                            <  1:3024991/36‑1 (MQ=255)
                                            tCTTAGATATCCTTATAAAAGACGATATTTCATTAt                                                            <  1:4468802/36‑1 (MQ=255)
                                            tCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTa                                                       <  1:590022/41‑1 (MQ=255)
                                            tCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTgg                                                 <  1:1641511/47‑1 (MQ=255)
                                            tCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCt                                               >  1:4344058/1‑49 (MQ=255)
                                            tCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCt                                               >  1:2250792/1‑49 (MQ=255)
                                             cTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTa                                                       <  1:229228/40‑1 (MQ=255)
                                             cTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTa                                                       <  1:5693997/40‑1 (MQ=255)
                                              ttAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAAt                                                    >  1:2868786/1‑42 (MQ=255)
                                              ttAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAAt                                                    >  1:5324664/1‑42 (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:5165286/1‑53 (MQ=255)
                                                              aaGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACg                         >  1:2608953/1‑53 (MQ=255)
                                                              aaGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACg                         >  1:3197803/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)
                                                                     |                                                                    
ATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTT  >  minE/2662784‑2662921

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: