Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,017,691 A→G E50G (GAA→GGA)  ydgA → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,017,6910AG95.1% 138.6 / ‑1.3 41E50G (GAA→GGA) ydgAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (13/26);  total (14/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTGTTG  >  minE/1017632‑1017758
                                                           |                                                                   
aaCCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTg                                                           >  1:2066694/1‑70 (MQ=255)
 aCCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTgg                                                          >  1:2059297/1‑70 (MQ=255)
    cATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGt                                                      <  1:1268780/71‑1 (MQ=255)
    cATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGt                                                      <  1:2377/71‑1 (MQ=255)
     aTCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGt                                                      >  1:2074397/1‑70 (MQ=255)
         cGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGt                                                  >  1:204373/1‑70 (MQ=255)
          gAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATcc                                                                 <  1:697060/54‑1 (MQ=255)
          gAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAg                                                   <  1:595869/68‑1 (MQ=255)
            aGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCaa                                                               >  1:202822/1‑54 (MQ=255)
                  ggTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTAt                                               <  1:2202351/64‑1 (MQ=255)
                   gTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACt                                        <  1:898538/70‑1 (MQ=255)
                   gTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACt                                        <  1:2232911/70‑1 (MQ=255)
                   gTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACt                                        <  1:1730951/70‑1 (MQ=255)
                    ttgcgcAGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCaaa                                           <  1:1137612/64‑1 (MQ=255)
                          aGGCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGaa                                                        >  1:1092540/1‑47 (MQ=255)
                            gCGAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGtt                                                 >  1:1774474/1‑52 (MQ=255)
                              gAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGtt                                                 <  1:539414/50‑1 (MQ=255)
                              gAACGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCaaaa                                          >  1:2378874/1‑57 (MQ=255)
                               aaCGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGt                           <  1:2132201/71‑1 (MQ=255)
                                aCGCGCAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTAt                                               >  1:1693290/1‑50 (MQ=255)
                                  gcgcAACTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTAtt                        <  1:1887345/71‑1 (MQ=255)
                                       aCTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTcag                     >  1:1925052/1‑69 (MQ=255)
                                       aCTCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTcag                     >  1:577762/1‑69 (MQ=255)
                                         tCAAACTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATc                                              <  1:1403556/42‑1 (MQ=255)
                                           aaaCTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACt                                        <  1:1722069/46‑1 (MQ=255)
                                           aaaCTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCagcagc                 >  1:1422031/1‑69 (MQ=255)
                                           aaaCTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCagcagc                 >  1:575847/1‑69 (MQ=255)
                                           aaaCTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGcc                <  1:1365124/70‑1 (MQ=255)
                                            aaCTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAg              <  1:1234547/71‑1 (MQ=255)
                                            aaCTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAg              <  1:952516/71‑1 (MQ=255)
                                             aCTGACAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGt             <  1:487590/71‑1 (MQ=255)
                                                 aCAGCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCa         <  1:2012538/71‑1 (MQ=255)
                                                    gCTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACt      <  1:579025/71‑1 (MQ=255)
                                                     cTCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTg     >  1:1029628/1‑71 (MQ=255)
                                                      tCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATcgc                                <  1:1842642/43‑1 (MQ=255)
                                                      tCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTcagca                   <  1:1903616/56‑1 (MQ=255)
                                                      tCCTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTcag                     <  1:888738/54‑1 (MQ=255)
                                                        cTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTca                      <  1:2133772/51‑1 (MQ=255)
                                                        cTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTGtt   <  1:1216765/70‑1 (MQ=255)
                                                        cTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTGt    <  1:2102787/69‑1 (MQ=255)
                                                        cTGGATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTGTTg  <  1:2185763/71‑1 (MQ=255)
                                                           |                                                                   
AACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACAGCTCCTGAATCCAACCTGGAAGTGAGTTATCAAAACTATCATCGCGGCGTATTCAGCAGCCAGTTGCAACTGTTG  >  minE/1017632‑1017758

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: