Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,344,624 G→A 63.0% E430E (GAG→GAA glnG → fused DNA‑binding response regulator in two‑component regulatory system with GlnL, nitrogen regulator I (NRI)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,344,6240GA63.0% 43.8 / 50.5 81E430E (GAG→GAAglnGfused DNA‑binding response regulator in two‑component regulatory system with GlnL, nitrogen regulator I (NRI)
Reads supporting (aligned to +/- strand):  ref base G (20/10);  new base A (23/28);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 6.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

CAGAGCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGAAGCGGCGCGGCTA  >  minE/2344567‑2344687
                                                         |                                                               
cagaGCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACt                                                      >  1:2496760/1‑69 (MQ=255)
cagaGCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACt                                                      >  1:3268410/1‑69 (MQ=255)
cagaGCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACt                                                      >  1:375088/1‑69 (MQ=255)
cagaGCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACt                                                      >  1:2611305/1‑69 (MQ=255)
 agagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:1620789/1‑69 (MQ=255)
 agagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:3013797/1‑69 (MQ=255)
 agagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:1823943/1‑69 (MQ=255)
 agagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:1664620/1‑69 (MQ=255)
 agagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:1620628/1‑69 (MQ=255)
 agagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:630583/1‑69 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:2613523/1‑68 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:1023943/1‑68 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:3226044/1‑68 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:444537/1‑68 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:536893/1‑68 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:1660872/1‑68 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:1236909/1‑68 (MQ=255)
  gagCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTg                                                     >  1:13786/1‑68 (MQ=255)
   agCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTgac                                                   <  1:3305352/69‑1 (MQ=255)
   agCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTgac                                                   <  1:1866022/69‑1 (MQ=255)
    gcgcTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTgac                                                   <  1:463234/68‑1 (MQ=255)
              ccGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGc                                        >  1:2200729/1‑69 (MQ=255)
                   catcaAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGCGACAt                                   >  1:1560341/1‑69 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCTACATACGCAg                             <  1:1630692/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:2354376/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:2748172/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:969726/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:2885643/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:2973585/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:3140000/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:3317342/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:367430/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:367765/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:380155/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:925064/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1008579/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:115513/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:2338347/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1858035/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1668475/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1869967/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1459121/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1421529/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1346274/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:22032/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:231259/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1134435/66‑1 (MQ=255)
                         aaccTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:2318412/66‑1 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:721587/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:67402/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:1800669/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:552629/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:1227469/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:73675/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:909930/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:943905/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:173355/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:1787032/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:110384/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:2401722/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:2714367/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:2718234/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:2814253/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:2078155/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:323963/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:3010608/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:2004081/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:3192603/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:3201138/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGc                                        >  1:2579592/3‑57 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:2237479/66‑1 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:49390/66‑1 (MQ=255)
                          accTGCTTTCCGAAGCACAGCCAGAGCTGGAACGGACGTTACTGACCACCGCGCTGCGACATACGCAg                             <  1:1613791/66‑1 (MQ=255)
                          accTGCTTACCGAAGCACAGCCAGAGCTGGAACGGACGTTACTg                                                     >  1:485720/3‑44 (MQ=37)
                                       aGCGCAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGa               <  1:836215/69‑1 (MQ=255)
                                       aGCGCAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGa               <  1:2898563/69‑1 (MQ=255)
                                       aGCGCAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGa               <  1:276397/69‑1 (MQ=255)
                                        gcgcAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGa               <  1:2559925/68‑1 (MQ=255)
                                                    cTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGAAGCGGCGCGGCTa  <  1:202313/69‑1 (MQ=255)
                                                    cTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGAAGCGGCGCGGCTa  <  1:3371763/69‑1 (MQ=255)
                                                    cTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGAAGCGGCGCGGCTa  <  1:1756436/69‑1 (MQ=255)
                                                         |                                                               
CAGAGCGCTGCGTTCCGGTCATCAAAATCTGCTTTCCGAAGCGCAGCCAGAGCTGGAGCGGACGTTACTGACGACCGCGTTGCGACATACGCAGGGGCATAAACAGGAAGCGGCGCGGCTA  >  minE/2344567‑2344687

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: