Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,519	C→T	T48T (ACG→ACA)	dcuS ←	sensor histidine kinase DcuS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351,519	0	C	T	96.5%	87.9 / ‑3.4	29	T48T (ACG→ACA)	dcuS	sensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (11/17);  total (11/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.42e-01

TAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATA‑GGGGAAGAGACTCTGGCA  >  NC_000913/4351284‑4351764                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       tAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGAt                                                                                                                                                                                                                                          <  2:168297/250‑1 (MQ=255)        gTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCgag                                                                                                                                                                                                                                   <  1:6868/250‑1 (MQ=255)                    aaaTGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATc                                                                                                                                                                                                                          >  2:230314/1‑247 (MQ=255)                     aaTGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATg                                                                                                                                                                                                                     <  1:99511/251‑1 (MQ=255)                                     gCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCAcc                                                                                                                                                                                                     >  2:111296/1‑251 (MQ=255)                                                     cGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGCACAGTACCGCAc                                                                                                                                                                                     >  2:90337/1‑251 (MQ=255)                                                      gAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACt                                                                                                                                                                                    <  1:230314/251‑1 (MQ=255)                                                           gCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGAGCTTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCTAGAAGTAAATCAGATGCACCGCCAATAGCACCGAGAACAGTACCGCACTGACCa                                                                                                                                                                               <  1:144534/251‑1 (MQ=255)                                                              acgACTTTGCATATCGGAAACGACAATAAACAGCAGATCGTTGAGTTTGCGTCCGGCGTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTa                                                                                                                                                                            <  1:222300/249‑1 (MQ=255)                                                                      gCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACt                                                                                                                                                                    >  2:13516/1‑251 (MQ=255)                                                                        atatCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgt                                                                                                                                                                  >  1:209811/1‑251 (MQ=255)                                                                                          aaCAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGAGCATTAAGATCACTGTGGTACTCAATTTCATCgg                                                                                                                                                <  1:180054/251‑1 (MQ=255)                                                                                           aCAGCCGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCgg                                                                                                                                                >  1:218686/1‑250 (MQ=255)                                                                                                agATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGttt                                                                                                                                          >  2:287842/1‑251 (MQ=255)                                                                                                       tgcgtttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTaa                                                                                                                                    <  1:1237/250‑1 (MQ=255)                                                                                                                            ggATGGCCTGAATGCCACTCTGCTGCGGGTTTTTCTGGAAGCCCTGACGGATTTCTGGCGAGTCGGCGGGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCGCATCTCGTGTCATAGCACTGATTTGCGAGAAGTAAATCAGATGCACCAGCAAGAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTATTCAATGTCGGCGAACGGTTACGTATCATGCGGTGGGGCAATGAATGt                                                                                                              <  2:276192/250‑1 (MQ=255)                                                                                                                             gatggccggtgggCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtc                                                                                                             <  2:6659/239‑1 (MQ=255)                                                                                                                              aTGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGCGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtct                                                                                                            >  1:91367/1‑251 (MQ=255)                                                                                                                                     ggATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGc                                                                                                     <  1:172026/251‑1 (MQ=255)                                                                                                                                      gATGCCACTCTCCGGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCt                                                                                                    <  1:74801/251‑1 (MQ=255)                                                                                                                                                 ccTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGa                                                                                         <  2:100288/251‑1 (MQ=255)                                                                                                                                                                                     gCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTGGTGTAACATGCGGTAGTGCAATGAATGTCTCATCAGCTTCCTTGTGTTACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTa                                                     <  1:153208/251‑1 (MQ=255)                                                                                                                                                                                              cGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAATCATTATCTCTGATGAGGCGGGTAATTCaaatg                                            <  1:157711/251‑3 (MQ=255)                                                                                                                                                                                                      cgcgCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCAGATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGa                                    <  2:179686/251‑1 (MQ=255)                                                                                                                                                                                                                       cTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGAGGAAGAAt                                  >  1:108725/1‑236 (MQ=255)                                                                                                                                                                                                                       cTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAAt                                  <  2:108725/236‑1 (MQ=255)                                                                                                                                                                                                                                 aGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCTAGAACAGTACCTCTCTGATCATTAAGATCACTGTTGTACTCAATTTCATCGGATTTTTTCGTAACATGCGTGAGGGCAATGAATGTCTCATCATCTTCCGTGTTTGACAAATTTCTTAAGCAGTATTTCTGTTGATGCTGGTTATTGAAAGTTAGTAAGAATGATTGGCTATATA‑GGGGAAGAGAct        <  1:23061/251‑1 (MQ=255)                                                                                                                                                                                                                                 aGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGGGGAAGAATGATTGGGTATATATGGGGagggacc         >  2:17684/1‑245 (MQ=255)                                                                                                                                                                                                                                        ctcGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGAGTAAGAATGATTTGTCATATA‑TGGGAagaaaatcctgga  >  2:190674/1‑240 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       TAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATA‑GGGGAAGAGACTCTGGCA  >  NC_000913/4351284‑4351764

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