Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 435,693 C→T 54.7% S137S (TCG→TCA ybfF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE435,6930CT54.7% 12.6 / 71.5 86S137S (TCG→TCAybfFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/18);  new base T (20/27);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

CTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACGGCGTACGTGATAGTCGACCGGCGC  >  minE/435638‑435761
                                                       |                                                                    
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGAtt                                                                  >  1:814097/1‑60 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGAt                                                                   >  1:1840360/1‑59 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:283873/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:343495/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:1497523/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:337376/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:1487041/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:1671225/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:1347215/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:2744700/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:181143/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:687698/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:539121/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:1029594/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                                         >  1:2226006/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGAct                                                          >  1:1490669/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGAct                                                          >  1:3321612/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGAct                                                          >  1:1410193/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGAct                                                          >  1:256404/1‑67 (MQ=255)
cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTAACTGAct                                                          >  1:1680903/1‑67 (MQ=255)
  tcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:939566/68‑1 (MQ=255)
   cATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACcgcg                                                       >  1:806109/1‑68 (MQ=255)
   cATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACcgcg                                                       >  1:2681675/1‑68 (MQ=255)
   cATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACcgcg                                                       >  1:2858807/1‑68 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGtt                                                     <  1:1617657/68‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:653625/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:2630284/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:2289266/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:1716926/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:1347190/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:1343416/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:462420/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:466449/66‑1 (MQ=255)
     ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:560281/66‑1 (MQ=255)
       aaGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGtt                                                     >  1:644266/1‑66 (MQ=255)
         gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                                                 <  1:360311/68‑1 (MQ=255)
         gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                                                 <  1:677383/68‑1 (MQ=255)
         gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                                                 <  1:838916/68‑1 (MQ=255)
                     ttattGATGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtct                                    >  1:225877/1‑69 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:3200257/1‑68 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:220209/1‑68 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:16481/1‑68 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:92441/1‑68 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:209260/1‑68 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:1821557/1‑68 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:1193360/1‑68 (MQ=255)
                       attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                                   >  1:2001632/1‑68 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:1895712/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:3327397/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:1526313/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:415967/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:3097338/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:3083196/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:1331095/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:2974239/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:1327194/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:2066711/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:1187917/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:2616984/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:2155526/67‑1 (MQ=255)
                            cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:2174042/67‑1 (MQ=255)
                             gctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:1778075/67‑1 (MQ=255)
                             gctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                              <  1:1584395/67‑1 (MQ=255)
                                  ttgctGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                                       <  1:3014410/37‑1 (MQ=25)
                                   tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:3248816/1‑66 (MQ=255)
                                   tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:3386499/1‑66 (MQ=255)
                                   tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:3172434/1‑66 (MQ=255)
                                   tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:545270/1‑66 (MQ=255)
                                   tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:207521/1‑66 (MQ=255)
                                   tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:1107987/1‑66 (MQ=255)
                                   tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:989951/1‑66 (MQ=255)
                                   tgctGACGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat                       >  1:1084782/1‑66 (MQ=255)
                                          gAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACGGCGTACGTGa                <  1:3412107/68‑1 (MQ=255)
                                          gAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACGGCGTACGTGa                <  1:1432032/68‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:2296815/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:958672/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:931137/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:551952/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:3121224/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:2345073/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:2296471/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:2121702/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:2057124/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:2031801/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:2023608/42‑1 (MQ=255)
                                                   cGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCtcat                                 <  1:1039455/42‑1 (MQ=255)
                                                       cGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACGGCGTACGTGATAGTCGACCGgcgc  >  1:2831147/1‑69 (MQ=255)
                                                       |                                                                    
CTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACGGCGTACGTGATAGTCGACCGGCGC  >  minE/435638‑435761

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: