Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,486,500 G→A 38.0% A271V (GCC→GTC)  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,5000GA38.0% 47.7 / 52.7 84A271V (GCC→GTC) nuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base G (30/22);  new base A (12/20);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTACC  >  minE/1486439‑1486562
                                                             |                                                              
aTACGGGCTGGACATCATGCTGACCTCGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:162840/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:3042800/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:44800/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:2019826/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:2118287/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:1602824/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:1542899/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:2811818/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:3141069/69‑1 (MQ=255)
aTACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:2999054/69‑1 (MQ=255)
 tACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:3374419/68‑1 (MQ=255)
 tACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:73291/68‑1 (MQ=255)
 tACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:1816595/68‑1 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCtt                                           >  1:1637704/1‑68 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          >  1:3350393/1‑69 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          >  1:329960/1‑69 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          >  1:1973631/1‑69 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          >  1:207274/1‑69 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          >  1:508934/1‑69 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          >  1:57734/1‑69 (MQ=255)
               cAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          >  1:1374227/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCGTACCGCTTcc                                         >  1:2545625/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:922561/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:811576/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:3229645/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:2758763/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:260481/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:2401529/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:1379639/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:660853/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTcc                                         >  1:1701228/1‑69 (MQ=255)
                aGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          <  1:2192368/68‑1 (MQ=255)
                      cccACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:2518080/2‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGa                                    >  1:2709375/1‑68 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:873524/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:2888443/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:670446/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:776909/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:2702029/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:378431/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:2475505/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:1407598/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:3347999/1‑69 (MQ=255)
                      aCCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGAc                                   >  1:2191737/1‑69 (MQ=255)
                        cACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACAt                                 <  1:1708889/69‑1 (MQ=255)
                        cACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACAt                                 <  1:2258235/69‑1 (MQ=255)
                                 cgcgCCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:535847/36‑1 (MQ=255)
                                  gcgcCGAGGCTGCTGAACAGATAACCGGCCAGGTa                                                         <  1:1074107/35‑1 (MQ=255)
                                          gcttctGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAgcgc               <  1:2341033/69‑1 (MQ=255)
                                          gctgctGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAgcgc               <  1:3093528/69‑1 (MQ=255)
                                          gctgctGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAgcgc               <  1:3171834/69‑1 (MQ=255)
                                           ctgctgAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCa              >  1:867730/1‑69 (MQ=255)
                                                gAACAGATAACCGGCCAGGTAAACCCCTACCGCTTc                                          <  1:1227652/36‑1 (MQ=255)
                                                 agcAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:820866/65‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCTCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:461560/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:746417/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:1764847/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:64387/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:630836/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:3028091/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:1854006/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:1883846/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:1932857/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:2089949/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:2153639/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:851913/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:2604918/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:2887152/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:1109322/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:1015571/67‑1 (MQ=255)
                                                 aaCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:936327/67‑1 (MQ=255)
                                                  aCAGATAACCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:1843572/66‑1 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTg                     >  1:853692/1‑50 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:744013/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:1736382/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:1763058/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:1154643/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:3393422/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:2299623/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:3227147/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:890582/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:3089751/1‑69 (MQ=255)
                                                       tAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTAcc  >  1:3088813/1‑69 (MQ=255)
                                                         aCCGACCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCa          <  1:2917448/59‑1 (MQ=255)
                                                             |                                                              
ATACGGGCTGGACATCAGGCTGACCACGCCGAACGCGCCGAGGCTGCTGAACAGATAACCGGCCAGGTAAACCCCTACCGCTTCCATCGACATCTCGCCGGTTTGCAGCGCAATCAGCGCTACC  >  minE/1486439‑1486562

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: