Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,137,778	A→G	10.0%	intergenic (‑535/‑124)	wza ← / → yegH	colanic acid export protein; outer membrane auxillary lipoprotein/inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,137,778	0	A	G	10.0%	118.9 / 4.6	40	intergenic (‑535/‑124)	wza/yegH	colanic acid export protein; outer membrane auxillary lipoprotein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (22/14);  new base G (2/2);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTAC  >  NC_000913/2137586‑2138047                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              gaggaAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATa                                                                                                                                                                                                                                                                               >  2:359154/1‑193 (MQ=255)     aaGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAAtgt                                                                                                                                                                                                                                                                           <  4:359149/193‑1 (MQ=255)                                               gTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAAt                                                                                                                                                                                                                                                                             >  5:243038/1‑149 (MQ=255)                                                      gCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAg                                                                                                                                                                                                                                                                        <  6:326365/147‑1 (MQ=255)                                                      gCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAg                                                                                                                                                                                                                                                                        >  5:326365/1‑147 (MQ=255)                                                          tttttAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTTTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGttatt                                                                                                                                                                                                                                                                  <  6:130898/149‑1 (MQ=255)                                                          tttttAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATa                                                                                                                                                                                                                                                                               <  8:263528/136‑1 (MQ=255)                                                          tttttAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATa                                                                                                                                                                                                                                                                               >  7:263528/1‑136 (MQ=255)                                                           ttttAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGttatta                                                                                                                                                                                                                                                                 >  7:142682/1‑149 (MQ=255)                                                                   aTTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCaaa                                                                                                                                                                                                                                                         <  8:75299/149‑1 (MQ=255)                                                                           aaGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATga                                                                                                                                                                                                                                                 <  8:142682/149‑1 (MQ=255)                                                                                       atatGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGtgc                                                                                                                                                                                                                                     >  8:34470/1‑149 (MQ=255)                                                                                            gTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCaa                                                                                                                                                                                                                                >  7:131392/1‑149 (MQ=255)                                                                                                ttAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACAtt                                                                                                                                                                                                                            <  6:201335/149‑1 (MQ=255)                                                                                                                    aCAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCa                                                                                                                                                                                                                         >  5:2109/1‑132 (MQ=255)                                                                                                                    aCAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCa                                                                                                                                                                                                                         <  6:2109/132‑1 (MQ=255)                                                                                                                    aCAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTAt                                                                                                                                                                                                        <  6:243038/149‑1 (MQ=255)                                                                                                                        tatTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCAt                                                                                                                                                                                                    >  8:367015/1‑149 (MQ=255)                                                                                                                          tttAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTATTTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGa                                                                                                                                                                                                  >  7:305864/1‑149 (MQ=255)                                                                                                                          tttAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATGATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGa                                                                                                                                                                                                  >  8:114748/1‑149 (MQ=255)                                                                                                                          tttAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATGATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGa                                                                                                                                                                                                  >  6:326043/1‑149 (MQ=255)                                                                                                                                  aGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACATGGTc                                                             >  2:287925/1‑274 (MQ=255)                                                                                                                                   gtgtCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGAt                                                                                                                                                                                                            >  5:140418/1‑130 (MQ=255)                                                                                                                                   gtgtCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGAt                                                                                                                                                                                                            <  6:140418/130‑1 (MQ=255)                                                                                                                                        tcCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACAt                                                                                                                                                                                    >  5:181166/1‑149 (MQ=255)                                                                                                                                           acaTGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTg                                                                                                                                                                                 >  5:394278/1‑149 (MQ=255)                                                                                                                                                 gATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTg                                                                                                                                                                                 <  6:273704/143‑1 (MQ=255)                                                                                                                                                 gATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTg                                                                                                                                                                                 >  5:273704/1‑143 (MQ=255)                                                                                                                                                         ttAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATa                                                                                                                                                                   <  8:131392/149‑1 (MQ=255)                                                                                                                                                          caaaTAATGTTTTATTATTACCACTTTTGTTTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATATAACGAATAg                                                                                                                                                                  <  8:305864/148‑1 (MQ=255)                                                                                                                                                               aatGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGtttt                                                                                                                                                             >  7:406674/1‑149 (MQ=255)                                                                                                                                                                 tGTTTTATTATTACCACTTTTAATTCAGGGATGATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGt                                                                                                                                                           <  5:326043/149‑1 (MQ=255)                                                                                                                                                                 tGTTTTATTATTACCACTTTTAATTCAGGGATGATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGt                                                                                                                                                           <  7:114748/149‑1 (MQ=255)                                                                                                                                                                       attattaCCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGtacttac                                                                                                                                                     <  7:367015/149‑1 (MQ=255)                                                                                                                                                                           ttaCCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTatg                                                                                                                                                 >  5:93829/1‑149 (MQ=255)                                                                                                                                                                             aCCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTatgga                                                                                                                                               >  7:376761/1‑149 (MQ=255)                                                                                                                                                                                aCTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTatggaatg                                                                                                                                            >  7:63619/1‑149 (MQ=255)                                                                                                                                                                                  ttttAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTg                                                                                                                                                            <  6:162633/131‑1 (MQ=255)                                                                                                                                                                                  ttttAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTg                                                                                                                                                            >  5:162633/1‑131 (MQ=255)                                                                                                                                                                                  ttttAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTatggaatgga                                                                                                                                          >  7:436518/1‑149 (MQ=255)                                                                                                                                                                                        ttCAGGGATAATGTGAGGTTATTAACCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCAATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCg                                                                                                                                    >  6:453264/1‑149 (MQ=255)                                                                                                                                                                                            gggATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGTGTTAATCACGCTGTTTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGCGACCGCGCACGGGTTAc  >  1:18427/1‑275 (MQ=255)                                                                                                                                                                                             ggATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCg                                                                                                                               >  8:269765/1‑149 (MQ=255)                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              GAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATAATGTGAGGTTATTACCCTCAAATAATATGAGATAAATAGTGCTGCAACATTGCATTTTTGCCCCGATTTATCCATGATCGAATTGTGACATTTGTCATACAACGAATAGGTTTTGTACTTACTATGGAATGGATTGCCGATCCGTCTATCTGGGCCGGGTTAATCACGCTGATTGTGATCGAACTGGTCCTCGGCATTGATAACCTGGTCTTTATTGCCATCCTCGCCGAAAAACTACCGCCGAAGCAGCGTGACCGCGCACGGGTTAC  >  NC_000913/2137586‑2138047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑