Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 435,666 T→C 55.0% A146A (GCA→GCG ybfF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE435,6660TC55.0% 10.9 / 46.1 60A146A (GCA→GCGybfFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/5);  new base C (23/10);  total (45/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AACAAAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTC  >  minE/435603‑435728
                                                               |                                                              
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGGTGCTGGCGCATTATtgctgctgc                                                           >  1:1821053/1‑69 (MQ=255)
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgctgc                                                           >  1:263422/1‑69 (MQ=255)
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgctgc                                                           >  1:990399/1‑69 (MQ=255)
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgctgc                                                           >  1:411683/1‑69 (MQ=255)
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgctgc                                                           >  1:3352355/1‑69 (MQ=255)
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgctgc                                                           >  1:3198036/1‑69 (MQ=255)
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgctgc                                                           >  1:2290082/1‑69 (MQ=255)
aacaaAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgc                                                              >  1:1410178/1‑66 (MQ=255)
     aaGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATtgctgct                                                             >  1:1831667/1‑62 (MQ=255)
     aaGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCtgcttgct                                                      >  1:2180823/1‑69 (MQ=255)
     aaGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCtgcttgct                                                      >  1:300334/1‑69 (MQ=255)
     aaGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCtgcttgct                                                      >  1:2781371/1‑69 (MQ=255)
     aaGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCtgcttgct                                                      >  1:637932/1‑69 (MQ=255)
     aaGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCtgcttgc                                                       >  1:2543251/1‑68 (MQ=255)
     aaGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCtgcttgc                                                       >  1:1835780/1‑68 (MQ=255)
              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)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCtgcttgct                                                      >  1:2676515/1‑39 (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:1029594/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:343495/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:283873/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:2744700/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:1347215/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:2226006/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:337376/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:539121/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:181143/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:687698/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGActg                        >  1:1671225/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:1490669/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGAct                         >  1:256404/1‑67 (MQ=255)
                                   cttcATTAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGAct                         >  1:1410193/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:1343416/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:653625/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:1716926/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:560281/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:466449/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:462420/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:1347190/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:2289266/66‑1 (MQ=255)
                                        ttAAGATGCTGGCGCATTATTGCCGCTGCTTGCTGGCGTGTTTGTGCGTCTGATTCACTGACtgcg                      <  1:2630284/66‑1 (MQ=255)
                                          aaGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGtt                    >  1:644266/1‑66 (MQ=255)
                                            gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                <  1:838916/68‑1 (MQ=255)
                                            gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                <  1:360311/68‑1 (MQ=255)
                                            gATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATa                <  1:677383/68‑1 (MQ=255)
                                                        ttattGATGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtct   >  1:225877/1‑69 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:1193360/1‑68 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:92441/1‑68 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:16481/1‑68 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:1821557/1‑68 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:2001632/1‑68 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:209260/1‑68 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:220209/1‑68 (MQ=255)
                                                          attGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAAtctc  >  1:3200257/1‑68 (MQ=255)
                                                               |                                                              
AACAAAAGATTTCAGCAGAAACTGAATCACCCCTTCTTCATTAAGATGCTGGCGCATTATTGCTGCTGCTTGCTGGCGAGTTTGTGCGTCCGATTCACTGACCGCGTTGATAGCCGCAAAAATCTC  >  minE/435603‑435728

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: