Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 271,637 C→T 31.6% P158P (CCG→CCA yajL ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE271,6370CT31.6% 33.4 / 17.1 41P158P (CCG→CCAyajLconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/17);  new base T (9/4);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCAT  >  minE/271569‑271697
                                                                    |                                                            
aCTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTAccc                                                              >  1:2105923/1‑69 (MQ=255)
  ttCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg                                                            >  1:2279723/1‑69 (MQ=255)
  ttCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg                                                            >  1:1490731/1‑69 (MQ=255)
  ttCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg                                                            >  1:3099124/1‑69 (MQ=255)
  ttCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg                                                            >  1:274009/1‑69 (MQ=255)
  ttCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCgg                                                            >  1:187890/1‑69 (MQ=255)
    cATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGc                                                           >  1:2283538/1‑68 (MQ=255)
        ggCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCctg                                                      >  1:86976/1‑69 (MQ=255)
                   gCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCa                                           <  1:1602088/69‑1 (MQ=255)
                   gCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCa                                           <  1:3067980/69‑1 (MQ=255)
                   gCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCa                                           <  1:2965995/69‑1 (MQ=255)
                       aaccaacAGGTCTATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAAttt                                       <  1:995366/69‑1 (MQ=255)
                       aaccaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAAttt                                       <  1:1938537/69‑1 (MQ=255)
                       aaccaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAAttt                                       <  1:1406772/69‑1 (MQ=255)
                       aaccaacAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAAttt                                       <  1:2611776/69‑1 (MQ=255)
                               ggTCGATAATTTTCAGACCATAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                               <  1:307461/69‑1 (MQ=255)
                               ggTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                               <  1:839483/69‑1 (MQ=255)
                               ggTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                               <  1:2113016/69‑1 (MQ=255)
                               ggTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                               <  1:2729762/69‑1 (MQ=255)
                               ggTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                               <  1:2395171/69‑1 (MQ=255)
                               ggTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                               <  1:2675016/69‑1 (MQ=255)
                                         tttCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTg                               <  1:2180007/59‑1 (MQ=255)
                                              taCCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAAttt                                       <  1:1266995/45‑1 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAAtt                                        >  1:2720001/1‑42 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:1471011/1‑69 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:1975870/1‑69 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:2892804/1‑69 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:1790325/1‑69 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:1752506/1‑69 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:3137229/1‑69 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:3145926/1‑69 (MQ=255)
                                                 cAAAGTCGATAGCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCt             >  1:3290581/1‑69 (MQ=255)
                                                          tAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGc     <  1:1547914/68‑1 (MQ=255)
                                                          tAGCTGTACCCGGCCCCTGGCTGATCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGc     <  1:258735/68‑1 (MQ=255)
                                                            gCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCTGCCAt  <  1:1789716/69‑1 (MQ=255)
                                                            gCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCTGCCAt  <  1:788678/69‑1 (MQ=255)
                                                            gCTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCTGCCAt  <  1:711320/69‑1 (MQ=255)
                                                            gCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCAt  >  1:1929537/1‑69 (MQ=255)
                                                            gCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCAt  >  1:28467/1‑69 (MQ=255)
                                                            gCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCAt  >  1:2230662/1‑69 (MQ=255)
                                                            gCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCAt  >  1:1196131/1‑69 (MQ=255)
                                                             cTGTACCTGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCTGCCAt  <  1:2627596/68‑1 (MQ=255)
                                                                    |                                                            
ACTTCATGGGCTTTTTCACGCCCAACCAACAGGTCGATAATTTTCAGACCAAAGTCGATAGCTGTACCCGGCCCCTGGCTGGTCAGCAATTTTACCCGTGCATCCCAGACGACGCGCTTGTCCAGCCAT  >  minE/271569‑271697

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: