Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,546,268 A→G 39.5% S93S (TCT→TCC rfaF ← ADP‑heptose:LPS heptosyltransferase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,546,2680AG39.5% 36.0 / 55.2 81S93S (TCT→TCCrfaFADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base A (33/16);  new base G (20/12);  total (53/28)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

CATCTCGCCGCGCCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGACCCAGTT  >  minE/2546211‑2546335
                                                         |                                                                   
cATCTCGCCGCGCCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGa                                                          >  1:3010507/1‑69 (MQ=255)
cATCTCGCCGCGCCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGa                                                          <  1:2672583/69‑1 (MQ=255)
cATCTCGCCGCGCCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGa                                                          >  1:3337014/1‑69 (MQ=255)
         gcgcCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTaa                                                 <  1:3280237/69‑1 (MQ=255)
          cgcCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTaa                                                 <  1:566884/68‑1 (MQ=255)
          cgcCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTaa                                                 <  1:1359888/68‑1 (MQ=255)
          cgcCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTaa                                                 <  1:764030/68‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:2495443/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:2067576/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:2282985/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:2304610/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:444313/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:1672431/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:364733/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:3141597/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:1438866/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:2963775/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:1319167/69‑1 (MQ=255)
                      cGATGAGGAATACCTGCGAAGAAAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCgg                                    <  1:3198593/69‑1 (MQ=255)
                      cGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCgg                                    >  1:2260415/1‑69 (MQ=255)
                      cGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCgg                                    >  1:1230937/1‑69 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGt                                                   >  1:1005668/1‑36 (MQ=37)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:725203/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:1330160/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:284598/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:1348814/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:2791755/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:833241/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:88031/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:926960/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:158937/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:1619191/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:699127/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:1917978/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:2203896/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:1156576/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:1974624/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:533359/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:2054081/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:2044317/1‑54 (MQ=255)
                                        aagaaAGGCACTAATGCGGATTTGAAGGAGTTAGGTAAGACGTAGGCGCGGTCg                                 >  1:1999762/1‑54 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGTTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2684642/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:1070667/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:3028320/59‑1 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:3216906/59‑1 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:3402219/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:3429566/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:359074/59‑1 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:624099/59‑1 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:496065/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:499897/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:1646361/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2123785/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:193456/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2380837/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:1735420/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2649284/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2703432/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2804313/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:1716011/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2871446/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:288970/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2913724/1‑59 (MQ=255)
                                                 aCTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:2949247/1‑59 (MQ=255)
                                                   tAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGgcg                                       <  1:1162615/37‑1 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:3044803/56‑1 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:1582359/1‑56 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:1235809/1‑56 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:1778541/56‑1 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:957984/56‑1 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:3301110/56‑1 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   <  1:325003/56‑1 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACg                   >  1:3129661/1‑56 (MQ=255)
                                                    aaTGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGAccc      >  1:3235742/1‑69 (MQ=255)
                                                      tGCAGATTTGAATGAGTTGGGTAAGACGTATGAGCGGTCGTAGCGCTTTTCACg                   <  1:982532/54‑1 (MQ=255)
                                                        cAGATTTGAAGGAGTTGGGTAAGACGTAGTCGCGGTCGTAgcgc                           >  1:2139303/1‑44 (MQ=255)
                                                        cAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGACCCAGtt  >  1:1938144/1‑69 (MQ=255)
                                                        cAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGACCCAGtt  >  1:2366402/1‑69 (MQ=255)
                                                        cAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGACCCAGtt  >  1:2556052/1‑69 (MQ=255)
                                                        cAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGACCCAGtt  >  1:2866412/1‑69 (MQ=255)
                                                        cAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGACCCAGtt  >  1:3242670/1‑69 (MQ=255)
                                                         |                                                                   
CATCTCGCCGCGCCAGCCGGTGCGATGAGGAATACCCGCGAAGAAAGGCACTAATGCAGATTTGAAGGAGTTGGGTAAGACGTAGGCGCGGTCGTAGCGCTTTTCACGCAGGCTATGACCCAGTT  >  minE/2546211‑2546335

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: