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	97.8%	131.1 / NA	46	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (24/21);  minor base G (1/0);  total (25/21)
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.91e-01

GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132648‑2132925                                                                                                                                            |                                                                                                                                           gATCTGGTAGATGGTAAACAGCTGATATGTTTTGCTGGCGCGGTCGAATTTGTCTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                         <  2:219137/143‑1 (MQ=255) gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                         <  1:559630/143‑1 (MQ=255)  aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                        >  2:416178/1‑143 (MQ=255)      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTTATATGCAGGCTCTGCGTc                                                                                                                                    >  2:490558/1‑143 (MQ=255)        tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTCAAAGTGGCAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTccc                                                                                                                                  >  2:789277/1‑141 (MQ=255)                  aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGGATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                        >  1:554635/1‑143 (MQ=255)                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttttcg                                                                                                                   >  1:443953/1‑140 (MQ=255)                         atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                 >  1:451822/1‑143 (MQ=255)                                 tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGCGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                         <  1:285408/143‑1 (MQ=255)                                   cTGGCGCGGTCGAATTTGTCTTTGTACTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGGGATTGATATGCAGCGTCTGCGTCGCCGCTTATACTCTCCACGGTCAGCCgga                                                                                                       >  2:647638/1‑141 (MQ=255)                                   cTGGCGCGGTCGAATTTGTCTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                       <  1:493755/143‑1 (MQ=255)                                         cgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                 <  1:801744/143‑1 (MQ=255)                                                           tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                               <  1:500691/143‑1 (MQ=255)                                                            gCTTTCGGTATAAGTGGAAGTATCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTCCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATg                                                                              <  2:505078/143‑1 (MQ=255)                                                            gCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCCTCTCCTTTCCACGGTTAGCCGTACTCCACCACCCTCCGGAGGGAGCTg                                                                              >  2:325620/1‑143 (MQ=255)                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTAGCCTATACTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                             >  1:284832/1‑143 (MQ=255)                                                               cgCGGTAAAAGTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                           <  1:98645/142‑1 (MQ=255)                                                                  ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                        <  1:733812/143‑1 (MQ=255)                                                                   gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt                                                                       <  1:862005/143‑1 (MQ=255)                                                                    tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATACAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAGGAAGATGTCGTAATc                                                                      >  2:444176/1‑143 (MQ=255)                                                                        aGTGGAAGTACCCGGAGACTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGCTGTCGTAATCCtgt                                                                  >  1:594145/1‑143 (MQ=255)                                                                         gTGGAAGCACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                 <  1:159421/143‑1 (MQ=255)                                                                                 aCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGACGTCGTAATCCTGCGCCGCttcc                                                         >  2:603192/1‑141 (MQ=255)                                                                                   ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTCATCCTGTGCCGCtctccg                                                       >  1:12543/1‑141 (MQ=255)                                                                                   ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtctacg                                                       >  2:136404/1‑138 (MQ=255)                                                                                     ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCGTCCGCAGGAGGAAGTCGTAGCCGCGTGCCGCTTTCAGtt                                                     >  2:764389/1‑143 (MQ=255)                                                                                      gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCATCTTCTACTTTCCACGGTTCGCCGTACTCAACCACCCTCCGCTGGGAGAGGTCGTAATCCTGTCCCGCTGTCAGTTc                                                    >  1:838164/1‑143 (MQ=255)                                                                                         aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCCTTCAGTTCGGt                                                 >  2:227511/1‑143 (MQ=255)                                                                                         aaTTTTTTCGGTGACGAGGTGATCTGCAATTCACCGTGTATGTTATTTAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTTGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAGCCTGTGCCGGTTTCAGTTCGGt                                                 >  1:845351/1‑143 (MQ=255)                                                                                         aaTTTTTTCGGAGGCGAGTTGATCTGCATTTCTCCGTTATTGATATGCAGGCTCTGCGACGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGt                                                 <  2:785262/143‑1 (MQ=255)                                                                                           ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCTCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCCGTGCCGCTTTCCGTTCGGTAt                                               >  2:442352/1‑143 (MQ=255)                                                                                            tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                              <  1:416178/143‑1 (MQ=255)                                                                                               ttCGGTGACGAGGTAATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                           <  2:456606/143‑1 (MQ=255)                                                                                                  ggTGACGAGGTGATCTGCATTTCATCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGCCGCTGTTGTCATCCTTTGCTGCTTTCAGTTAGGTAttaaaaaa                                        >  1:13991/1‑136 (MQ=255)                                                                                                        gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCATGTGCCGCTTTCAGTTAGGTAGCGAACAAACActa                                  >  2:128237/1‑142 (MQ=255)                                                                                                          ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                <  1:636412/143‑1 (MQ=255)                                                                                                          ggTGATCTGCACTTCACCGTGATTGACATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                <  2:149432/143‑1 (MQ=255)                                                                                                          ggTGAGCTCCTTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                <  2:339310/143‑1 (MQ=255)                                                                                                               tCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                           <  1:226877/143‑1 (MQ=255)                                                                                                                  gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCATACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGc                        <  1:695662/143‑1 (MQ=255)                                                                                                                      gtacccGTGAGTGACATGCAGGCTCTCCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCATGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  2:605470/139‑1 (MQ=255)                                                                                                                       tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTCTCAGCTCGGTATCGCACAAACACACTTGCAcacgccacgc                   >  2:447028/1‑143 (MQ=255)                                                                                                                         aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCAGACAAACGATAACGGAGTTGTGTGCCGCTTTCAGTTCGGTATCGCACAACCACTCTTTGCACCGCCATGccc                 >  1:409783/1‑143 (MQ=255)                                                                                                                          ccgcgcttgctATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACTGTTCGCCGTCCTCCACCCCCATCCGCACGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCTGTATCGACCAAACACTCTTTGAAACGCCATGCCCa                <  2:743782/137‑1 (MQ=255)                                                                                                                           cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGGTGTCGGAATCCTGTGCCGCTTGCAGTTCGGTATCGAACAAACGCTCTTTGAAACGCCATGCCCAg               >  1:211782/1‑143 (MQ=255)                                                                                                                                     atGCCGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt     <  1:497953/143‑1 (MQ=255)                                                                                                                                       gCAGGCTCTGCGTCGCCGCGTCTACTTTCCACGGTCCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACTCTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg   <  1:9705/143‑1 (MQ=255)                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGt  >  2:675710/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                           GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132648‑2132925

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: ‑