Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,073,671 C→T 41.9% G50S (GGC→AGC)  ydhQ ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,073,6710CT41.9% 14.9 / 34.8 31G50S (GGC→AGC) ydhQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (5/13);  new base T (2/11);  total (7/24)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

CGATAGCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGTA  >  minE/1073603‑1073735
                                                                    |                                                                
cGATAGCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAAccagccag                                                                >  1:927150/1‑71 (MQ=255)
    aGCGAGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGcccc                                                            >  1:503442/1‑71 (MQ=255)
     gCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGcccc                                                            <  1:888011/70‑1 (MQ=255)
       gTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGt                                                         <  1:2193253/71‑1 (MQ=255)
        tGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGt                                                         <  1:2036432/70‑1 (MQ=255)
         gCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGa                                                       >  1:740913/1‑71 (MQ=255)
          ccATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGAc                                                      <  1:219601/71‑1 (MQ=255)
               gTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAGCCAGCCAGCCCCCGTGACCACgg                                                 >  1:75508/1‑71 (MQ=255)
                ttACCGTACTGTGTGCATTCGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCACgg                                                 <  1:77254/70‑1 (MQ=255)
                   ccGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCACGGTAtt                                             <  1:924777/71‑1 (MQ=255)
                    cGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTg                                                        <  1:1919451/59‑1 (MQ=255)
                        cTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCACCCAGCTAGCCACCGTGACCACGGTATTTTCgg                                        >  1:1504816/1‑71 (MQ=255)
                        cTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:756330/66‑1 (MQ=255)
                         tgtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGt                                                         <  1:2042167/53‑1 (MQ=255)
                           tgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGct                                      <  1:489336/70‑1 (MQ=255)
                           tgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGc                                       <  1:866077/69‑1 (MQ=255)
                                aTTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtttt                                           <  1:40310/60‑1 (MQ=255)
                                aTTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtttt                                           <  1:1485545/60‑1 (MQ=255)
                                 tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGAc                                                      <  1:2145835/48‑1 (MQ=255)
                                 tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTc                                          <  1:1837104/60‑1 (MQ=255)
                                       gATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCtt                         >  1:2059271/1‑71 (MQ=255)
                                           ccACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTcaac                     <  1:352193/71‑1 (MQ=255)
                                              ccGGTTGCCACTTTCAACCAGCCAGCCCCCGTGAACACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAAt                  <  1:1503538/71‑1 (MQ=255)
                                               cGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATc                 <  1:1141735/71‑1 (MQ=255)
                                                     ccACTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCaa               <  1:1798116/67‑1 (MQ=255)
                                                     ccACTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTc            <  1:1743103/70‑1 (MQ=255)
                                                     ccACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTc            <  1:925347/70‑1 (MQ=255)
                                                          ttCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCtt       >  1:1335621/1‑70 (MQ=255)
                                                          ttCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCtt       >  1:2134004/1‑70 (MQ=255)
                                                             aaCCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGt   <  1:1547721/71‑1 (MQ=255)
                                                              aCCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGTa  <  1:409790/71‑1 (MQ=255)
                                                              aCCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGTa  <  1:152866/71‑1 (MQ=255)
                                                                    |                                                                
CGATAGCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGTA  >  minE/1073603‑1073735

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: