Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 277,701 T→C 54.8% P470P (CCA→CCG cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE277,7010TC54.8% 5.2 / 30.9 42P470P (CCA→CCGcyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (5/14);  new base C (15/8);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGG  >  minE/277638‑277759
                                                               |                                                          
gAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATCCCCATAAAGCCCAACGCATACAGCGGCa                                                        >  1:2457894/1‑68 (MQ=255)
gAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATCCCCATAAAGCCCAACGCATACAGCGGCa                                                        >  1:294297/1‑68 (MQ=255)
gAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGACATCCCCATAAAGCCCAACGCATACAGCGGCa                                                        >  1:2228543/1‑68 (MQ=255)
  aCTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATaa                                                     >  1:839857/1‑69 (MQ=255)
  aCTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATaa                                                     >  1:2755710/1‑69 (MQ=255)
  aCTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATaa                                                     >  1:1584434/1‑69 (MQ=255)
                tgctgGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACgaagaa                                        >  1:1504726/1‑68 (MQ=255)
                    ggCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGcc                                     <  1:3232909/67‑1 (MQ=255)
                    ggCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGcc                                     <  1:1347073/67‑1 (MQ=255)
                    ggCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGcc                                     <  1:1785909/67‑1 (MQ=255)
                    ggCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGcc                                     <  1:601612/67‑1 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:1920311/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:288183/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:2479776/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:1238561/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:2414349/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:502257/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:3075000/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:3124969/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:1883551/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:3184512/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaagaa                                        >  1:1581516/1‑55 (MQ=255)
                             gacgGGTCATCCCCATAAAGCCCAACGCATACAGCGGCATAAAGGCAACgaaga                                         >  1:3270767/1‑54 (MQ=255)
                                  gTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATccagaaccag                     <  1:912561/69‑1 (MQ=255)
                                  gTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATccagaaccag                     <  1:2495186/69‑1 (MQ=255)
                                          cATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAAcgcgcg             <  1:1970500/69‑1 (MQ=255)
                                          cATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAAcgcgcg             <  1:1411905/69‑1 (MQ=255)
                                          cATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAAcgcgcg             <  1:1316688/69‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCTGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:797998/66‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCGGCATAACGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:2600628/66‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:1886031/66‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:172430/66‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:1529480/66‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:3290157/66‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:1142164/66‑1 (MQ=255)
                                           ataaaGCCCAACGCATACAGCGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGt            <  1:70002/66‑1 (MQ=255)
                                           aTGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAAcgcgcg             <  1:3425754/68‑1 (MQ=255)
                                            tGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGtt           <  1:124443/69‑1 (MQ=255)
                                              aaGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTa         >  1:1017869/1‑69 (MQ=255)
                                                     gcgcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAgg  <  1:3346134/69‑1 (MQ=255)
                                                     gcgcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAgg  <  1:884966/69‑1 (MQ=255)
                                                     gcgcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAgg  <  1:2145948/69‑1 (MQ=255)
                                                               |                                                          
GAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGG  >  minE/277638‑277759

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: