Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 110,575 G→A 55.3% P47P (CCG→CCA aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE110,5750GA55.3% 5.5 / 36.0 47P47P (CCG→CCAaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (8/13);  new base A (16/10);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGC  >  minE/110509‑110640
                                                                  |                                                                 
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:1013154/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:918112/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:914082/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:819563/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:73921/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:683246/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:554147/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:3031856/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:2725584/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:2663022/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:2019370/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:1938202/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:173519/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:1415393/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:1401380/1‑69 (MQ=255)
cGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTACCAtc                                                                 >  1:1185292/1‑69 (MQ=255)
                       cGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTa                                             >  1:1490026/1‑66 (MQ=255)
                       cGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTa                                             >  1:1301574/1‑66 (MQ=255)
                       cGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTa                                             >  1:3023321/1‑66 (MQ=255)
                               aCCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa                                  >  1:963238/1‑69 (MQ=255)
                               aCCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa                                  >  1:1125461/1‑69 (MQ=255)
                               aCCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCaaa                                  >  1:3340893/1‑69 (MQ=255)
                                    aGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGtctc                             <  1:3414991/69‑1 (MQ=255)
                                     gAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTctct                            <  1:1663411/69‑1 (MQ=255)
                                     gAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTctct                            <  1:2557225/69‑1 (MQ=255)
                                     gAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTctct                            <  1:854161/69‑1 (MQ=255)
                                     gAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTctct                            <  1:3285055/69‑1 (MQ=255)
                                     gAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTctct                            <  1:1216884/69‑1 (MQ=255)
                                     gAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTctct                            <  1:413188/69‑1 (MQ=255)
                                                    tctATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAcc             <  1:577087/69‑1 (MQ=255)
                                                    tctATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAcc             <  1:2263724/69‑1 (MQ=255)
                                                    tctATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAcc             <  1:119481/69‑1 (MQ=255)
                                                    tctATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAcc             <  1:1312221/69‑1 (MQ=255)
                                                    tctATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAcc             <  1:3078005/69‑1 (MQ=255)
                                                    tctATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAAcc             <  1:3014741/69‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:525385/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:1258319/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:474318/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:673307/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:6950/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:2528587/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:786945/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:233784/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:1329964/67‑1 (MQ=255)
                                                          gAAGTACCATCTCCGAAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGa         <  1:100659/67‑1 (MQ=255)
                                                               tCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgcgc  >  1:886005/1‑69 (MQ=255)
                                                               tCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgcg   >  1:1523150/1‑68 (MQ=255)
                                                                  |                                                                 
CGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGC  >  minE/110509‑110640

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: