Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,319,078	C→A	D267Y (GAT→TAT)	rcsC ←	histidine kinase RcsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,319,078	0	C	A	96.8%	102.5 / ‑4.3	31	D267Y (GAT→TAT)	rcsC	histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (16/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATCGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTAT  >  NC_000913/2318839‑2319314                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             gtCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAAtt                                                                                                                                                                                                                                    <  2:6991/249‑1 (MQ=255)     gccgGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTGCCTGCATCCAGCGAGGATAGCCCTTAATTTTACt                                                                                                                                                                                                                               >  2:82492/1‑251 (MQ=255)                       cGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTGCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGa                                                                                                                                                                                                             <  2:381620/251‑1 (MQ=255)                          ccGTACAACCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACAGGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCGGCATCCAGCGAGGATCGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAAt                                                                                                                                                                                                          <  2:2072/250‑1 (MQ=255)                                     tGCGAGAGTAAACAATGAAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCAAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGAc                                                                                                                                                                                               >  1:92657/1‑251 (MQ=255)                                                            cGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAAcc                                                                                                                                                                        >  1:342151/1‑251 (MQ=255)                                                                aaGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAAccccc                                                                                                                                                                     >  1:94738/1‑250 (MQ=255)                                                                  gCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATc                                                                                                                                                                  <  2:202941/251‑1 (MQ=255)                                                                             aaaaTTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTAcc                                                                                                                                                       <  2:238413/251‑1 (MQ=255)                                                                                                  aTGCGAATGCGTTCCAGCCCCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCAt                                                                                                                                  <  2:370993/251‑1 (MQ=255)                                                                                                    gCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATcc                                                                                                                                <  2:72960/251‑1 (MQ=255)                                                                                                     cGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCg                                                                                                                               <  2:361198/251‑1 (MQ=255)                                                                                                            gTTCCAGCACCTTATCCACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGACCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGAAACCCCCATCGGCAACGTACACGATAAGAAAAAGTTCTCCATCAGGATGGTc                                                                                                                        >  2:364737/1‑251 (MQ=255)                                                                                                             ttCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCt                                                                                                                       >  1:166618/1‑251 (MQ=255)                                                                                                                       ttATCAACCGGCACAGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACAcc                                                                                                             >  1:138477/1‑251 (MQ=255)                                                                                                                            aaCCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCaaa                                                                                                         >  1:326276/1‑250 (MQ=255)                                                                                                                               cGGCCCCGAATACACGATGCTTAGCGATGAGTGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCAGTTTCATCAAGAATGGTAAGCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCACCCGGATGGTCTGCTCTACATCCAAAAgcg                                                                                                     <  2:29056/251‑1 (MQ=255)                                                                                                                                   cccGAATACACGAGGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCGTCCGTATAGCCAAACCAGGAGCGTGCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAATGTTCTCCATCCGGATGGTCTGCTCGACTCCCAAAAGCGCCTg                                                                                                 <  2:410432/250‑1 (MQ=255)                                                                                                                                                          tGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTg                                                                          >  1:170552/1‑251 (MQ=255)                                                                                                                                                             gggTGGCAGCTTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTACGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCACCTGCAAcccgatcggacgtgcgtcgtgt                                                                        >  2:365420/1‑230 (MQ=255)                                                                                                                                                             gggTGGCAGAGTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTTCAAcc                                                                                            <  1:365420/230‑1 (MQ=255)                                                                                                                                                                                aGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATcccc                                                    >  2:242326/1‑251 (MQ=255)                                                                                                                                                                                    ccAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACAcc                                                                                                             >  1:195459/1‑190 (MQ=255)                                                                                                                                                                                    ccAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACAcc                                                                                                             <  2:195459/190‑1 (MQ=255)                                                                                                                                                                                     cAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGc                                               >  1:5898/1‑251 (MQ=255)                                                                                                                                                                                                            aGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATc                        <  2:209046/251‑1 (MQ=255)                                                                                                                                                                                                              ccAAACCAGGAGCGTTCCTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATcc                       >  1:355656/1‑250 (MQ=255)                                                                                                                                                                                                                               cTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGAt                         <  2:299027/231‑1 (MQ=255)                                                                                                                                                                                                                               cTGCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGAt                         >  1:299027/1‑231 (MQ=255)                                                                                                                                                                                                                                 gCATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCt    <  1:171510/250‑1 (MQ=255)                                                                                                                                                                                                                                  cATCCAGCGAGGATAGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTat  >  1:352229/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             TTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCAGCGAGGATCGCCCTTAATTTTACTTTCTGGTCCGGTAAGCGAAATCAGGGTATGACCATTTTCATCAAGAATGGTAACCCCCATCGGCAACGTACCCGGTAAGAAAAAGTTCTCCATCCGGATGGTCTGCTCGACACCCAAAAGCGCCTGCAACCGGTTCGCCAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTAT  >  NC_000913/2318839‑2319314

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