Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 720,741 A→G 42.5% V961A (GTC→GCC)  rne ← fused ribonucleaseE endoribonuclease and scaffold for formation of degradosome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE720,7410AG42.5% 18.2 / 50.9 69V961A (GTC→GCC) rnefused ribonucleaseE endoribonuclease and scaffold for formation of degradosome
Reads supporting (aligned to +/- strand):  ref base A (7/32);  new base G (17/12);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 7.80e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

CGACTACTGGCGCGGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACT  >  minE/720676‑720806
                                                                 |                                                                 
cGACTACTGGCGCGGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACa                                                                 <  1:202668/68‑1 (MQ=255)
cGACTACTGGCGCGGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACa                                                                 <  1:816411/68‑1 (MQ=255)
cGACTACTGGCGCGGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACa                                                                 <  1:1874826/68‑1 (MQ=255)
cGACTACTGGCGCGGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACa                                                                 <  1:2382218/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCTACAACTTCTTCAat                                                      <  1:1539985/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGCCAACTTCTTCAat                                                      <  1:162527/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:100355/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:508703/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:2207445/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:2977031/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:1839677/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:2272585/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:2296945/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:2315060/68‑1 (MQ=255)
           gcgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:79463/68‑1 (MQ=255)
            cgGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:1369979/67‑1 (MQ=255)
                            aacTTCAGGTTCAGCAACTTCAACCTCAGCAGTTTCGACATCTTCTTCAat                                                      <  1:1537038/51‑1 (MQ=255)
                                 cAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAActtc                                                            >  1:3160881/1‑40 (MQ=255)
                                 cAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCgg                                 >  1:908073/1‑67 (MQ=255)
                                 cAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGt                                >  1:737882/1‑68 (MQ=255)
                                 cAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGt                                >  1:1193390/1‑68 (MQ=255)
                                   ggTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAat                                                      <  1:2052594/44‑1 (MQ=255)
                                   ggTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTc                              <  1:621707/68‑1 (MQ=255)
                                   ggTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTc                              <  1:1351700/68‑1 (MQ=255)
                                   ggTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTc                              <  1:3083223/68‑1 (MQ=255)
                                   ggTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTc                              <  1:3346841/68‑1 (MQ=255)
                                    gTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTc                              <  1:2745595/67‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCTCCg                      <  1:2162500/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:3239556/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:207003/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:1959192/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:566674/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:2252695/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:1791205/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:1767009/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:716435/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:280935/69‑1 (MQ=255)
                                          cAACCACAACCTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:2800564/69‑1 (MQ=255)
                                           aaCTACAACCTCAGCTGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCgg                     <  1:1964885/69‑1 (MQ=255)
                                           aaCTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCgg                     <  1:739118/69‑1 (MQ=255)
                                           aaCTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCgg                     <  1:2394705/69‑1 (MQ=255)
                                           aaCTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCgg                     <  1:2756832/69‑1 (MQ=255)
                                           aaCTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCAACgg                     <  1:942005/69‑1 (MQ=255)
                                               aCAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1204328/1‑69 (MQ=255)
                                               aCAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:2012560/1‑69 (MQ=255)
                                                    cTCAGCAGTTTCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCg                      <  1:1982305/59‑1 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCACCCACCGGGTCt                 >  1:2215482/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:400494/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1118124/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:896367/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1205229/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1229419/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1430470/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1476721/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1513170/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1905131/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:1968394/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:39724/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:3119507/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:2220511/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:2814148/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:2337941/1‑54 (MQ=255)
                                                              tCGGCAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCt                 >  1:2791673/1‑54 (MQ=255)
                                                              tCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCa     <  1:2639344/66‑1 (MQ=255)
                                                              tCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCa     <  1:563586/66‑1 (MQ=255)
                                                              tCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCa     <  1:1589694/66‑1 (MQ=255)
                                                              tCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCa     <  1:2014107/66‑1 (MQ=255)
                                                              tCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCa     <  1:2631808/66‑1 (MQ=255)
                                                              tCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACt  >  1:392383/1‑69 (MQ=255)
                                                                 |                                                                 
CGACTACTGGCGCGGCAGGTTGAGCAACAACTTCAGGTTCAGCAACTACAACCTCAGCAGTTTCGACAACTTCTTCAATATCTGCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACT  >  minE/720676‑720806

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: