Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,486,679 C→T 54.9% V211V (GTG→GTA 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,6790CT54.9% 9.5 / 57.4 71V211V (GTG→GTAnuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base C (8/24);  new base T (24/15);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 3.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGTCTGGCTCAGCGCCATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATA  >  minE/1486614‑1486745
                                                                 |                                                                  
ttGGTCTGGCTCAGCGCCATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGcaccac                                                                  >  1:2942329/1‑68 (MQ=255)
              cgcCATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAAtcgcttcgc                                                    <  1:3274625/68‑1 (MQ=255)
                ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                 <  1:1090080/69‑1 (MQ=255)
                ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                 <  1:2269758/69‑1 (MQ=255)
                ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                 <  1:1997835/69‑1 (MQ=255)
                ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                 <  1:2534684/69‑1 (MQ=255)
                ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                 <  1:1778785/69‑1 (MQ=255)
                ccATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                 <  1:572873/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)
                 cATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGt                                                 <  1:2446702/68‑1 (MQ=255)
                     aGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTcaccc                                            >  1:474708/1‑69 (MQ=255)
                     aGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTcaccc                                            >  1:637840/1‑69 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:1177886/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:2350852/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:2356199/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:2391279/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:2402086/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:2567048/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:2932892/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:3052767/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:3274182/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:654045/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:2187304/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:1528746/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:1796666/68‑1 (MQ=255)
                       gTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTcaccca                                           <  1:1373866/68‑1 (MQ=255)
                         tACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCCCCACGCGAATCGCTTCGCTGTCACCCACCg                                        >  1:1630808/1‑69 (MQ=255)
                              aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                   <  1:376933/69‑1 (MQ=255)
                              aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                   <  1:2653188/69‑1 (MQ=255)
                              aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                   <  1:2366285/69‑1 (MQ=255)
                              aaGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCg                                   <  1:2203364/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: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)
                                                gCGATAATCGCCAGCACTACGCGAAtcgcttcact                                                   >  1:2881043/1‑35 (MQ=25)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGt                                                 >  1:1757429/1‑37 (MQ=37)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGGACAAACGCATc                 >  1:981785/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCAt                  >  1:662874/1‑68 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCAt                  >  1:2325328/1‑68 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:297999/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:231331/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:2400080/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:738223/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:2522969/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:1937003/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:1491609/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:1497475/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:1826388/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:505729/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:2607136/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:459161/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:2667625/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:2820034/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:173632/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:3237005/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:3208336/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACTACGCGAATCGCTTCACTGTCACCCACCGGTGCGTAAAGGAACAAACGCATc                 >  1:29555/1‑69 (MQ=255)
                                                gCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAgg                             <  1:547344/57‑1 (MQ=255)
                                                gCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAgg                             <  1:1740785/57‑1 (MQ=255)
                                                  gATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcac               <  1:1377205/69‑1 (MQ=255)
                                                     aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac            >  1:2626965/1‑69 (MQ=255)
                                                     aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac            >  1:787501/1‑69 (MQ=255)
                                                     aaTCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATcaccac            >  1:814620/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)
                                                                 |                                                                  
TTGGTCTGGCTCAGCGCCATCAGGTTACCGAAGATGATGGAGGCAAAGGCGATAATCGCCAGCACCACGCGAATCGCTTCGCTGTCACCCACCGGTGCGTAGAGGAACAGACGCATCACCACACCGAAGATA  >  minE/1486614‑1486745

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