Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,751,849	G→A	D19N (GAT→AAT)	recN →	DNA repair protein RecN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,751,849	0	G	A	93.8%	97.2 / ‑1.0	32	D19N (GAT→AAT)	recN	DNA repair protein RecN
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (13/17);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

CGTCGCTGTGATTACCATCTGAATCTGATTCATCCGAAAGATTACAGTTATTTCAACACATTAAGCACCAAGCTCGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGTG  >  NC_000913/2751611‑2752082                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           cgtcgCTGTGATTACCATCTGAATCTGATTCATCCGAAAGATTACAGTTATTTCAACACATTAAGCACCAAGCTCGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGc                                                                                                                                                                                                                                 <  1:239335/250‑1 (MQ=255) cgtcgCTGTGATTACCATCTGAATCTGATTCATCCGAAAGATTACAGTTATTTCAACACATTAAGCACCAAGCTCGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCg                                                                                                                                                                                                                                >  2:340556/1‑251 (MQ=255)                                      aaGATTACAGTTATTTCAACACATTAAGCACCAAGCTCGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTa                                                                                                                                                                                            >  2:91007/1‑250 (MQ=255)                                       aGATTACCGTTATTTCAACACATTAAGCCCCAAGCTAGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTaa                                                                                                                                                                                           >  2:359189/1‑250 (MQ=255)                                           tACAGTTATTTCAACACATTAAGCACCAAGCTCGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCt                                                                                                                                                                                       >  2:173170/1‑250 (MQ=255)                                                  aTTTCAACACATTAAGCACCAAGCTCTGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCTGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTACTGCAAt                                                                                                                                                                               <  2:82238/251‑1 (MQ=255)                                                            aTTAAGCACCAAGCTCGTCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGAACGGCGCGGGTAAATCTATTGCAATAGATGCCCTc                                                                                                                                                                     >  2:126373/1‑251 (MQ=255)                                                              tAAGCACCAAGCTCGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCgg                                                                                                                                                                   >  1:358470/1‑251 (MQ=255)                                                              tAAGCACCAAGCTCGGCTGGTCAAAAAAAATATTCTAAGTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTGTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACGGTGGTATCGTTCGTGAGCGTGAGATTAGTTTTCATAGCGGCATTACCGGAAGAACTGGCGAGGCCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCg                                                                                                                                                                    <  2:34988/250‑1 (MQ=255)                                                                            ggCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCggtg                                                                                                                                                     <  2:241639/251‑1 (MQ=255)                                                                                ggTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTc                                                                                                                                                  <  2:273383/250‑1 (MQ=255)                                                                                        aaaTTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAg                                                                                                                                         <  1:29022/251‑1 (MQ=255)                                                                                                                                 aaCCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGGTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACGGGCGAGACCGTCGCGGGTAAATCGATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGGCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTgc                                                                                               <  1:14526/251‑1 (MQ=255)                                                                                                                                       tttATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATAGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGt                                                                                         >  2:270747/1‑251 (MQ=255)                                                                                                                                                      aTAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCT‑AACCTTAACATAGTGCATACCGTCGCTTCTAGCGCTGACCTGAGCGCCCGTTTTtcgctagaaggt                                                                          >  2:43350/1‑241 (MQ=255)                                                                                                                                                        aaCAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATATATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCTGCGCTGCTGTCGGTGACCTGTGCGCCCGTTTTTCTCTGAAAGATa                                                                         <  1:359189/250‑1 (MQ=255)                                                                                                                                                        aaCAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATa                                                                         <  1:270747/250‑1 (MQ=255)                                                                                                                                                             tAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGcc                                                                     >  2:83967/1‑249 (MQ=255)                                                                                                                                                                 ggTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACGTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGAGGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCgg                                                                <  1:166933/250‑1 (MQ=255)                                                                                                                                                                     tttaatgcAGGAAAACGACTATGTTGGCACAAATGAACATCAGCAACTTTGCAATCGTTCGTGAGATTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAAAACGGCGAGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGAC‑GACACGGTGCGTACCGGCTCTGCTCGAGCTGCCCTGTGCGCCCGTTTTTATCTGAAAGATACGCCAGCGGctcgg                                                           >  1:295933/8‑249 (MQ=255)                                                                                                                                                                     tttCATTCAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTg                                                           >  1:99773/1‑251 (MQ=255)                                                                                                                                                                     tttCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTg                                                           >  1:295928/1‑251 (MQ=255)                                                                                                                                                                           aCAGGAAAACGAATATGTATGCACAACTAACCATCAGCAACGTTGCTATCGTTCGTGAGATTGAGATTATTTTTCATAGCGGCAAGACCGTAAAAACTGGCGAGACCGGCGCGGGTAAATCTATTTCAATAGATGCCCTCGGTCTTTGTCGCGGTGGTCGCGCTGAAGCC‑GACATGGCGCGTACCGGCGCCGCTCGCGCTGACCGGTGCGCCCGTTTTTCTCGGAAAGATACGCCAGCGGCTCTGCgctgg                                                     <  2:166968/251‑1 (MQ=255)                                                                                                                                                                               gAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTAACACT‑GAAAAGGATCGTCCCGGAAATGCGCGAGCTGCCCTGTGCGTACGTTTTTCACTGAAAGATACGCCAGCGGCTCTACGcgggctgg                                                 >  2:298335/1‑251 (MQ=255)                                                                                                                                                                                  aaCGACTAGGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGgaag                                              <  2:68480/251‑1 (MQ=255)                                                                                                                                                                                                   aaCTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGAc                             <  1:141584/251‑1 (MQ=255)                                                                                                                                                                                                       gACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACggg                          <  1:83967/250‑1 (MQ=255)                                                                                                                                                                                                          cATCAGCATCGTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCAGAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCAt                       <  2:198041/250‑1 (MQ=255)                                                                                                                                                                                                          cATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCAt                       >  1:198041/1‑250 (MQ=255)                                                                                                                                                                                                                  aCTTTGCTATCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTg              <  1:24045/251‑1 (MQ=255)                                                                                                                                                                                                                          aTCGTTCGTGAGCTTGAGATTAATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTtcgt        <  1:24688/249‑1 (MQ=255)                                                                                                                                                                                                                               tcgtGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGAGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGTg  <  1:91007/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           CGTCGCTGTGATTACCATCTGAATCTGATTCATCCGAAAGATTACAGTTATTTCAACACATTAAGCACCAAGCTCGGCTGGTCAAAAAAATTATTCTAATTTTACGCCAGCCTCTTTACTGTATATAAAACCAGTTTATACTGTACACAATAACAGTAATGGTTTTTCATACAGGAAAACGACTATGTTGGCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCC‑GACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACGCCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGTG  >  NC_000913/2751611‑2752082

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