Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 788,075 G→A 28.6% P461S (CCA→TCA)  ycgB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE788,0750GA28.6% 61.0 / 26.9 42P461S (CCA→TCA) ycgBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/18);  new base A (6/6);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCA  >  minE/788009‑788143
                                                                  |                                                                    
cccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAat                                                                   <  1:548767/70‑1 (MQ=255)
 ccACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAatat                                                                 <  1:187626/71‑1 (MQ=255)
       gCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCa                                                           <  1:673941/71‑1 (MQ=255)
           tGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGAGt                                                       <  1:1441794/71‑1 (MQ=255)
              acaTGCTTCAGGACTTCTTTGCGCCCCCGATCTAGCGGTGCGCGATTATGTgg                                                                      <  1:3431034/53‑1 (MQ=255)
                aTGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCg                                                  <  1:1584187/71‑1 (MQ=255)
                 tGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGa                                                 <  1:3335757/71‑1 (MQ=255)
                   cTTCAGTACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGaa                                                <  1:533717/70‑1 (MQ=255)
                           cttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGt                                            <  1:1123299/66‑1 (MQ=255)
                           cttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTcgccg                                       >  1:156816/1‑71 (MQ=255)
                            ttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTcgccgc                                      >  1:1106077/1‑71 (MQ=255)
                            ttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTcgccgc                                      >  1:2192647/1‑71 (MQ=255)
                            ttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTcgccgc                                      >  1:1950734/1‑71 (MQ=255)
                            ttcttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTcgccg                                       >  1:2175364/1‑70 (MQ=255)
                              cttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCcgcg                                     >  1:230876/1‑70 (MQ=255)
                              cttTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCcgcg                                     >  1:2605247/1‑70 (MQ=255)
                                 tGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCaa                                  <  1:854101/70‑1 (MQ=255)
                                  gcgcCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTc                                                      >  1:2591487/1‑49 (MQ=255)
                                  gcgcCCCCGATCCAGCGGTGAGCGATTATGTGGAATATAACGCAGCGTc                                                      >  1:1508299/1‑49 (MQ=255)
                                    gcCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATcc                              <  1:1182818/71‑1 (MQ=255)
                                       cccTATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGt                                            <  1:2856702/54‑1 (MQ=255)
                                       cccGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCAc                            >  1:3502551/1‑70 (MQ=255)
                                         cGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGtt                         >  1:3847548/1‑71 (MQ=255)
                                          gATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCa                                   <  1:1748347/60‑1 (MQ=255)
                                           aTCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGt                                            <  1:3092849/50‑1 (MQ=255)
                                           aTCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTcc                       <  1:2894773/71‑1 (MQ=255)
                                            tCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTcgccgc                                      <  1:2031779/55‑1 (MQ=255)
                                            tCCAGCGGTGCGCGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTcc                       <  1:2979154/70‑1 (MQ=255)
                                              cAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAAc                                               >  1:3235093/1‑44 (MQ=255)
                                              cAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAg                     <  1:2354668/70‑1 (MQ=255)
                                               aGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTc                        >  1:340937/1‑66 (MQ=255)
                                                 cGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCt                 <  1:3056880/71‑1 (MQ=255)
                                                  ggggcgcGATTATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCt                 >  1:1141065/4‑70 (MQ=255)
                                                  ggTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTg                <  1:2609896/71‑1 (MQ=255)
                                                   gTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTg                >  1:186860/1‑70 (MQ=255)
                                                      cgcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTc                                           >  1:830455/1‑40 (MQ=255)
                                                      cgcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAata            <  1:926751/71‑1 (MQ=255)
                                                      cgcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAata            <  1:3437630/71‑1 (MQ=255)
                                                      cgcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGACCTGAata            <  1:3291789/71‑1 (MQ=255)
                                                       gcgATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAatat           >  1:3238249/1‑71 (MQ=255)
                                                           ttATGTGAAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCg        >  1:3184303/1‑70 (MQ=255)
                                                               gtgGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTc    <  1:3719746/70‑1 (MQ=255)
                                                                tgGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCa  <  1:3811341/71‑1 (MQ=255)
                                                                 ggAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTc    >  1:1054525/1‑68 (MQ=255)
                                                                  |                                                                    
CCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCA  >  minE/788009‑788143

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: