Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 278,904 C→T 64.2% L69L (TTG→TTA cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE278,9040CT64.2% 41.3 / 39.1 67L69L (TTG→TTAcyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base T (28/15);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 7.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

GAGCCTGCTGGCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACGG  >  minE/278858‑278959
                                              |                                                       
gagCCTGCTGGCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACgatga                                   <  1:3295949/69‑1 (MQ=255)
gagCCTGCTGGCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACgatga                                   <  1:2075149/69‑1 (MQ=255)
gagCCTGCTGGCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACAACACAATCGCCACgatga                                   <  1:922717/69‑1 (MQ=255)
   ccTGCTGGCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGata                                 >  1:928628/1‑68 (MQ=255)
   ccTGCTGGCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGata                                 >  1:1770453/1‑68 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:2642316/67‑1 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:999159/67‑1 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:3115307/67‑1 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:2968220/67‑1 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:2928116/67‑1 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:2583941/67‑1 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:1114356/67‑1 (MQ=255)
         ggCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:1200972/67‑1 (MQ=255)
            tACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATATATGATa                        >  1:2245009/1‑68 (MQ=255)
            tACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATa                        >  1:2704193/1‑68 (MQ=255)
            tACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATa                        >  1:2437405/1‑68 (MQ=255)
            tACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATa                        >  1:3153137/1‑68 (MQ=255)
            tACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATa                        >  1:2148368/1‑68 (MQ=255)
            tACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATa                        >  1:1783802/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAACCCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:2624016/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACgatgat                                  >  1:2650699/1‑50 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAgg                  >  1:3029980/1‑66 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:897355/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1740130/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:2766094/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1406489/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:3030597/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:3108862/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:3147027/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:736943/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:796815/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:3285032/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:781851/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1334822/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:451692/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:561616/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:606915/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1546592/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1663470/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1747734/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1632872/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1892522/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1975872/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:246991/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:2626674/1‑68 (MQ=255)
                    atAATGGCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:1030168/1‑68 (MQ=255)
                    atAATGCCGTCAGCAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCg                >  1:5731/1‑68 (MQ=255)
                       aTGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:778482/53‑1 (MQ=255)
                       aTGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACAt                            <  1:64424/53‑1 (MQ=255)
                       aTGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTa            >  1:1258025/1‑69 (MQ=255)
                       aTGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTa            >  1:2652101/1‑69 (MQ=255)
                            gTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTc       <  1:2211906/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:2789557/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:1460579/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:2922444/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:1419033/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:1972855/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:750535/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:3020340/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:200126/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:2036302/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:542311/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:2120399/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:2371379/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:3141656/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:1374916/69‑1 (MQ=255)
                                 cAAAACCACGCAGTAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACgg  <  1:2547975/69‑1 (MQ=255)
                                              |                                                       
GAGCCTGCTGGCTACGCATCATAATGGCGTCAGCAAAACCACGCAGCAACATCACAATCGCCACGATGATATACATGATACCGAGGCGTTTATGGTCGACGG  >  minE/278858‑278959

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