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.0%	145.7 / NA	51	I204S (ATC→AGC)	wcaA	putative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (26/24);  minor base G (0/1);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132650‑2132926                                                                                                                                          |                                                                                                                                            tCTGGTGGATGTTATACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGTGGTAAAAGTGGAAGTACCCTGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  1:777177/143‑1 (MQ=255) tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  2:442455/143‑1 (MQ=255)    ggTAGAGTGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                     <  1:777805/143‑1 (MQ=255)     gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                    >  1:791842/1‑143 (MQ=255)      tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  1:936459/143‑1 (MQ=255)              aaaCAGTTGATCTTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           <  2:698913/143‑1 (MQ=255)                 cAGCTGATGTTTTCGGCTGGCTCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATCTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                        <  1:332604/143‑1 (MQ=255)                  aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGCGAATTTTTTCGGTGACGAGGTGATCTGCATTTAACCGTGATCGATATACAGGCACTGCTTCGCCTCTTCTACttt                                                                                                                       >  1:61095/1‑143 (MQ=255)                   gcggATATTGTTGGGTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                      <  1:909532/140‑1 (MQ=255)                     tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    >  2:954138/1‑143 (MQ=255)                     tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    >  2:239784/1‑143 (MQ=255)                     tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    >  2:723969/1‑143 (MQ=255)                       gtgttgtttGCTGGCGCGGTCGAATTTGTCTTTGGTCTGGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                  <  1:919910/137‑1 (MQ=255)                           tttttGCTGGCGCGTTCGAATTTGTCTTTGTGGTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTc                                                                                                              <  1:310913/143‑1 (MQ=255)                            ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCCGGCCCTGCGTCGGACCATCCACTTTCCCCGGTTCg                                                                                                             >  1:386540/1‑143 (MQ=255)                                gCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTAATATGCAGGCTCTGCGTCGCCTCGTCTACTTTTCACGGTTCGCCGt                                                                                                         >  2:451375/1‑143 (MQ=255)                                 cTGGCGCGGTCGAAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                        <  2:400450/143‑1 (MQ=255)                                  tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATATGACTTTCAGCTTGATGGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTGCACGGGTCGCCGCAc                                                                                                       >  1:859442/1‑143 (MQ=255)                                        ggTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGCTATGCAGGCTGTGCGTCGCCTGTTCTACTTTCCCGGGTTCGCCGTTCTgcacc                                                                                                 >  2:243420/1‑143 (MQ=255)                                                 tttctccgtccctGCGCTAAACGTGCAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                        <  2:64995/131‑1 (MQ=255)                                                  gTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                       <  2:190200/143‑1 (MQ=255)                                                                ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGGTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACAATCCGCAGGCAGATTTCGTat                                                                         >  2:47328/1‑142 (MQ=255)                                                                  tAAAAGTGGAAGTACTGGGAGAATGTTTTTGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                       <  1:776248/143‑1 (MQ=255)                                                                  tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGGTATTGATATGCAGGCTCTGCGGCGCCTCATCTACTTTCCACGATTAGCCGGACTAGACCCCCACCCGCAGGCAGATGCCGtcacc                                                                       >  1:46329/1‑139 (MQ=255)                                                                      aGTGGAGGTACCCGGAGAATCTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                   <  2:328182/143‑1 (MQ=255)                                                                                  cggacgaTTTTTCCGGTGCCGTGGTCATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTATTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGt                                                       <  1:456996/137‑1 (MQ=255)                                                                                   ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:58325/143‑1 (MQ=255)                                                                                      gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                   >  2:973974/1‑143 (MQ=255)                                                                                      gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAGCCTGTGCCGCTTTCAGTGCgg                                                   >  1:266228/1‑143 (MQ=255)                                                                                       gaTTTTTTCGGTGACGGGGTCGTCGGCATTTCGCCGTGGTTGATATGCAGGGCCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTTCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGt                                                  <  1:284849/142‑1 (MQ=255)                                                                                       aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCTCCACCCTCCGCCGTAACCTGTCGTCATCCTGTGCCGACTTCCGTTCGGt                                                  >  1:984111/1‑143 (MQ=255)                                                                                          tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                               <  2:666632/143‑1 (MQ=255)                                                                                             tgcGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                            <  2:497286/141‑1 (MQ=255)                                                                                               cGGTGACGACGTTTTCTTCTTTTCTCCTTTATTTATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                                          <  2:216918/143‑1 (MQ=255)                                                                                                 gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGACGTCGTAATCCTCTGCCGCTTCCAGTTCCGTACCGaacaac                                        >  1:627342/1‑142 (MQ=255)                                                                                                    aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac                                     >  2:230218/1‑143 (MQ=255)                                                                                                         gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCACCTTTCAGCTCGGTATCGAACCAACACTCttt                                >  2:654811/1‑143 (MQ=255)                                                                                                          tGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTg                               >  2:253675/1‑143 (MQ=255)                                                                                                           gATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCCGGAAGATATCCTAATCCTGTGCCGCTGTCAGTTCGGTATCGAACAGAGACCCTTTGa                              >  1:797325/1‑143 (MQ=255)                                                                                                            aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                             >  1:654742/1‑143 (MQ=255)                                                                                                             ccggcattgcggtgTGATTGAGGGGGAGGCTCTGCGTCGCCTCTTCTGCTTTGCGCGGGTTGCCGTACCGCAGCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                            <  1:65488/130‑1 (MQ=255)                                                                                                                gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGCGCCGCTTTCAGTTCGCGATGGAAACAACACCCTTTGGAACGc                         >  1:834923/1‑143 (MQ=255)                                                                                                                 cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCAGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                        <  1:478951/143‑1 (MQ=255)                                                                                                                 cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCAACACCATCCGCAGGAAGATGGCGTACTCCTGTGCCGCTTTCAGCGCGGTATCGAACAAACACTCTTTGAAACGcg                        >  1:373047/1‑142 (MQ=255)                                                                                                                   tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCAACCACCCTCCGCAGGAAGATGTCGTAACCCTGTGCCGCTTTCAGTTCGGCATCGTCCAAACCCCACTTGACACGCCAt                      >  1:242825/1‑143 (MQ=255)                                                                                                                     tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGc                    <  1:954138/143‑1 (MQ=255)                                                                                                                       gccgttgttccgatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACTGTTCGCCGTACTCCACCTCCATCCGCAGGAAGATGTCGTACTCCTGTGCCGCTTTCAGTTCGGTATCGATCAAACACTCTTTGAAACGCCATGccc                  <  1:93325/131‑1 (MQ=255)                                                                                                                         cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTCTCCACGGTTCGCCGTACTCCACCCCCATCCGCCGGCAGATGTCGTAATCCTGTGCCGCTTTCAGCTCGGTATCGCACCAACCCTCTTTGAAACGCCATGccggg                >  1:216357/1‑140 (MQ=255)                                                                                                                              ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa           <  2:359585/143‑1 (MQ=255)                                                                                                                                 atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCGCGACGATCCGCAGGAAGATGTCGTACTCCTGTGCCGCTGTCAGTGCGGTATGGGACAAACCCTGTTTGTAAGGCCCTGCACCGGtggaggc        >  2:203198/1‑137 (MQ=255)                                                                                                                                       aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGCAACGAACACACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGAtt  >  1:554091/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                            TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132650‑2132926

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