Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,681,918 A→G 25.0% R598R (CGT→CGC yhgF ← predicted transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,681,9180AG25.0% 119.1 / 45.2 68R598R (CGT→CGCyhgFpredicted transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base A (31/20);  new base G (10/7);  total (41/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CTTTGATGATGTCGGTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGAATGCGTT  >  minE/2681851‑2681982
                                                                   |                                                                
ctTTGATGATGTCGGTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACg                                                                >  1:948753/1‑70 (MQ=255)
ctTTGATGATGTCGGTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACg                                                                >  1:2076590/1‑70 (MQ=255)
   tgatgaTGTCGGTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCa                                                              <  1:636126/69‑1 (MQ=255)
   tgatgaTGTCGGTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGt                                                            >  1:93500/1‑71 (MQ=255)
         tgTCGGTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCa                                                              >  1:2189147/1‑63 (MQ=255)
             ggTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGt                                                  >  1:77709/1‑71 (MQ=255)
                 aCCGTCTGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTAcc                                              >  1:2160855/1‑71 (MQ=255)
                 aCCGTCGGCACAGCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTAcc                                              >  1:790541/1‑71 (MQ=255)
                 aCCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTAcc                                              >  1:2070284/1‑71 (MQ=255)
                 aCCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTAcc                                              >  1:497712/1‑71 (MQ=255)
                        gCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTACCCAt                                           >  1:682563/1‑67 (MQ=255)
                        gCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTACCCATCAg                                        <  1:1924205/70‑1 (MQ=255)
                        gCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTACCCATCAGa                                       <  1:6214/71‑1 (MQ=255)
                        gCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCAt                                           <  1:719296/67‑1 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAg                                                             <  1:1082963/47‑1 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCa                                                              >  1:1520501/1‑46 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCa                                                              >  1:1666224/1‑46 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCa                                                              >  1:1881391/1‑46 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATc                                     >  1:421416/1‑71 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATc                                     >  1:995069/1‑71 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATc                                     >  1:819134/1‑71 (MQ=255)
                          acacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATc                                     >  1:1510788/1‑71 (MQ=255)
                           cacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTAcc                                              <  1:361095/61‑1 (MQ=255)
                           cacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTAcc                                              <  1:734022/61‑1 (MQ=255)
                           cacCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCt                                    <  1:1064387/71‑1 (MQ=255)
                            acCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCgct                                                       <  1:1046910/51‑1 (MQ=255)
                            acCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCgct                                                       <  1:1804547/51‑1 (MQ=255)
                            acCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCgct                                                       <  1:1090134/51‑1 (MQ=255)
                            acCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCgct                                                       <  1:1366975/51‑1 (MQ=255)
                             ccGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCGCTGCTGTTACCCATCAGATCttt                                  >  1:1349892/1‑71 (MQ=255)
                             ccGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCttt                                  <  1:2142827/71‑1 (MQ=255)
                             ccGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCttt                                  <  1:1402195/71‑1 (MQ=255)
                               gAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCgctgct                                                    >  1:247493/1‑51 (MQ=255)
                                aaTTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTAc                                                                 <  1:2202514/37‑1 (MQ=255)
                                aaTTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTAc                                                                 <  1:133103/37‑1 (MQ=255)
                                aaTTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTAc                                                                 <  1:183541/37‑1 (MQ=255)
                                    tctcATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTc                                 >  1:1616467/1‑65 (MQ=255)
                                        atcaGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCttt                                  >  1:725343/1‑60 (MQ=255)
                                         tcaGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCctgctg                      >  1:1879055/1‑71 (MQ=255)
                                          caGTAAAGTCAGACGCTTTCAGGTTGCGCAGTTCgc                                                        <  1:249669/36‑1 (MQ=255)
                                           aGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGAGCCTGCtgtg                    >  1:17788/1‑71 (MQ=255)
                                             tAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATc                                     <  1:679593/52‑1 (MQ=255)
                                              aaaGTCAGACGCTTTCAGGTTACGCAGTTCgctga                                                     >  1:1977767/1‑34 (MQ=255)
                                              aaaGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGt                              >  1:670909/1‑58 (MQ=255)
                                               aaGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCtgtg                    >  1:2036756/1‑67 (MQ=255)
                                                  tCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTtgctgc             <  1:1389911/71‑1 (MQ=255)
                                                    aGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCa                                >  1:583835/1‑50 (MQ=255)
                                                     gACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCtgt                     <  1:265106/60‑1 (MQ=255)
                                                     gACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCtgt                     <  1:480038/60‑1 (MQ=255)
                                                      aCGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAg          <  1:2077270/70‑1 (MQ=255)
                                                       cGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGaa        >  1:2176164/1‑71 (MQ=255)
                                                        gCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAg                               >  1:1776155/1‑47 (MQ=255)
                                                        gCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAg                               >  1:77567/1‑47 (MQ=255)
                                                        gCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTtgc                >  1:909213/1‑62 (MQ=255)
                                                        gCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGaa        >  1:198640/1‑70 (MQ=255)
                                                          tttCAGGTTACGCAGTTCGCTGCTGTTACCCATCAg                                        >  1:179431/1‑36 (MQ=255)
                                                          tttCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGa                                       >  1:1594411/1‑37 (MQ=255)
                                                          tttCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGaa        >  1:1501396/1‑68 (MQ=255)
                                                            tCAGGTTGCGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAg                               >  1:243219/1‑43 (MQ=255)
                                                             cAGGTTGCGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGAATGCGtt  >  1:1957752/1‑71 (MQ=255)
                                                             cAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGt                              >  1:1182380/1‑43 (MQ=255)
                                                             cAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGAATGCg    >  1:1624570/1‑69 (MQ=255)
                                                             cAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGAATGCg    >  1:2086580/1‑69 (MQ=255)
                                                             cAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTAAGt                              >  1:1576239/1‑43 (MQ=255)
                                                               ggTTACGCAGTTCGCTGCTGTTACCCATCAGATCttt                                  <  1:1080742/37‑1 (MQ=255)
                                                               ggTTACGCAGTTCGCTGCTGTTACCCATCAGATCttt                                  <  1:182069/37‑1 (MQ=255)
                                                                 ttACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAg          <  1:771837/59‑1 (MQ=255)
                                                                 ttACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAg          <  1:155397/59‑1 (MQ=255)
                                                                   |                                                                
CTTTGATGATGTCGGTTACCGTCGGCACACCGAATTTCTCATCAGTAAAGTCAGACGCTTTCAGGTTACGCAGTTCGCTGCTGTTACCCATCAGATCTTTCAGTGCCTGCTGTGTTGCTGCCAGAATGCGTT  >  minE/2681851‑2681982

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: