Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 828,872 G→A 35.7% D142D (GAC→GAT narX ← sensory histidine kinase in two‑component regulatory system with NarL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE828,8720GA35.7% 54.2 / 47.1 81D142D (GAC→GATnarXsensory histidine kinase in two‑component regulatory system with NarL
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (21/8);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 3.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

AATACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGACAC  >  minE/828813‑828931
                                                           |                                                           
aaTACCGCCATTACCCGATGTACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:544883/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGGCCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:372156/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:959331/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:1031696/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:1043842/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:910841/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:906150/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:728700/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:712874/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:1582046/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:2324725/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:2001434/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:1025059/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:2149395/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:23023/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:2951080/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:456633/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:498404/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:1301016/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:1274542/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:1091536/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGTACCACCGTCTCGATGCGCATTTCCGTGGTGCGGTCaaa                                                        >  1:1198166/1‑65 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:590732/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:3410403/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:818402/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:2758094/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:2393523/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:1490870/69‑1 (MQ=255)
aaTACCGCCATTACCCGATGGACCAGCACCACTGTCTCAATGCGCATTTCCGTGGTGCGATCAAAACCa                                                    <  1:156445/69‑1 (MQ=255)
  tACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCAGa                                                  >  1:842809/1‑69 (MQ=255)
                            ccACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGcaacaa                        >  1:1261007/1‑69 (MQ=255)
                            ccACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGcaacaa                        >  1:2761735/1‑69 (MQ=255)
                                  tctcGATGCGCATTTCCGTGGTGCGGTCAAAACCa                                                    <  1:3385662/35‑1 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGc                  >  1:387007/3‑66 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:2189712/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:395204/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:390171/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:1661839/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:652/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:1237648/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:481498/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:3226876/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:2956849/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:1462493/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:280329/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:2214000/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:520867/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:1914636/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:907139/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:1968924/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:2488755/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:196919/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:2370810/3‑67 (MQ=255)
                                     caaTGCGCATTTCCGTGGTGCGATCAAAACCAGATACCAGCTGATCAAGCCCGGCAACAAACTGGCt                 >  1:1614491/3‑67 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGAc              >  1:1818002/1‑59 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACa             >  1:2907663/1‑60 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACa             >  1:1418012/1‑60 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACa             >  1:2340030/1‑60 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACa             >  1:2665768/1‑60 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACa             >  1:1884907/1‑60 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGac    >  1:169605/1‑69 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGac    >  1:649190/1‑69 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGac    >  1:1235025/1‑69 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGac    >  1:761797/1‑69 (MQ=255)
                                                tCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGac    >  1:2284958/1‑69 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:832267/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:640613/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:1019945/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:534655/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:462992/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:3406964/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:3395031/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:3327939/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:3021768/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:2945437/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:2795371/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:2661363/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:2535423/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:1492151/68‑1 (MQ=255)
                                                   gtggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:1704726/68‑1 (MQ=255)
                                                    tggtgCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGacac  <  1:2128502/67‑1 (MQ=255)
                                                           |                                                           
AATACCGCCATTACCCGATGGACCAGTACCACTGTCTCGATGCGCATTTCCGTGGTGCGGTCAAAACCAGATACCAGTTGATCAAGCCCGGCAACAAACTGGCTGACATCCGCTGACAC  >  minE/828813‑828931

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: