Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 276,399 T→C 53.4% G39G (GGA→GGG cyoD ← cytochrome o ubiquinol oxidase subunit IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE276,3990TC53.4% 8.8 / 62.3 75G39G (GGA→GGGcyoDcytochrome o ubiquinol oxidase subunit IV
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base C (30/10);  total (48/27)
Fisher's exact test for biased strand distribution p-value = 5.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCA  >  minE/276335‑276460
                                                                |                                                             
gCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCt                                                            >  1:2734635/1‑68 (MQ=255)
gCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCt                                                            >  1:2655000/1‑68 (MQ=255)
gCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCt                                                            >  1:3208602/1‑68 (MQ=255)
   ccAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:1808277/1‑68 (MQ=255)
   ccAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:1636406/1‑68 (MQ=255)
   ccAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:86200/1‑68 (MQ=255)
                   tcccATTGCCAGGATTGTTCCCAGAATTACGTCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACgg                                        <  1:1120567/67‑1 (MQ=255)
                   tGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCATAACgg                                        <  1:2350410/69‑1 (MQ=255)
                   tGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACgg                                        <  1:70532/69‑1 (MQ=255)
                   tGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACgg                                        <  1:541959/69‑1 (MQ=255)
                   tGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACgg                                        <  1:3127611/69‑1 (MQ=255)
                   tGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACgg                                        <  1:2382931/69‑1 (MQ=255)
                   tGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACgg                                        <  1:909637/69‑1 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:1338941/1‑51 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:3328155/1‑51 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:1583307/1‑51 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:867325/1‑51 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGt                                                          >  1:1022151/1‑50 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:2382829/69‑1 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:903774/69‑1 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:1840326/69‑1 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:126686/69‑1 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:1621119/69‑1 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:626487/69‑1 (MQ=255)
                    gCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:2708635/69‑1 (MQ=255)
                     ccATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGa                                       <  1:27233/68‑1 (MQ=255)
                          gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGtc                                                         >  1:1661302/1‑45 (MQ=255)
                          gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACgg                                        >  1:2575906/1‑62 (MQ=255)
                          gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCAcc                                 >  1:3419621/1‑69 (MQ=255)
                          gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCAcc                                 >  1:1089207/1‑69 (MQ=255)
                          gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCAcc                                 >  1:2636070/1‑69 (MQ=255)
                          gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCAcc                                 >  1:195024/1‑69 (MQ=255)
                          gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCAcc                                 >  1:1823938/1‑69 (MQ=255)
                                        cAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCg                    <  1:816387/68‑1 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCa                                   >  1:218287/1‑42 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGt                               >  1:1368793/1‑46 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:96836/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:2855937/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:3147481/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:854943/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:3389861/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:453422/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:465164/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:525227/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:1682272/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:2013220/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:2747997/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:2124505/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:1766632/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:2201190/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:226480/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:1435717/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:2340572/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:1139794/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAgg               >  1:1087054/1‑62 (MQ=255)
                                                   ccGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATAGACAg                >  1:3065345/1‑61 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGcc       >  1:434664/1‑69 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGcc       >  1:1278651/1‑69 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGcc       >  1:1227188/1‑69 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGcc       >  1:55077/1‑69 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGcc       >  1:1981522/1‑69 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGcc       >  1:419309/1‑69 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGcc       >  1:411417/1‑69 (MQ=255)
                                                    cGGAGAGGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGc        >  1:836827/1‑68 (MQ=255)
                                                       agagGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGt    <  1:2016061/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:2101951/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:3404555/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:2318978/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:152300/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:2322743/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:1320383/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:3087224/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:2954455/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:2360527/69‑1 (MQ=255)
                                                         agGCAGCCCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCa  <  1:2780193/69‑1 (MQ=255)
                                                                |                                                             
GCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCA  >  minE/276335‑276460

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: