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	98.1%	154.5 / NA	52	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (31/20);  minor base G (0/1);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132647‑2132924                                                                                                                                             |                                                                                                                                          ggatgggggagagGGGAAGCAGCTGATATACTGTGCTGGCGGGGGGGACTGCGTCTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                         <  1:401092/134‑1 (MQ=255)    tCTGGTAGAGGGTAAACAGCTGATATTTCTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                      <  2:710197/143‑1 (MQ=255)      tcgTAGAGGGTAAACAGCTGATATTTCTTGCTGGCGCGGTCGCCTTTGTCTTCGTGCTTGCGGTAAGNGCGGAAGTACCCGGACAATTTCCTCGGCGACCAGGTGATCTGCATTCCACCGTCACTGACCTGCACGCTCCGCGt                                                                                                                                    >  1:575599/3‑143 (MQ=255)       ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                   <  2:676824/143‑1 (MQ=255)              ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                            <  1:672792/143‑1 (MQ=255)                   aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGTTCTGCATTTCACCGTGCTTGATATGCAGGCTCTACGTCGActcgcctagt                                                                                                                       >  2:84744/1‑136 (MQ=255)                   aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATCTGCCGGCTCTGCGTCGCCTCGTCTGCt                                                                                                                       >  2:746880/1‑143 (MQ=255)                   aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACt                                                                                                                       >  2:708153/1‑143 (MQ=255)                    catgtgctatTTTTTGCTGGCGCGGGCGAATTTGGCTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                      <  1:288458/136‑1 (MQ=255)                        tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                  >  2:396632/1‑143 (MQ=255)                             ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                             >  2:343253/1‑143 (MQ=255)                             ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                             <  1:889712/143‑1 (MQ=255)                              tttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTc                                                                                                            >  2:491030/1‑143 (MQ=255)                                           ggTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                               >  2:778488/1‑143 (MQ=255)                                             tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGCATTTTTTCGGTGCCGAGGTGAGCTGCATTTCACCGTGATTGATATGCAGGCGCTGAGTCGCCTCTTCTACTGTCCACGGTTCGCCCTACTctacgac                                                                                             >  1:736078/1‑143 (MQ=255)                                                         ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGATGTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                 <  1:604035/143‑1 (MQ=255)                                                           gtgCTCGCGGTAAAAGGGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                               <  2:895816/143‑1 (MQ=255)                                                                tgtggtaggaGTGGAAGTCTCCGGATGATTTGTTCGGTGACGTGGTGATCTGCGTTTCACCGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                          <  1:594263/134‑1 (MQ=255)                                                                 gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGAAGGCTCTGCGTCGCGTCTTATACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGt                                                                         >  1:212770/1‑143 (MQ=255)                                                                 gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCATGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGGACTCCCCCCACATCCGCAGGAAGCTGTCGt                                                                         >  1:326805/1‑143 (MQ=255)                                                                   ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCATCTTCTACTTTCCACGGTTCGCAGTACTCCACAAACATCCGCAAGAAGATGtagaaa                                                                       >  1:287426/1‑138 (MQ=255)                                                                    gcagatggggAAGTGCCCGGCGGATTTTGTCGGTGGCGTGGTGCTCTGCATTTCAGCGGTATTGAGATGTAGGGTCTGCGTCGCCTCTTCGACTTTCCGCGGTTCGCCGTACTCCACCACCGTCCGCAGGAAGATGTCGTAAt                                                                      <  1:494635/135‑1 (MQ=255)                                                                    gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGGTCGTAAt                                                                      >  1:641619/1‑143 (MQ=255)                                                                     tAAAAGTGGAAGTACCCGGAGACTTTTTACGGTGACGAGGTGATCTGAATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTACTACTATCGACGGATCGCGGTACTGCACCACACCCCGCACGACGATGGCGCAATc                                                                     >  1:403005/1‑143 (MQ=255)                                                                     tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTACTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGGTGTCGTAATc                                                                     >  1:717246/1‑143 (MQ=255)                                                                         aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                 <  1:708153/143‑1 (MQ=255)                                                                          gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCCCTGCGTCGCCTCTTCTACTCCCCCCGGTCCGCCGTACTCCACCACCACCCGCCGGACGCCGCCGCAATCCtgtg                                                                >  2:941208/1‑143 (MQ=255)                                                                               aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                           <  1:343253/143‑1 (MQ=255)                                                                                   cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                       <  2:477912/143‑1 (MQ=255)                                                                                    ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGCTGTCGTCATCCTGTGCCGCTTTCAg                                                      >  2:636068/1‑143 (MQ=255)                                                                                       gttcaTTTTTTCGTTGACGAGGGGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                   <  1:579520/139‑1 (MQ=255)                                                                                       gagaATTTTTTCGGTGACGAGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                   <  1:809637/143‑1 (MQ=255)                                                                                             tttttCGGTGACGAGGTGATCTTCATTTCACCGTGATTGATATGCAGGCTCTGCGTTGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGGAGGAGGATGTCGTAATCCTGTACCGCTTTCCATTCGATATc                                             >  1:768188/1‑143 (MQ=255)                                                                                                           ggTGATCTGCATTGGAGCGGGCTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCAGCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                               <  2:837860/143‑1 (MQ=255)                                                                                                            gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTTTCGAACAAACAGTCttt                              >  1:1014193/1‑143 (MQ=255)                                                                                                             tGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTg                             >  2:249511/1‑143 (MQ=255)                                                                                                                tgtGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                          <  1:743243/141‑1 (MQ=255)                                                                                                                tCTGCATTTCACCGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCTTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                          <  1:819461/143‑1 (MQ=255)                                                                                                                 cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                <  2:393062/72‑1 (MQ=255)                                                                                                                 cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                >  1:393062/1‑72 (MQ=255)                                                                                                                 cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGCAGATGTCGTACTCCTGTGCCGCCCTCAGCTCGGTCCCGACCAACCCCCCtttcagac                         >  2:447702/1‑138 (MQ=255)                                                                                                                   gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCCCCCTACTATATCTGCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGc                       >  2:637899/1‑143 (MQ=255)                                                                                                                     aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                         <  2:811112/107‑1 (MQ=255)                                                                                                                     aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                         >  1:811112/1‑107 (MQ=255)                                                                                                                          aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGACACGCCATGccc                >  2:677370/1‑143 (MQ=255)                                                                                                                                attgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACACACACTCTTTGAAACGACATTCCCAGGTaa          >  1:922060/1‑143 (MQ=255)                                                                                                                                attgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGCAACGCCATGCCCCGGTaa          >  1:154610/1‑143 (MQ=255)                                                                                                                                  tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGACCAAACACTCTTTGAAACGCCATGCCCAGGTAAAg        >  2:650061/1‑143 (MQ=255)                                                                                                                                   gatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGa       >  1:184838/1‑143 (MQ=255)                                                                                                                                       tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTCTGCCGCTTTCAGTTCGGTATCGAACAAACACCCTTTGAAACGCCCAGCCCAGGCaacgccccg   >  1:183638/1‑136 (MQ=255)                                                                                                                                       tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg   >  2:561433/1‑143 (MQ=255)                                                                                                                                        gCAGGCGCTGCGTCGCCTCTGCCACTTTCCTCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGGAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg  <  1:396632/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132647‑2132924

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