Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	2,175,305	T→C	13.0%	L253L (TTA→TTG)	U712_10880 ←	putative glycosyltransferase ypjH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,175,305	0	T	C	13.0%	109.9 / 8.3	46	L253L (TTA→TTG)	U712_10880	putative glycosyltransferase ypjH
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (3/3);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGCACACGTTTGACCTTTC  >  CP006881/2175169‑2175438                                                                                                                                         |                                                                                                                                      aaGGCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTaa                                                                                                                                     <  4:255288/139‑1 (MQ=255)   ggCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATc                                                                                                                                   >  1:45334/1‑139 (MQ=255)       cTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCata                                                                                                                               >  1:97905/1‑139 (MQ=255)                ccATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGAc                                                                                                                      <  2:45334/139‑1 (MQ=255)                    cGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGc                                                                                                                  >  3:155089/1‑139 (MQ=255)                       cTCAAGAACGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTca                                                                                                               <  1:234198/139‑1 (MQ=255)                         cAAGCAGGACAAGGCCAAAGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTcg                                                                                                               >  4:141010/1‑136 (MQ=255)                         cAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTca                                                                                                               >  4:191414/1‑137 (MQ=255)                          aagaagGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTAGGTCTTCTAATCCATATTTTCTGACAAGCTCACAg                                                                                                            >  6:111662/1‑139 (MQ=255)                                  cAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGa                                                                                                    <  7:103762/139‑1 (MQ=255)                                  cAAGGCCGAAGCGTTCTTTTTCAGCCAGCAGAAGCTTTAAATCGCTAACGGAATAACGCTCCTCAACCAGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTTCTCCATATTTTCTGACAAGCTCACAGGCTGTCGa                                                                                                    <  8:124362/139‑1 (MQ=255)                                   aaGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGAt                                                                                                   <  6:243990/139‑1 (MQ=255)                                        cGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTtctc                                                                                              >  2:212714/1‑139 (MQ=255)                                          aaGCTTTCTTTTTCAGATAGCAGCAGCTTCAAATCGCTAATAGAATAAAGATCTTCAACACGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTtctccg                                                                                            >  8:153889/1‑137 (MQ=255)                                                    tttCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGtctc                                                                                   <  2:69176/138‑1 (MQ=255)                                                    tttCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGtctc                                                                                   <  2:64290/138‑1 (MQ=255)                                                            agcagcagcTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGTTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTaaaa                                                                          >  3:80649/1‑139 (MQ=255)                                                                              cTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGttt                                                        >  6:309026/1‑139 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGa                                                                                                    >  3:15337/1‑78 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGa                                                                                                    <  4:15337/78‑1 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTTCCGGCGATATTGCgg                                       >  6:145120/1‑139 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTTCg                                        >  3:292149/1‑138 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCg                                        <  6:46650/138‑1 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCg                                        <  4:155089/138‑1 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTGTCGGTTTTTCGGCGAGAtttggg                                       >  2:332025/1‑135 (MQ=255)                                                                                               ctcAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCGTTCGTTTTTCCGGCGCTATTGCgg                                       >  2:225926/1‑139 (MQ=255)                                                                                                 cAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGaa                                     <  5:54821/139‑1 (MQ=255)                                                                                                  aaCACGTTCTTGATTTCCAAGAATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGTACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGaa                                     <  5:261879/138‑1 (MQ=255)                                                                                                      cGGTCTTGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGTTTTGCCAGCGATATTGCGGAacacac                                <  3:141010/139‑1 (MQ=255)                                                                                                      cGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGGCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAacacac                                <  5:201431/139‑1 (MQ=255)                                                                                                         tCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGa                             >  7:67002/1‑139 (MQ=255)                                                                                                           ttGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTcgca                                                                                                             >  6:62387/1‑57 (MQ=255)                                                                                                           ttGATTTCCAAGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTcgca                                                                                                             <  5:62387/57‑1 (MQ=255)                                                                                                               tttCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACg                        <  3:191414/138‑1 (MQ=255)                                                                                                                   cAAGCATTAAAACTTTGTCTTCTAATCCATATTTTCTGACAACCTCACAGGCTGTCGCGTTCACTTGACCGTCTCCCACTAAAAGCAGCTTCGCTTTCTGTTTGACCGCGATATTACGCTGCACACGGATTAAGTCCtc                   >  2:91938/1‑138 (MQ=255)                                                                                                                   cAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGTACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTTCCGGCGATATTGCCGAAAACACGGAAGACGTCCTg                   >  7:12028/1‑139 (MQ=255)                                                                                                                   cAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTg                   >  5:273937/1‑139 (MQ=255)                                                                                                                   cAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCCTCTCAAACAAAAATCAGCTTCGCTTTCTTTTTTCCGCCGATAATACGGAACACACGGAAGACGTCCtc                   >  8:23673/1‑138 (MQ=255)                                                                                                                    aaGCATTAAGACTTGGTCTTCCAAGCCATATTTTCTGATAAGCTCGCAGGCTGTCGATTTCTCCGGACCGTCTCCGACTAAAAGCAGCTTCGCTTTCGTTTTGCCAGCGATATTGCGGAACACACGGATGACATCCTGc                  <  7:153889/139‑1 (MQ=255)                                                                                                                     aGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGca                 <  8:67002/139‑1 (MQ=255)                                                                                                                     aGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGca                 <  2:97905/139‑1 (MQ=255)                                                                                                                         ttAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTGGTTTTGGCGGCGATATTGCGGCACACACCGAAGTCGTCCTGcccccg             >  3:244500/1‑134 (MQ=255)                                                                                                                               cTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGCACACGTTTGAc       <  5:309026/139‑1 (MQ=255)                                                                                                                               cTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGCACACGTTTGAc       <  5:145120/139‑1 (MQ=255)                                                                                                                               cTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGCACACGTTTGAc       <  5:111662/139‑1 (MQ=255)                                                                                                                               cTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGCACACGTTTGAc       <  1:332025/139‑1 (MQ=255)                                                                                                                                    tcttctAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGCACACGTTTGACCTTTc  >  4:322144/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                      AAGGCACTCCACAAGCCATCGCCTCAAGAAGGACAAGGCCGAAGCTTTCTTTTTCAGACAGCAGCAGCTTTAAATCGCTAATGGAATAAAGCTCCTCAACACGGTCTTGATTTCCAAGCATTAAAACTTGGTCTTCTAATCCATATTTTCTGACAAGCTCACAGGCTGTCGATTTCTCTGGACCGTCTCCAACTAAAAGCAGCTTCGCTTTCGTTTTGCCGGCGATATTGCGGAACACACGGATGACGTCCTGCACACGTTTGACCTTTC  >  CP006881/2175169‑2175438

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑