Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,486,723 G→A 61.9% L197L (CTG→TTG)  nuoN ← NADH:ubiquinone oxidoreductase, membrane subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,486,7230GA61.9% 37.1 / 49.1 76L197L (CTG→TTG) nuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base G (8/21);  new base A (22/25);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTGGAA  >  minE/1486658‑1486778
                                                                 |                                                       
aaaGGCGATAATCGCCAGCACCACGCGAATCGCTTCTCTGTCACCCACCGGTGCGTAGAGGAACAGACg                                                      >  1:2617519/1‑69 (MQ=255)
aaaGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACg                                                      <  1:712642/69‑1 (MQ=255)
aaaGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACg                                                      <  1:15284/69‑1 (MQ=255)
 aaGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACAGGTGCGTAGAGGAACAGACg                                                      <  1:69177/68‑1 (MQ=255)
    gCGATATTCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:438766/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGGACAAACGCATc                                                  >  1:981785/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCAt                                                   >  1:2325328/1‑68 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCAt                                                   >  1:662874/1‑68 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:1497475/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:231331/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:2400080/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:2522969/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:2607136/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:2667625/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:2820034/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:29555/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:1937003/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:1826388/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:297999/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:3208336/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:173632/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:3237005/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:459161/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:505729/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:738223/1‑69 (MQ=255)
    gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                                                  >  1:1491609/1‑69 (MQ=255)
      gATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcac                                                <  1:1377205/69‑1 (MQ=255)
         aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac                                             >  1:787501/1‑69 (MQ=255)
         aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac                                             >  1:814620/1‑69 (MQ=255)
         aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac                                             >  1:2626965/1‑69 (MQ=255)
         aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac                                             >  1:2170444/1‑69 (MQ=255)
               cAGCACACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGaa                                       <  1:1964516/68‑1 (MQ=255)
                   accacGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATa                                   >  1:3127515/1‑69 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:600953/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:3068188/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:3207294/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:17444/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:2302196/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:3273639/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:1592695/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:3239046/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:1561000/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:1496257/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:1289592/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:115691/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:962193/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:110586/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:1812172/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:2424066/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:2435626/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:247200/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:2216542/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:1081768/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:2774316/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:279011/69‑1 (MQ=255)
                      acGCGAATCGCTTCACTGTCACCCACCGGTGCGAAAAGGATCAAACGCATCACCACGCCGAAGATAGCg                                <  1:1732979/69‑1 (MQ=255)
                      acGCGAATCGCATCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:3426052/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGt      >  1:3431148/1‑68 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGt      >  1:2012329/1‑68 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:2788790/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:1295097/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:2024425/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:483735/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:809528/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:1676854/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:3402853/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:3401745/69‑1 (MQ=255)
                                                 ggTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:3044969/69‑1 (MQ=255)
                                                  gTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:766368/68‑1 (MQ=255)
                                                  gTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:2727173/68‑1 (MQ=255)
                                                  gTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:2523899/68‑1 (MQ=255)
                                                  gTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:1648819/68‑1 (MQ=255)
                                                  gTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:3285515/68‑1 (MQ=255)
                                                   tGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTg     <  1:508886/67‑1 (MQ=255)
                                                    gCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTGGaa  <  1:1970127/69‑1 (MQ=255)
                                                     cGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg                                <  1:112152/38‑1 (MQ=25)
                                                                 |                                                       
AAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCGATTTTGCTCGCCGTCGCCAGGAAAGTGGAA  >  minE/1486658‑1486778

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