Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 435,681 A→T 61.3% T141T (ACT→ACA ybfF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE435,6810AT61.3% 50.7 / 73.0 80T141T (ACT→ACAybfFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/7);  new base T (22/27);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 5.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACGGCGTACGTGA  >  minE/435617‑435747
                                                                |                                                                  
gcagAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGtt                                                                 <  1:555535/68‑1 (MQ=255)
     aaCTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGttt                                                                >  1:2506266/1‑64 (MQ=255)
     aaCTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCg                                                           >  1:343916/1‑69 (MQ=255)
     aaCTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCg                                                           >  1:2697338/1‑69 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGtt                                                                 >  1:869965/1‑47 (MQ=39)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTctg                                                       >  1:1775197/1‑55 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGAtt                                                    >  1:814097/1‑60 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGAt                                                     >  1:1840360/1‑59 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:1671225/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:1029594/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:181143/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:2226006/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:2744700/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:283873/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:337376/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:687698/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:343495/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:539121/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:1347215/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:1497523/1‑67 (MQ=255)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                                           >  1:1487041/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)
                     cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGAct                                            >  1:1490669/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:2858807/1‑68 (MQ=255)
                        cATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACcgcg                                         >  1:2681675/1‑68 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGtt                                       <  1:1617657/68‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:653625/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:560281/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:466449/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:2289266/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:462420/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:1716926/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:2630284/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:1343416/66‑1 (MQ=255)
                          ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:1347190/66‑1 (MQ=255)
                            aaGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGtt                                       >  1:644266/1‑66 (MQ=255)
                              gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                                   <  1:838916/68‑1 (MQ=255)
                              gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                                   <  1:677383/68‑1 (MQ=255)
                              gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                                   <  1:360311/68‑1 (MQ=255)
                                          ttattGATGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtct                      >  1:225877/1‑69 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:1821557/1‑68 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:220209/1‑68 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:16481/1‑68 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:209260/1‑68 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:3200257/1‑68 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:92441/1‑68 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:2001632/1‑68 (MQ=255)
                                            attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc                     >  1:1193360/1‑68 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:2155526/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:1526313/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:2174042/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:2066711/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:415967/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:1187917/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:1327194/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:2616984/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:1331095/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:3327397/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:1895712/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:3097338/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:3083196/67‑1 (MQ=255)
                                                 cgctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:2974239/67‑1 (MQ=255)
                                                  gctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:1584395/67‑1 (MQ=255)
                                                  gctgctTGCTGGCGTGTTTGTGCGTCTGATTCACTGACTGCGTTGATGGCCGCGAAAATCTcatcat                <  1:1778075/67‑1 (MQ=255)
                                                       ttgctGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                                         <  1:3014410/37‑1 (MQ=25)
                                                        tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat         >  1:989951/1‑66 (MQ=255)
                                                        tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat         >  1:545270/1‑66 (MQ=255)
                                                        tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat         >  1:3386499/1‑66 (MQ=255)
                                                        tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat         >  1:3248816/1‑66 (MQ=255)
                                                        tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat         >  1:3172434/1‑66 (MQ=255)
                                                        tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat         >  1:207521/1‑66 (MQ=255)
                                                        tgctGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACggat         >  1:1107987/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)
                                                                |                                                                  
GCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTCATCATGACGGCGTACGTGA  >  minE/435617‑435747

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: