Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,073,671 C→T 24.4% 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,6710CT24.4% 83.3 / 25.8 45G50S (GGC→AGC) ydhQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/20);  new base T (6/5);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AGCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGTA  >  minE/1073607‑1073735
                                                                |                                                                
aGCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGccc                                                             >  1:2625274/1‑70 (MQ=255)
 gCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGcccc                                                            <  1:3333390/70‑1 (MQ=255)
 gCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGcccc                                                            <  1:1677669/70‑1 (MQ=255)
 gCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGcccc                                                            <  1:670067/70‑1 (MQ=255)
          tGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACg                                                  >  1:3595290/1‑71 (MQ=255)
           gTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACgg                                                 >  1:3331846/1‑71 (MQ=255)
           gTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACgg                                                 >  1:2833418/1‑71 (MQ=255)
               ccGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCATGGTAtt                                             <  1:2573750/71‑1 (MQ=255)
                       tgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGctc                                     >  1:2271957/1‑71 (MQ=255)
                         tgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGctct                                    <  1:2173261/70‑1 (MQ=255)
                             tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGAc                                                      <  1:2576868/48‑1 (MQ=255)
                             tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:1337583/57‑1 (MQ=255)
                             tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:3829716/57‑1 (MQ=255)
                             tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:3257521/57‑1 (MQ=255)
                             tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:2906338/57‑1 (MQ=255)
                             tttGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGc                               <  1:2782318/71‑1 (MQ=255)
                               tGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGa                                                       >  1:2351355/1‑45 (MQ=255)
                               tGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGa                                                       >  1:1126862/1‑45 (MQ=255)
                               tGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGa                                                       >  1:2608161/1‑45 (MQ=255)
                               tGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGa                                                       >  1:659245/1‑45 (MQ=255)
                               tGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtttt                                           >  1:2501768/1‑57 (MQ=255)
                               tGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtttt                                           >  1:1432946/1‑57 (MQ=255)
                                   gATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:2750098/51‑1 (MQ=255)
                                   gATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:342550/51‑1 (MQ=255)
                                   gATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:2378052/51‑1 (MQ=255)
                                   gATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:2289192/51‑1 (MQ=255)
                                     tCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTc                        >  1:1832922/1‑70 (MQ=255)
                                           cGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtt                                             <  1:2284530/43‑1 (MQ=255)
                                           cGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATc                 <  1:3514359/71‑1 (MQ=255)
                                                    cTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:909830/1‑71 (MQ=255)
                                                    cTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:1363319/1‑71 (MQ=255)
                                                    cTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:565291/1‑71 (MQ=255)
                                                    cTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:3641074/1‑71 (MQ=255)
                                                    cTTTCAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:2257858/1‑71 (MQ=255)
                                                    cTTTCAACCAGCCAGCCCCCGTGAGCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:3440647/1‑71 (MQ=255)
                                                    cTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:4102344/1‑71 (MQ=255)
                                                    cTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCt        >  1:4085024/1‑71 (MQ=255)
                                                       tCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGctct                                    <  1:712814/40‑1 (MQ=255)
                                                       tCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGctct                                    <  1:2413712/40‑1 (MQ=255)
                                                       tCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGctct                                    <  1:2052061/40‑1 (MQ=255)
                                                       tCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCtt       <  1:3332576/69‑1 (MQ=255)
                                                        cAACCAGCTAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCtt       <  1:1033640/68‑1 (MQ=255)
                                                        cAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCtt       <  1:3949419/68‑1 (MQ=255)
                                                        cAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCtt       <  1:3776220/68‑1 (MQ=255)
                                                          aCCAGCTAGCCCCCGTGACCACGCTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGTa  <  1:2557705/71‑1 (MQ=255)
                                                            cagctagcCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATc                 <  1:2493597/54‑1 (MQ=255)
                                                            cagctagcCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATc                 <  1:2799182/54‑1 (MQ=255)
                                                            cagctagcCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATc                 <  1:634307/54‑1 (MQ=255)
                                                            cagctagcCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATc                 <  1:2455946/54‑1 (MQ=255)
                                                            cagctagcCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATc                 <  1:2705264/54‑1 (MQ=255)
                                                            cagctagcCCCCGTGACCACGGTATTTACGGCTCTTCCGCCAGCTTCAACAATc                 <  1:3452078/54‑1 (MQ=255)
                                                                |                                                                
AGCGTGCCATTGTTACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTCGGCTCTTCCGCCAGCTTCAACAATCAACTCGCCTTCAGTA  >  minE/1073607‑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: