Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,489,398 A→G 100% Y545Y (TAT→TAC nuoL ← NADH:ubiquinone oxidoreductase, membrane subunit L

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,489,3980AG77.9% 45.0 / 8.5 36Y545Y (TAT→TACnuoLNADH:ubiquinone oxidoreductase, membrane subunit L
Reads supporting (aligned to +/- strand):  ref base A (0/8);  new base G (19/9);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 8.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGCCAGTCAAATCCCCAGGCGTTGTACCACCAGGTGCCCAGCAGACGGCCCGGCGCACTGTTGGC  >  minE/1489332‑1489466
                                                                  |                                                                    
tGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCAt                                                                     <  1:2934272/68‑1 (MQ=255)
tGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCAt                                                                     <  1:1849073/68‑1 (MQ=255)
tGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCAt                                                                     <  1:2652996/68‑1 (MQ=255)
tGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCAt                                                                     <  1:2540482/68‑1 (MQ=255)
tGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCAt                                                                     <  1:770548/68‑1 (MQ=255)
     tCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGcc                                                               <  1:195361/69‑1 (MQ=255)
      cAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGcc                                                               <  1:1348585/68‑1 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCa                                                              >  1:3154870/1‑65 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:593957/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:348130/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:710791/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:336304/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:766871/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:2717257/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:2704502/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:2589998/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:2473164/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:2244493/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:815172/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:1991072/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:191426/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:1909529/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:1793406/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:1647925/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:1577902/1‑67 (MQ=255)
          ggATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGt                                                            >  1:1501830/1‑67 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTAGAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:2736385/69‑1 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:826679/69‑1 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:1254922/69‑1 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:945168/69‑1 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:3402798/69‑1 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:23020/69‑1 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:2051408/69‑1 (MQ=255)
                            cAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:1276519/69‑1 (MQ=255)
                                                       aaCACTTTGTCGTACAGCCAGTCAAAGCCCCAGGCGTTGTac                                        <  1:3324798/42‑1 (MQ=37)
                                                         cacTTTGTCATACAGCCAGTCAAATCCCCAGGCGTTGTACCACCAGGTGCCCAGCAGACGGCCCGgcgc           <  1:174886/69‑1 (MQ=255)
                                                                  aTACAGCCAGTCAAATCCCCAGGCGTTGTACCACCAGGTGCCCAGCAGACGGCCCGGCGCACTGTTGGc  <  1:538903/69‑1 (MQ=255)
                                                                  |                                                                    
TGAGTTCAGCGGATCGCGTTTCAGCAACCAGGCAATACCCAGGAACGGCTTGACGAACACTTTGTCATACAGCCAGTCAAATCCCCAGGCGTTGTACCACCAGGTGCCCAGCAGACGGCCCGGCGCACTGTTGGC  >  minE/1489332‑1489466

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: