Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,493,288 G→A 63.1% Y47Y (TAC→TAT nuoH ← NADH:ubiquinone oxidoreductase, membrane subunit H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,493,2880GA63.1% 22.9 / 27.1 46Y47Y (TAC→TATnuoHNADH:ubiquinone oxidoreductase, membrane subunit H
Reads supporting (aligned to +/- strand):  ref base G (6/11);  new base A (9/20);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

AACATTTTGATCATGTCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCGCAGGTGACA  >  minE/1493226‑1493354
                                                              |                                                                  
aaCATTTTGATCATGTCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACgg                                                              <  1:1664660/69‑1 (MQ=255)
 aCATTTTGATCATGTCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACgg                                                              <  1:3341502/68‑1 (MQ=255)
 aCATTTTGATCATGTCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACgg                                                              <  1:2734104/68‑1 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:2754968/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:962431/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:954771/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:2767180/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:2716391/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:2266703/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:1917553/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:191301/1‑69 (MQ=255)
  cATTTTGATCATGTCCTCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCATAACGGtt                                                            >  1:155987/1‑69 (MQ=255)
               tCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGAccc                                               <  1:825354/69‑1 (MQ=255)
               tCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGAccc                                               >  1:2273047/1‑69 (MQ=255)
               tCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGAccc                                               <  1:366604/69‑1 (MQ=255)
               tCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGAccc                                               <  1:127405/69‑1 (MQ=255)
               tCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGAccc                                               <  1:2827237/69‑1 (MQ=255)
               tCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGAccc                                               >  1:3073139/1‑69 (MQ=255)
               tCCGCAACCAGCTGGAGAGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCt                                                          >  1:2122215/1‑58 (MQ=255)
                ccGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGAccc                                               <  1:1523239/68‑1 (MQ=255)
                                          cAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATg                    <  1:1184871/69‑1 (MQ=255)
                                          cAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATg                    <  1:1512946/69‑1 (MQ=255)
                                          cAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATg                    <  1:593048/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:1260344/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:1850053/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:3404406/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:371625/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:405278/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:439410/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:622235/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:623959/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:1188458/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:3332184/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:3250976/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:3209761/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:3207521/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:1345456/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:1356532/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:2415172/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:2180312/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:2159548/69‑1 (MQ=255)
                                                       tAGGTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:1557010/69‑1 (MQ=255)
                                                         ggTCCATAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCACAgg       <  1:1034254/67‑1 (MQ=255)
                                                             cGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCGCAGGTGACa  >  1:3088021/1‑68 (MQ=255)
                                                             cGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCGCAGGTGACa  >  1:3021693/1‑68 (MQ=255)
                                                             cGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCGCAGGTGACa  >  1:585965/1‑68 (MQ=255)
                                                              |                                                                  
AACATTTTGATCATGTCCGCAACCAGCTGGAGCGAACCGCCCCAGCCAACACGGTTAGGTCCGTAACGGTTCTGGAACAGACCCAGCAGGCGACGTTCGCCAAAGCTCATGAATGCCCCGCAGGTGACA  >  minE/1493226‑1493354

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: