Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,132,787	A→C	I204S (ATC→AGC)	wcaA ←	putative glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,132,787	0	A	C	94.1%	90.1 / NA	34	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (17/15);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACC  >  NC_000913/2132647‑2132921                                                                                                                                             |                                                                                                                                       ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGCCTTTTTTCGGTGACGAGGTCCTCTGCATATCAACGTGAATAATAAGCAAGAta                                                                                                                                      >  1:32841/1‑142 (MQ=255) agATATGGTAGGGGGTAAACAGCTGATATGCGTTGCTGGCGCGGTCGAATTTGTCCTTGTGCTTGCTGTTAAAGTGGAAGTACCCGGACAATTTTTTCGGTGACGAGGAGATCGGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                      <  1:296650/142‑1 (MQ=255)  gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGCTCTGCATTTCACCGTGATTGATATGCAGGAGCt                                                                                                                                     >  1:658398/1‑143 (MQ=255)       ggTAGAGGTTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                <  1:84555/143‑1 (MQ=255)       ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGTTCTGCATTCCACCGTGATCGATATGCAGGCCCTGCGTc                                                                                                                                >  2:150323/1‑143 (MQ=255)        gTAGAGGGTAAACGGCGGATCGTGTTTGCTGGCGCGGTCGAATGTGTCTTTGTTCTTGCGGTAAAAGGGGAAGTTCTCGGTGAATTTTTTCGGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                               <  2:489867/143‑1 (MQ=255)        gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                               <  1:622407/143‑1 (MQ=255)            agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTTAACTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                           >  1:570175/1‑143 (MQ=255)              ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                         <  1:404192/143‑1 (MQ=255)                    cggcggagagtgtgtGCTGGCGCGGTCGAATTTGTCTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                   <  2:693998/129‑1 (MQ=255)                     aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTGTACttt                                                                                                                  >  2:546320/1‑143 (MQ=255)                                      ggCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTATTCGGTGACGAGGTGATCTGCATTTCACCGTGATGGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGtccc                                                                                                 >  1:206625/1‑140 (MQ=255)                                        cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGtcctcc                                                                                               >  2:278925/1‑143 (MQ=255)                                            gTCTAATTTGTCTTTGTGCTTTCTGGAAAAGTGCAAGTACCCGGAGAATTTTTTCGGTGACGATGTGATCTGCCTTTCACTGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                           <  2:291393/143‑1 (MQ=255)                                             tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTccaccac                                                                                          >  1:367799/1‑143 (MQ=255)                                               gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGACGTACTccaccacca                                                                                        >  2:96506/1‑143 (MQ=255)                                               gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAAGCTCTGCGTCGCCTCTTCTACTGTCCCCGGTTCGCCGTACTccaccacga                                                                                        >  2:101264/1‑141 (MQ=255)                                                    tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                   >  1:233489/1‑143 (MQ=255)                                                        tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGa                                                                               >  1:211205/1‑143 (MQ=255)                                                         ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGTTCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAggcc                                                                              >  2:663633/1‑141 (MQ=255)                                                                       aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                <  1:184276/143‑1 (MQ=255)                                                                         aggggACGTCCCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                              <  1:187253/140‑1 (MQ=255)                                                                           tgctaGTACTCGGAGAATTTTTTCGGTGACGATGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCTACTGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                            <  2:109651/139‑1 (MQ=255)                                                                                gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                       >  1:595965/1‑143 (MQ=255)                                                                                   cccGGAGAATTTTTTCGGTGACGAGTTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTGCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                    <  1:481625/143‑1 (MQ=255)                                                                                   cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCCTGCAGCTCTCGTAATCCCGcgccgcctcca                                                    >  2:719281/1‑138 (MQ=255)                                                                                         gaATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                              <  1:730447/143‑1 (MQ=255)                                                                                             tttgtgggTGAGGAGGTGATCTGCATTTCACCTGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                          <  2:19948/137‑1 (MQ=255)                                                                                                    gTGACGAGGTGATCTGCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCTACCACCCTCCGCCGGACGATGTTGTGCTCCTTTCCCGCTTTCAGTTCGGTATCGaacaaa                                   <  2:54079/143‑1 (MQ=255)                                                                                                    gTGACGAGGTGATCTGCATTTAACCGAGATTGCTATGCAGGCTCTGCGACGCCTCTTCTACTTTCAACGGTTCGCCGTACTGCACCACCATCCGCAGGAAGATGGCGAAATCCGATGCCGCTGTCAGTTCGGTATGGaacacc                                   >  2:381881/1‑141 (MQ=255)                                                                                                      gagcgggTGATCTGCATTTCACGGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCGTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                 <  2:382548/138‑1 (MQ=255)                                                                                                       aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTGGCGAACAAacac                                >  1:399636/1‑143 (MQ=255)                                                                                                         gAGGTGATCTGCATTTCACAGTGATTGATATGCAGGCTCTGCGTCGCCTCATCTACTTTCCACGGTTCACCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCGGCTTTCAATGCGGAAACGaacaaacggac                              >  2:249062/1‑139 (MQ=255)                                                                                                                                     tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGCAAAGatc  >  1:193338/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                       GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACC  >  NC_000913/2132647‑2132921

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