Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,760,795 A→G 52.1% A62A (GCA→GCG yrdC → dsRNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,760,7950AG52.1% 9.5 / 108.8 123A62A (GCA→GCGyrdCdsRNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (39/25);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

Only 100 of 124 total aligned reads displayed.

ACAGCAGTGATGCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGCAGCGTG  >  minE/2760730‑2760856
                                                                 |                                                             
aCAGCAGTGTTGCGACTGTTGGAGTTCAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCgcagca                                                            <  1:3354897/69‑1 (MQ=255)
        gATGCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACg                                                     >  1:3195780/1‑68 (MQ=255)
           gCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGagc                                                  <  1:1013786/68‑1 (MQ=255)
           gCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGagc                                                  <  1:1910215/68‑1 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGGTGATTTTAATCGCGGCAAATTa                                                       >  1:1320051/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:1934005/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:1030398/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:2275766/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:326893/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:2479856/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:169633/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:80157/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:3150712/1‑61 (MQ=255)
             gACTGTTGGAGTTAAAACGGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTa                                                       >  1:2813039/1‑61 (MQ=255)
                tGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:3314133/69‑1 (MQ=255)
                tGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:304916/69‑1 (MQ=255)
                    ggggTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTaaa                                          <  1:1571319/64‑1 (MQ=255)
                    ggAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGagc                                                  <  1:1453672/59‑1 (MQ=255)
                      aGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAAccc                                       <  1:2136160/68‑1 (MQ=255)
                      aGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAAccc                                       <  1:3212645/68‑1 (MQ=255)
                       gTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCta                                     <  1:3208709/69‑1 (MQ=255)
                          aaaaCAGCGTCCGGTTGAGAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCtata                                   >  1:390986/1‑68 (MQ=255)
                            aaCAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATAtt                                 <  1:2610808/68‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2718633/69‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2530608/69‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2354809/69‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:978926/69‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:330632/69‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:2159048/69‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:578063/69‑1 (MQ=255)
                              cAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:1052111/69‑1 (MQ=255)
                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:2783391/2‑68 (MQ=255)
                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:206033/2‑68 (MQ=255)
                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:1944204/2‑68 (MQ=255)
                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:126409/2‑68 (MQ=255)
                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:2990435/2‑68 (MQ=255)
                               ggCGTCCGGTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATtgat                              >  1:945391/2‑68 (MQ=255)
                                   tCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGacac                         <  1:1814937/69‑1 (MQ=255)
                                   tCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGacac                         <  1:455096/69‑1 (MQ=255)
                                   tCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGacac                         <  1:1631637/69‑1 (MQ=255)
                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:3133682/68‑1 (MQ=255)
                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:775690/68‑1 (MQ=255)
                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:133840/68‑1 (MQ=255)
                                       gTTGATAAGGGGCTGATTTTAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCAt                      <  1:1192586/68‑1 (MQ=255)
                                          gATATGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtg                  >  1:1489191/1‑69 (MQ=255)
                                          gATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAAccc                                       <  1:44542/48‑1 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgatga                            >  1:1921638/1‑55 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgtctg              >  1:18554/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:1175092/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:1603630/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:1599749/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:2884431/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:36580/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:3369985/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:2701093/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:909443/1‑69 (MQ=255)
                                              aGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTtgactg              >  1:2407560/1‑69 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:783531/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:1028413/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:398459/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:3408756/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:3397746/37‑1 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:12298/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:2955407/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:1889122/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:969439/1‑37 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            <  1:3079634/37‑1 (MQ=255)
                                                gggCTGATTTTAATCGCAGCAAATTACGAGCAGCTTa                                            >  1:2427949/1‑37 (MQ=255)
                                                    tGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATtgat                              <  1:192822/47‑1 (MQ=255)
                                                    tGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGc        >  1:2718308/1‑69 (MQ=255)
                                                    tGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGc        >  1:411601/1‑69 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2363331/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2458767/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:1368092/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:3215605/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:218035/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:3042222/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2492263/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2087600/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:3386591/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2880236/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2513726/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:2573513/1‑68 (MQ=255)
                                                         ttAATCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCg    >  1:363324/1‑68 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:459331/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1022107/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:203824/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:3414437/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:3079589/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2851150/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2850958/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2838599/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1649086/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2793584/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1766700/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2663355/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2558602/68‑1 (MQ=255)
                                                           aaTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2054368/68‑1 (MQ=255)
                                                           aaTCGCGGCAAAATACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:2397528/68‑1 (MQ=255)
                                                            aTCGCGGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGCTGACTGACGCGCAGCGTg  <  1:1892762/67‑1 (MQ=255)
                                                                 |                                                             
ACAGCAGTGATGCGACTGTTGGAGTTAAAACAGCGTCCGGTTGATAAGGGGCTGATTTTAATCGCAGCAAATTACGAGCAGCTTAAACCCTATATTGATGACACCATGTTGACTGACGTGCAGCGTG  >  minE/2760730‑2760856

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: