Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 10,173 C→T 62.8% V108I (GTA→ATA)  yaaH ← conserved inner membrane protein associated with acetate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE10,1730CT62.8% 30.0 / 36.0 59V108I (GTA→ATA) yaaHconserved inner membrane protein associated with acetate transport
Reads supporting (aligned to +/- strand):  ref base C (16/6);  new base T (25/12);  total (41/18)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTC  >  minE/10111‑10237
                                                              |                                                                
gAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGccc                                                              <  1:1107016/67‑1 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGc                                                                >  1:3159935/1‑36 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:1022340/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:996881/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:716791/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:697834/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:400655/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:371084/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:2896434/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:287016/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:2512604/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:2431423/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:2381442/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:2354664/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:230538/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:1854823/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:1265770/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:1131765/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCcaca                                                          >  1:1696296/1‑42 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCca                                                            >  1:2519463/1‑40 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACtgtg                               >  1:1826372/1‑69 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACtgtg                               >  1:2932861/1‑69 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACtgtg                               >  1:2970159/1‑69 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACtgtg                               >  1:524670/1‑69 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACtgtg                               >  1:883331/1‑69 (MQ=255)
                             tCAGCGTGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACtgtg                               >  1:2093726/1‑69 (MQ=255)
                              cAGCGTGCCTAAGAACATAAACAGCGTAAATACGCCCCACAGCCCCAGGTa                                                >  1:248283/1‑51 (MQ=39)
                              cAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGc                              >  1:1054963/1‑69 (MQ=255)
                              cAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGc                              >  1:1050950/1‑69 (MQ=255)
                              cAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGc                              >  1:738552/1‑69 (MQ=255)
                              cAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGc                              >  1:1187677/1‑69 (MQ=255)
                              cAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGc                              >  1:2760571/1‑69 (MQ=255)
                              cAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGc                              >  1:2770324/1‑69 (MQ=255)
                                gCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCa                             <  1:1914279/68‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATCTGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:2194833/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:588004/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:2843564/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:2906336/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:1363072/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:1854639/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:2204124/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:2689221/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:2300181/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:176482/69‑1 (MQ=255)
                                   tGCCGAAGAACATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:1400429/69‑1 (MQ=255)
                                     ccGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTgg                        >  1:1308750/1‑68 (MQ=255)
                                     ccGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTgg                        >  1:3141064/1‑68 (MQ=255)
                                     ccGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTgg                        >  1:2642267/1‑68 (MQ=255)
                                     ccGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTgg                        >  1:24701/1‑68 (MQ=255)
                                     ccGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACAAAGGAACTGTGCa                             >  1:680969/1‑63 (MQ=255)
                                        aagaaCATAAACAGCGTAAATATGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTg                         <  1:2881305/64‑1 (MQ=255)
                                        aagaaCATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCa                    <  1:138647/69‑1 (MQ=255)
                                        aagaaCATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCa                    <  1:3347993/69‑1 (MQ=255)
                                        aagaaCATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCa                    <  1:2117607/69‑1 (MQ=255)
                                        aagaaCATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCAGTTGGCGCa                    <  1:100621/69‑1 (MQ=255)
                                                          aaTACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTc  >  1:1604328/1‑69 (MQ=255)
                                                          aaTACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTc  >  1:627085/1‑69 (MQ=255)
                                                          aaTACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTc  >  1:408221/1‑69 (MQ=255)
                                                          aaTACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTc  >  1:2762827/1‑69 (MQ=255)
                                                              |                                                                
GAAAACGAATTGCAGAACGCGTGCGCCTTTCAGCGTGCCGAAGAACATAAACAGCGTAAATACGCCCCACAGACCCAGGTAGACACCAAGGAACTGTGCATTTGGCGCATCGGTCAGACCCAGTTTC  >  minE/10111‑10237

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: