Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,869,590 T→C 39.3% G425G (GGT→GGC rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,869,5900TC39.3% 56.6 / 84.5 122G425G (GGT→GGCrnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base T (30/44);  new base C (15/33);  total (45/77)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-04

Only 100 of 123 total aligned reads displayed.

TAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGT  >  minE/2869527‑2869653
                                                               |                                                               
tAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTggg                                                              <  1:1582199/67‑1 (MQ=255)
tAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTggg                                                              <  1:783846/67‑1 (MQ=255)
tAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTggg                                                              <  1:2302205/67‑1 (MQ=255)
tAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTggg                                                              <  1:614968/67‑1 (MQ=255)
tAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTggg                                                              <  1:2530817/67‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:1222832/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:2703456/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:137895/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:2879571/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:2951362/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:3287603/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:3373778/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:2050331/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:1237159/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:695021/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:793540/68‑1 (MQ=255)
    cATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:793702/68‑1 (MQ=255)
        tcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCAt                                                       >  1:2291832/1‑66 (MQ=255)
        tcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCAt                                                       >  1:2299101/1‑66 (MQ=255)
        tcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCAt                                                       >  1:2814573/1‑66 (MQ=255)
        tcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCAt                                                       >  1:216352/1‑66 (MQ=255)
             gagTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGAtctc                                                         <  1:713458/59‑1 (MQ=255)
                 tGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCgaag                                           <  1:150396/69‑1 (MQ=255)
                 tGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCgaag                                           <  1:894516/69‑1 (MQ=255)
                       aCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                     >  1:2424533/1‑69 (MQ=255)
                       aCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                     >  1:1926868/1‑69 (MQ=255)
                       aCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                     >  1:3038893/1‑69 (MQ=255)
                       aCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                     >  1:592728/1‑69 (MQ=255)
                       aCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                     >  1:823712/1‑69 (MQ=255)
                       aCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                     >  1:1452588/1‑69 (MQ=255)
                       aCCTCTATAAAGTCCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGa                                     >  1:1716197/1‑69 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGaagaag                                        >  1:3320427/1‑59 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                               >  1:1308331/1‑68 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                               >  1:1588773/1‑68 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAt                              >  1:353624/1‑69 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAt                              >  1:2884680/1‑69 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAt                              >  1:3205369/1‑69 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAt                              >  1:1016558/1‑69 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAt                              >  1:1021067/1‑69 (MQ=255)
                              taAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAt                              >  1:1988325/1‑69 (MQ=255)
                                aaGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                               <  1:3307441/66‑1 (MQ=255)
                                aaGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                               >  1:716173/1‑66 (MQ=255)
                                aaGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                               <  1:3387492/66‑1 (MQ=255)
                                aaGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtt                             >  1:1737831/1‑68 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:2548154/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:765191/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:1038790/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:1439690/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:307993/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:1065382/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:268170/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:2638645/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:3165069/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:699137/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:1859578/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:2153271/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:1963595/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:1090545/68‑1 (MQ=255)
                                  gTGCTGGATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGAGAAGAAGCGAAGTTCAtttt                           <  1:112503/68‑1 (MQ=255)
                                    gCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCaa                        <  1:1269001/69‑1 (MQ=255)
                                     cTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                               <  1:2612216/61‑1 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACg                      >  1:736279/1‑67 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACg                      >  1:2886604/1‑67 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:1743196/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:3007625/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:3023564/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:2709341/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:2632937/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:2135834/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:1824635/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:2386613/1‑68 (MQ=255)
                                        gATAAAGCCCGTGAAGAACGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGc                     >  1:499615/1‑68 (MQ=255)
                                              gCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCaa                        <  1:3179664/59‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCgaag                                           <  1:3254070/38‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg            <  1:1069960/69‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg            <  1:735143/69‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg            <  1:1139742/69‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg            <  1:694659/69‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg            <  1:604104/69‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg            <  1:3159869/69‑1 (MQ=255)
                                                ccGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCg            <  1:1307599/69‑1 (MQ=255)
                                                     aagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTg       <  1:1551982/69‑1 (MQ=255)
                                                     aagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTg       >  1:2082315/1‑69 (MQ=255)
                                                     aagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTg       >  1:3324846/1‑69 (MQ=255)
                                                     aagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTg       <  1:1302235/69‑1 (MQ=255)
                                                     aagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTg       >  1:1687435/1‑69 (MQ=255)
                                                       gaaCGCGGTGGGATCTTATTTGAGAGCGAAGAAGCGAAGTTCa                               <  1:219570/43‑1 (MQ=255)
                                                       gaaCGCGGCGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttt                           <  1:123917/47‑1 (MQ=255)
                                                          cgctGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:1246257/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:1969058/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:1745516/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:1808144/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:1663973/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:1189467/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:374970/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:2355507/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg   <  1:2497491/68‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTAAACGCCGTATTGAACGt  <  1:1463388/69‑1 (MQ=255)
                                                          cgcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTAAACGCCGTATTGAACGt  <  1:2885151/69‑1 (MQ=255)
                                                           gcggcgGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCAtttg                           <  1:2300780/43‑2 (MQ=255)
                                                               |                                                               
TAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGT  >  minE/2869527‑2869653

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: