Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,534,055 A→G 55.3% L66L (TTA→TTG yfcN → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,534,0550AG55.3% 14.2 / 51.4 67L66L (TTA→TTGyfcNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/24);  new base G (21/16);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 2.83e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCG  >  minE/1533987‑1534120
                                                                    |                                                                 
gTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttatta                                                                   <  1:2261356/69‑1 (MQ=255)
 tGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGttatta                                                                   <  1:1273006/68‑1 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                                 >  1:2022318/1‑69 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                                 >  1:2666129/1‑69 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                                 >  1:3252062/1‑69 (MQ=255)
  gAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTaaa                                                                 <  1:1283595/69‑1 (MQ=255)
              ccAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTcc                                                     <  1:3146584/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:1069770/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:1791028/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:1856532/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:1771809/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:908151/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:970670/69‑1 (MQ=255)
               cAGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:2957544/69‑1 (MQ=255)
                aGGAACAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCg                                                    <  1:2234424/68‑1 (MQ=255)
                    acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               <  1:3106075/68‑1 (MQ=255)
                    acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               <  1:1514270/68‑1 (MQ=255)
                    acAGGCCGATGCCAGCCACTATTTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               <  1:2584306/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:2190577/69‑1 (MQ=255)
                        gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                           <  1:1551542/69‑1 (MQ=255)
                        gCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAatat                                           <  1:1593718/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:2224488/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:2661651/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:973003/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:2199969/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:86578/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:583138/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:2218377/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:383993/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:3095200/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:2438598/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:2484031/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:3034801/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:2965009/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:1072854/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:1149314/1‑49 (MQ=255)
                                        attTTTCGGACGAATTTCAGCCATTATTGAATACCGAAGGTCCGGTCaa                                               >  1:1729705/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)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                           >  1:1433596/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                           >  1:1188298/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                           >  1:3201637/1‑69 (MQ=255)
                                        attTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCa                           >  1:2737347/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:2565434/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:2696780/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:1346116/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:863651/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:1429718/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:1619204/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:3074387/68‑1 (MQ=255)
                                                           gCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:1615601/68‑1 (MQ=255)
                                                            ccGTTTTTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa        <  1:1154672/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa        <  1:2504684/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa        <  1:1779427/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa        <  1:2607978/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaagaa        <  1:1510733/68‑1 (MQ=255)
                                                            ccGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCgaaga         <  1:1113906/67‑1 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:774439/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:1788819/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:709485/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:520871/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:1200591/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:2282599/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:3137782/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:1859499/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:2940162/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:281739/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:1438039/1‑69 (MQ=255)
                                                                 attaAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTgcg  >  1:259949/1‑69 (MQ=255)
                                                                    |                                                                 
GTGAAACGCTTGATCCAGGAGCAGGCTGATGCCAGCCATTATTTCTCCGATGAGTTTCAGCCGTTATTAAATACCGAAGGTCCGGTGAAATATGTTCGCCCGGATGTCAGCCATTTTGAGGCGAAGAAACTGCG  >  minE/1533987‑1534120

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: