Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,417,873 T→C 46.4% intergenic (+154/‑179) proL → / → narP tRNA‑Pro/DNA‑binding response regulator in two‑component regulatory system with NarQ or NarX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,417,8730TC46.4% 5.0 / 57.7 69intergenic (+154/‑179)proL/narPtRNA‑Pro/DNA‑binding response regulator in two‑component regulatory system with NarQ or NarX
Reads supporting (aligned to +/- strand):  ref base T (21/16);  new base C (19/13);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAG  >  minE/1417806‑1417936
                                                                   |                                                               
tGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCtt                                                                <  1:182712/69‑1 (MQ=255)
tGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCtt                                                                <  1:961976/69‑1 (MQ=255)
tGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCtt                                                                <  1:255666/69‑1 (MQ=255)
tGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCtt                                                                <  1:1507783/69‑1 (MQ=255)
    tATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTtg                                                               >  1:3040197/1‑66 (MQ=255)
    tATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTtg                                                               >  1:3036001/1‑66 (MQ=255)
            gTTAACTGTTTAAAGAAATGCACAAGTTTTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     <  1:2651289/68‑1 (MQ=255)
            gTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     <  1:997512/68‑1 (MQ=255)
            gTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     <  1:378965/68‑1 (MQ=255)
            gTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     <  1:1543678/68‑1 (MQ=255)
            gTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     <  1:149982/68‑1 (MQ=255)
                   gTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCaaa                                             >  1:2388106/1‑69 (MQ=255)
                   gTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCaaa                                             >  1:241349/1‑69 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGc                                                      >  1:1425833/1‑59 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:2327422/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:1454744/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:2644543/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:1206197/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:444086/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:2497470/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:1659325/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:744757/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:1421690/1‑60 (MQ=255)
                    tttAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCt                                                     >  1:1447189/1‑60 (MQ=255)
                             aTGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCtt                                   >  1:3170337/1‑69 (MQ=255)
                             aTGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCtt                                   >  1:3088228/1‑69 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:2883561/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1062048/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:2381909/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:2767082/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:2570565/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:2280032/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1894688/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:448439/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1742784/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1704250/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1538498/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:710039/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1469810/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1231192/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1156899/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtata                               >  1:1130388/1‑68 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTtat                                >  1:2537528/1‑67 (MQ=255)
                                  caAGTATTGTGATTGATTTTTTAGTTGTTTTTACTGATGAGaa                                                        >  1:2598475/1‑43 (MQ=38)
                                  caAGTATTGTGATTGATTTTTTAGTTCTTTTTCCTGATGAGAAGCTGATGCAAAATTc                                         >  1:1295429/1‑58 (MQ=255)
                                       aTTGTGATTGATTTTTTAGTTGTTTTTCTTCATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGaaa                         <  1:2299206/69‑1 (MQ=255)
                                          gtgATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAatg                      >  1:2928899/1‑69 (MQ=255)
                                          gtgATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAatg                      >  1:1766278/1‑69 (MQ=255)
                                          gtgATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAat                       <  1:2778686/68‑1 (MQ=255)
                                          gtgATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAat                       <  1:953011/68‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAg             >  1:2664576/1‑68 (MQ=255)
                                                    ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAg             >  1:891167/1‑68 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:1077149/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:838108/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:1526760/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:708440/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:637016/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:1769365/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:3238542/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:2326614/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:2953859/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:2863493/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:2859707/60‑1 (MQ=255)
                                                    ttttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:282662/60‑1 (MQ=255)
                                                     tttAGTTGTTTTTCCTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAGatga                     <  1:2313789/59‑1 (MQ=255)
                                                           tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGAc     <  1:2148687/69‑1 (MQ=255)
                                                           tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGAc     <  1:3220619/69‑1 (MQ=255)
                                                           tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGAc     <  1:312482/69‑1 (MQ=255)
                                                           tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGAc     <  1:1363575/69‑1 (MQ=255)
                                                           tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGAc     <  1:833865/69‑1 (MQ=255)
                                                           tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGAc     <  1:2880603/69‑1 (MQ=255)
                                                               tttCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAg  <  1:2462104/68‑1 (MQ=255)
                                                               tttCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGACAAAGCGAACGACAAg  <  1:483228/68‑1 (MQ=255)
                                                                   |                                                               
TGTTTATTTTCTGTTAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAG  >  minE/1417806‑1417936

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: