Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,635,706	C→A	T62T (ACC→ACA)	creA →	PF05981 family protein CreA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,635,706	0	C	A	94.3%	100.7 / ‑3.6	35	T62T (ACC→ACA)	creA	PF05981 family protein CreA
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (16/17);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCG  >  NC_000913/4635461‑4635949                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    cAATATGTTATTTACAGTGACGAGCTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACAGTATTTAAAATGCTCGGACAGGATCACAAAATT‑GTTGTGGAAACCTTTGCCGATCCCGATGGGAAAAATGTCAACTGTTATGTACGCGGGGCGAAAAACggtgg                                                                                                                                                                                                                                                >  2:186682/1‑251 (MQ=255)     atGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTAtt                                                                                                                                                                                                                                            >  2:317383/1‑251 (MQ=255)     atGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTAt                                                                                                                                                                                                                                             <  1:428036/250‑1 (MQ=255)     atGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAAGATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGTTGGTAt                                                                                                                                                                                                                                             <  2:303382/250‑1 (MQ=255)                 gTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGGTGTGAAAAATGTCACCTGTTATGTGGGCCGGGCGAAAACAGGTGGTATTAAAGGGGGAtcg                                                                                                                                                                                                                                >  2:221698/1‑249 (MQ=255)                      gAACTAATTTGTCTTGTAACAGATAAAAATGGTAACAATATAAAATGCAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGTCCAGTGGCTCATTCAGAAGAGATTGGTTCGGTCGACTCCGGATTTAAAATGATCGGCCCTGATCACAAAATT‑GTTGTGGAGGCCTTTGTCGATCCTGATTTGAAAAATGTCGCCTGGTATGTCTGCCGTGCGAAAATAGTTGGTATGAAAGAGTGATTGGGTCt                                                                                                                                                                                                                           <  2:473917/251‑1 (MQ=255)                       aaCTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTg                                                                                                                                                                                                                          <  1:114570/251‑1 (MQ=255)                       aaCTAATTGCTCGTGGAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTg                                                                                                                                                                                                                          <  2:87778/251‑1 (MQ=255)                                       aaTAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCg                                                                                                                                                                                                          >  1:15686/1‑251 (MQ=255)                                       aaTAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTAGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCAGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCg                                                                                                                                                                                                          >  1:466636/1‑251 (MQ=255)                                                  aTGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGCTGCGGCTAtt                                                                                                                                                                                               >  1:321813/1‑251 (MQ=255)                                                     gTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGACCTTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGACGATACCTACGATGCGGctttttgt                                                                                                                                                                                            >  2:64416/1‑244 (MQ=255)                                                         cAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTCAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                                                                        >  1:443460/1‑251 (MQ=255)                                                               gAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCGCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGTTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAg                                                                                                                                                                                  <  2:141413/251‑1 (MQ=255)                                                                          aTTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                                                                        <  2:238206/234‑1 (MQ=255)                                                                          aTTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                                                                        >  1:238206/1‑234 (MQ=255)                                                                                                                     tCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAgg                                                                                                                            >  1:99261/1‑251 (MQ=255)                                                                                                                            gaagaGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCg                                                                                                                                                                                                                       <  2:157659/153‑1 (MQ=255)                                                                                                                            gaagaGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCg                                                                                                                                                                                                                       >  1:157659/1‑153 (MQ=255)                                                                                                                                                     tATTTAAAATGATCGGCCCGGATCACAAATTT‑GTTGTGGCAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCACCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCAGATTGAACTGGAGGAACGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAATATTCAAATAACGCACGGCCCTGGTCttt                                                                                            >  1:478109/1‑251 (MQ=255)                                                                                                                                                      aTTTAAAATGATCGGCCCGGCTCACAAAATT‑GTTGTGGAAGCCTTTAACGCTCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAGCAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGCTACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGGCGGATCGTATTATAAACGGCAAAGCACAGGGCGCGGAAGTATTCAAAAACCGCACGACCCTGGCCttt                                                                                            >  2:92319/1‑250 (MQ=255)                                                                                                                                                      aTTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGa                                                                                                                                                                                                         >  2:91146/1‑141 (MQ=255)                                                                                                                                                      aTTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGa                                                                                                                                                                                                         <  1:91146/141‑1 (MQ=255)                                                                                                                                                           aaaTGATCGGACCGGAACACAAAATTTGTTG‑GGAAGCATTTGACGATCCCGATGTGAAAAATGGCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCg                                                                                      <  2:310169/251‑1 (MQ=255)                                                                                                                                                                            cacaAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCtt                                                                     <  1:402848/251‑1 (MQ=255)                                                                                                                                                                                      gttgtGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTatgacg                                                             <  1:221698/250‑3 (MQ=255)                                                                                                                                                                                              aGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAaacgca                                                    <  1:430948/251‑1 (MQ=255)                                                                                                                                                                                                ccTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTGCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAaacgcaac                                                  <  2:466636/251‑1 (MQ=255)                                                                                                                                                                                                  tttGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTCCAGGTCGTGCGCTTTTAAGACGCAAaacgcaccgc                                                >  2:466181/1‑251 (MQ=255)                                                                                                                                                                                                     gatgatCCCCATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGcgct                                             >  1:258327/1‑251 (MQ=255)                                                                                                                                                                                                               aTGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTACGATGCGGCTATTTCTTGTCAGCACGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCCAAGCTCATGGCGAGGTAGTATTCAAAAAACGAACGTACCTGGTCTTTAAGTCGTTACAGGTAGTGCGCTTTTATGACGCCAAACGAAACGCGCTCGCTTATCTg                                   >  1:203630/1‑251 (MQ=255)                                                                                                                                                                                                               aTGTGAAAAAGGTCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGTGGATTGGGTCTGGCGGAATATACCTCCGAGGCGGCTATTTCTGGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTAGTCAAAAAACGCACGTCCCTGGTCGTTAAGTCGTTACAGGTCGGGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTg                                   <  2:468747/251‑1 (MQ=255)                                                                                                                                                                                                                                 gTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAg                   <  1:317383/249‑1 (MQ=255)                                                                                                                                                                                                                                         aGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGCTTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTACGACAAAGTTGTAGAAgg         >  1:467829/1‑251 (MQ=255)                                                                                                                                                                                                                                                cGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGTCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGGCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCg  <  2:439403/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCG  >  NC_000913/4635461‑4635949

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