Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 236,687 A→G 13.2% S14S (TCA→TCG yaiE → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE236,6870AG13.2% 104.7 / 5.6 60S14S (TCA→TCGyaiEconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (2/6);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

ACCAAACAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGAA  >  minE/236628‑236747
                                                           |                                                            
aCCAAACAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCttt                                                     >  1:1989919/1‑69 (MQ=255)
aCCAAACAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCttt                                                     >  1:917440/1‑69 (MQ=255)
aCCAAACAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCttt                                                     >  1:835743/1‑69 (MQ=255)
aCCAAACAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCttt                                                     <  1:81387/69‑1 (MQ=255)
      cAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCcagc                                               >  1:2448507/1‑69 (MQ=255)
      cAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCcagc                                               >  1:467015/1‑69 (MQ=255)
      cAGGCCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCcagc                                               >  1:2134595/1‑69 (MQ=255)
             aGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCttt                                                     <  1:506561/56‑1 (MQ=255)
                 tATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCGATCGGCTTTTCCAGCAGCAGCACTgg                                    >  1:650227/1‑69 (MQ=255)
                 tATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCGATCGGCTTTTCCAGCAGCAGCACTgg                                    >  1:1843872/1‑69 (MQ=255)
                               aTGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCcagcagc                                            >  1:151732/1‑47 (MQ=255)
                               aTGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTg                       >  1:3229955/1‑68 (MQ=255)
                               aTGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTg                       >  1:2903828/1‑68 (MQ=255)
                                 gAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:3280520/1‑69 (MQ=255)
                                 gAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:1754091/1‑69 (MQ=255)
                                 gAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:1093168/1‑69 (MQ=255)
                                 gAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:1953182/1‑69 (MQ=255)
                                 gAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:940021/1‑69 (MQ=255)
                                  aGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:2748365/1‑68 (MQ=255)
                                  aGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:881657/1‑68 (MQ=255)
                                  aGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:2048043/1‑68 (MQ=255)
                                  aGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgt                    >  1:1921090/1‑68 (MQ=255)
                                    tACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgtgt                  <  1:809951/68‑1 (MQ=255)
                                    tACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGtt                 >  1:453761/1‑69 (MQ=255)
                                    tACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGtt                 >  1:868212/1‑69 (MQ=255)
                                    tACTTTTCCGGCAAAGTGAAATCAATCGGCTGTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGgtgt                  <  1:1290882/68‑1 (MQ=255)
                                      cTTTTCCTGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:779633/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:2746213/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:5736/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:1887137/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:3139255/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:2419339/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:3391113/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:454858/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:397191/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:2206004/43‑1 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTAt               >  1:1072428/1‑69 (MQ=255)
                                      cTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTAt               >  1:1927282/1‑69 (MQ=255)
                                      cTTTTCCGCAAAGTGAAATCAATCGGCTTTTCCagcagcagc                                         <  1:9137/42‑1 (MQ=37)
                                            ctgCAAAGTGAAATCGATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTg          <  1:1379075/66‑1 (MQ=255)
                                            cGGCAAAGTGAAATGGATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTg          <  1:368456/68‑1 (MQ=255)
                                            cGGCAAAGTGAAATCGATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTg          <  1:1112251/68‑1 (MQ=255)
                                            cGGCAAAGTGAAATCGATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTg          <  1:529136/68‑1 (MQ=255)
                                            cGGCAAAGTGAAATCGATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTg          <  1:1052921/68‑1 (MQ=255)
                                            cGGCAAAGTGAAATCGATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTg          <  1:1814808/68‑1 (MQ=255)
                                                 aaGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCg    <  1:1082514/69‑1 (MQ=255)
                                                 aaGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCg    <  1:1361910/69‑1 (MQ=255)
                                                 aaGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCg    <  1:216982/69‑1 (MQ=255)
                                                 aaGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCg    <  1:2862630/69‑1 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  <  1:870319/69‑1 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  <  1:1039683/69‑1 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  <  1:544252/69‑1 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  >  1:3411656/1‑69 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  >  1:1659029/1‑69 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  >  1:2860622/1‑69 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  >  1:250408/1‑69 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  >  1:2131775/1‑69 (MQ=255)
                                                   gTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGa   >  1:1039437/1‑68 (MQ=255)
                                                    tGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  <  1:3352117/68‑1 (MQ=255)
                                                    tGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGaa  <  1:2634900/68‑1 (MQ=255)
                                                           |                                                            
ACCAAACAGGGCCAGTTTATGCTTCAAAGTAATGAGTACTTTTCCGGCAAAGTGAAATCAATCGGCTTTTCCAGCAGCAGCACTGGTCGCGCCAGCGTGGGTGTTATGGTTGAAGGCGAA  >  minE/236628‑236747

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: