Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 277,857 C→T 50.0% A418A (GCG→GCA cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE277,8570CT50.0% ‑3.4 / 64.7 72A418A (GCG→GCAcyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (23/13);  new base T (18/18);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCA  >  minE/277799‑277924
                                                          |                                                                   
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:45656/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:1562227/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:1611324/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:2121952/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:2489311/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:1486266/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:3102999/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:475348/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:708903/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:831914/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:1087168/69‑1 (MQ=255)
aGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:962801/69‑1 (MQ=255)
 ggTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAgg                                                            >  1:1059342/1‑67 (MQ=255)
 ggTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAgg                                                            >  1:812902/1‑67 (MQ=255)
 ggTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAgg                                                            >  1:2547693/1‑67 (MQ=255)
 ggTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAgg                                                            >  1:813074/1‑67 (MQ=255)
 ggTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAgg                                                            >  1:1878886/1‑67 (MQ=255)
 ggTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAgg                                                            >  1:1371483/1‑67 (MQ=255)
 ggTCATCCCTGCGAAGCAGCCGAAGACCACGACGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:2538575/68‑1 (MQ=255)
  gTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:1551257/67‑1 (MQ=255)
  gTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:1204541/67‑1 (MQ=255)
    cATCCATGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAg                                                       >  1:1807732/1‑69 (MQ=255)
        ccTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCt                                                    >  1:3260174/1‑68 (MQ=255)
        ccTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCt                                                    >  1:2148004/1‑68 (MQ=255)
          tGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGt                                                  <  1:281594/68‑1 (MQ=255)
          tGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGt                                                  <  1:265764/68‑1 (MQ=255)
          tGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGt                                                  <  1:1112483/68‑1 (MQ=255)
          tGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGt                                                  <  1:574522/68‑1 (MQ=255)
                gcagcCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAg                                           <  1:1596178/69‑1 (MQ=255)
                gcagcCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAg                                           <  1:155202/69‑1 (MQ=255)
                gcagcCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAg                                           <  1:1861836/69‑1 (MQ=255)
                     cGAAGACCACGCCGCCGATGATCACGTTATTGAAGTGTGCAATCAGGa                                                           <  1:415804/48‑1 (MQ=255)
                      gAAGACCACGCCGCCGATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:2826236/47‑1 (MQ=255)
                        aGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGt                                   <  1:572929/69‑1 (MQ=255)
                        aGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGt                                   <  1:481044/69‑1 (MQ=255)
                        aGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGt                                   <  1:477085/69‑1 (MQ=255)
                        aGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGt                                   <  1:1722204/69‑1 (MQ=255)
                        aGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGt                                   <  1:359242/69‑1 (MQ=255)
                                  gccgATGATCACGTTATGGAAGTGTGCAATCAGGa                                                           <  1:1837810/35‑1 (MQ=255)
                                  gccgATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCgg                         <  1:3323640/69‑1 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:3106964/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:583993/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:848715/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:3401228/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:366559/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:391034/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:1095651/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:472939/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:3258190/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:2335646/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:2346713/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:1529577/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:2189749/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:153421/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:2626286/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:2682345/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAct                     >  1:3214500/1‑68 (MQ=255)
                                      atgatCACGTTATGGAAGTGTGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCTGCGCCAGGTAc                      >  1:1542650/1‑68 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCTGTACGGCcagca              >  1:2262033/1‑64 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:2090756/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:1129242/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:638641/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:1655294/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:943973/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:1021026/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:1309900/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:251656/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:1499833/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:291489/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:2970617/1‑65 (MQ=255)
                                                  tGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagcag             >  1:1380933/1‑65 (MQ=255)
                                                       gTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGt    >  1:456715/1‑69 (MQ=255)
                                                         gcgcAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCa  <  1:335552/69‑1 (MQ=255)
                                                         gcgcAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCa  <  1:327100/69‑1 (MQ=255)
                                                          |                                                                   
AGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCA  >  minE/277799‑277924

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