Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,534,049 G→A 50.7% P64P (CCG→CCA yfcN → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,534,0490GA50.7% ‑6.5 / 67.4 73P64P (CCG→CCAyfcNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (9/27);  new base A (21/16);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 8.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAA  >  minE/1533987‑1534114
                                                              |                                                                 
gTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttatta                                                             <  1:2261356/69‑1 (MQ=255)
 tGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttatta                                                             <  1:1273006/68‑1 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                           >  1:2666129/1‑69 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                           >  1:2022318/1‑69 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                           <  1:1283595/69‑1 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                           >  1:3252062/1‑69 (MQ=255)
              ccAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTcc                                               <  1:3146584/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:1856532/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:2957544/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:1771809/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:1791028/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:1069770/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:908151/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:970670/69‑1 (MQ=255)
                aGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                              <  1:2234424/68‑1 (MQ=255)
                    acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         <  1:1514270/68‑1 (MQ=255)
                    acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         <  1:2584306/68‑1 (MQ=255)
                    acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         <  1:3106075/68‑1 (MQ=255)
                     cAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGaaa                                        <  1:2199187/69‑1 (MQ=255)
                     cAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGaaa                                        <  1:1333011/69‑1 (MQ=255)
                      aGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         <  1:337820/67‑1 (MQ=255)
                        gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                     <  1:1551542/69‑1 (MQ=255)
                        gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                     <  1:1593718/69‑1 (MQ=255)
                        gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                     <  1:2190577/69‑1 (MQ=255)
                        gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                     <  1:2186767/69‑1 (MQ=255)
                         cTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                     <  1:1869473/68‑1 (MQ=255)
                                     aatattTTTCGGACGAATTTCTGCCATTATTGAATACCGAAGGTCCg                                              <  1:242573/45‑1 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:973003/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:2965009/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:86578/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:583138/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:383993/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:2438598/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:2484031/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:3095200/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:2661651/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:2199969/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:3034801/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:2218377/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:1072854/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:1149314/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:128879/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:1348984/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:1365489/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:1412495/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:1485954/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:2224488/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                         >  1:1729705/1‑49 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                     >  1:1188298/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                     >  1:2737347/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                     >  1:1433596/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                     >  1:3201637/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                     >  1:162614/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGGTATTAAATACCGAAGGTCCGGTg                                           >  1:2533740/1‑47 (MQ=255)
                                             tcCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCa                 <  1:3212765/68‑1 (MQ=255)
                                                     gTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCg                            <  1:1695137/49‑1 (MQ=255)
                                                       ttCAGCCATTATTGAATACCGAAGGTCCGGTCAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCg       <  1:2707855/68‑1 (MQ=255)
                                                       ttCAGCCATTATTGAATACCGAAGGTCCGGTCAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCg       <  1:3063179/68‑1 (MQ=255)
                                                       ttCAGCCATTATTGAATACCGAAGGTCCGGTCAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCg       <  1:1362953/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:1619204/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:863651/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:1346116/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:3074387/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:1429718/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:2696780/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:1615601/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:2565434/68‑1 (MQ=255)
                                                            ccGTTTTTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa  <  1:1154672/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa  <  1:2504684/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa  <  1:1510733/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa  <  1:2607978/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa  <  1:1779427/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga   <  1:1113906/67‑1 (MQ=255)
                                                              |                                                                 
GTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAA  >  minE/1533987‑1534114

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: