Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,486,694 G→A 66.6% S206S (AGC→AGT 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,6940GA66.6% 68.9 / 50.6 94S206S (AGC→AGTnuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base G (8/23);  new base A (23/39);  total (32/62)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CGCCATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCG  >  minE/1486628‑1486748
                                                                  |                                                      
cgcCATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAAtcgcttcgc                                                       <  1:3274625/68‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:2269758/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:1778785/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:1997835/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:572873/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:2534684/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:1474035/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:783160/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:1179071/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:784708/69‑1 (MQ=255)
  ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:1090080/69‑1 (MQ=255)
   cATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                    <  1:2446702/68‑1 (MQ=255)
       aGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTcaccc                                               >  1:637840/1‑69 (MQ=255)
       aGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTcaccc                                               >  1:474708/1‑69 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:2402086/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:2391279/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:2567048/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:2187304/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:2350852/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:2932892/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:3052767/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:2356199/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:1796666/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:3274182/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:1177886/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:1373866/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:654045/68‑1 (MQ=255)
         gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:1528746/68‑1 (MQ=255)
           tACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCCCCACGCGAATCGCTTCGCTGTCACCCACCg                                           >  1:1630808/1‑69 (MQ=255)
                aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                      <  1:2203364/69‑1 (MQ=255)
                aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                      <  1:2366285/69‑1 (MQ=255)
                aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                      <  1:2653188/69‑1 (MQ=255)
                aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                      <  1:376933/69‑1 (MQ=255)
                  gatgatGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                              <  1:1983953/59‑1 (MQ=39)
                  gatgatGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTa                                    <  1:894640/69‑1 (MQ=255)
                              aaaGGCGATAATCGCCAGCACCACGCGAATCGCTTCTCTGTCACCCACCGGTGCGTAGAGGAACAGACg                        >  1:2617519/1‑69 (MQ=255)
                              aaaGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACg                        <  1:15284/69‑1 (MQ=255)
                              aaaGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACg                        <  1:712642/69‑1 (MQ=255)
                               aaGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACAGGTGCGTAGAGGAACAGACg                        <  1:69177/68‑1 (MQ=255)
                                  gCGATATTCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:438766/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAAtcgcttcact                                                      >  1:2881043/1‑35 (MQ=25)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGt                                                    >  1:1757429/1‑37 (MQ=37)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGGACAAACGCATc                    >  1:981785/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCAt                     >  1:2325328/1‑68 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCAt                     >  1:662874/1‑68 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:3237005/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:738223/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:1491609/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:1497475/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:505729/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:459161/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:173632/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:2400080/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:3208336/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:297999/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:29555/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:1826388/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:1937003/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:2820034/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:2667625/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:231331/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:2522969/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                    >  1:2607136/1‑69 (MQ=255)
                                  gCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAgg                                <  1:1740785/57‑1 (MQ=255)
                                  gCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAgg                                <  1:547344/57‑1 (MQ=255)
                                    gATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcac                  <  1:1377205/69‑1 (MQ=255)
                                       aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac               >  1:814620/1‑69 (MQ=255)
                                       aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac               >  1:787501/1‑69 (MQ=255)
                                       aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac               >  1:2170444/1‑69 (MQ=255)
                                       aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac               >  1:2626965/1‑69 (MQ=255)
                                             cAGCACACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGaa         <  1:1964516/68‑1 (MQ=255)
                                                 accacGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATa     >  1:3127515/1‑69 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:247200/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:1496257/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:2435626/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:2302196/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:2424066/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:1289592/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:115691/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:110586/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:962193/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:1081768/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:600953/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:1561000/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:1592695/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:279011/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:2216542/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:17444/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:3273639/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:3239046/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:3207294/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:3068188/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:1812172/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:2774316/69‑1 (MQ=255)
                                                    acGCGAATCGCTTCACTGTCACCCACCGGTGCGAAAAGGATCAAACGCATCACCACGCCGAAGATAGCg  <  1:1732979/69‑1 (MQ=255)
                                                    acGCGAATCGCATCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATCACCACGCCGAAGATAGCg  <  1:3426052/69‑1 (MQ=255)
                                                                  |                                                      
CGCCATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATAGCG  >  minE/1486628‑1486748

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: