Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,044,626	A→C	50.4%	E90A (GAG→GCG)	rcsC_5 →	Sensor histidine kinase RcsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,044,626	0	A	C	50.4%	‑2.3 / 9.3	29	E90A (GAG→GCG)	rcsC_5	Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base A (11/3);  new base C (13/2);  total (24/5)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CCCATACGCCGCGTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATG  >  ACB122/3044495‑3044754                                                                                                                                    |                                                                                                                                  cccATACGCCGCGTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACg                                                                                                                               >  7:330963/1‑136 (MQ=255)             gTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGCCGGCCcccccgct                                                                                                                    >  4:258674/1‑134 (MQ=255)             gTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACCGCCccccaaca                                                                                                                    >  2:333272/1‑135 (MQ=255)               tCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGTCGGCCcacccgcccc                                                                                                                  >  4:309146/1‑135 (MQ=255)               tCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGCCGGCCcccccgcccc                                                                                                                  >  3:6420/1‑135 (MQ=255)               tCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGACGGCCccccagcacc                                                                                                                  >  3:110401/1‑135 (MQ=255)                              tgCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTG‑ATTa                                                                                                  >  3:377529/1‑135 (MQ=255)                                 tGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGACGGCCCACCAGCACCCGATTCGCCTG‑ATTACCg                                                                                               >  2:95436/1‑135 (MQ=255)                                   gCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGACGGCCCCCCCGCACCAGATTCGCCTG‑ATTACCGAc                                                                                             >  7:115830/1‑135 (MQ=255)                                   gCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCCCCAGCACCAGCTTCGCCGG‑ATTACCgc                                                                                              >  3:108118/1‑133 (MQ=255)                                           ggTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCACCCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCCCGCCGGCCCACCACCACCAGAGTCGCCGG‑AGTACCGACACCgggcc                                                                                     >  5:355003/1‑131 (MQ=255)                                                gcgAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCACGGCCGTCACCCGCCTGCTCAGCGCGGCCCGGCCGCGGGCCGAGACGGCCCCCCAGCACCCGCGTGGCCTG‑CTTACCGGAGACGGCCCGGCACt                                                                               >  7:277698/1‑136 (MQ=255)                                                                     aCCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGACGGCCCCCCCGCACCAGATTCGCCAG‑AATAACGACAACGGCCCCGCCCTGCACGCGTACCAGAACGcc                                                            >  7:251691/1‑134 (MQ=255)                                                                       cAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGCCGGCCCCCCCGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCACTGCTCGCTTACCTGAACGCCgct                                                         >  4:192979/1‑133 (MQ=255)                                                                       cAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCCAGCCGGGCCCCCCCCCCCAGACTCCCCAG‑CTTCGCGACCACGGCCCGGCACGGCTCGCTGTCCTGAACGccc                                                           >  3:247820/1‑132 (MQ=255)                                                                         ggCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGACGGCCCCCCCGCACCAGAGTCGCCGG‑ATTACCGACAACCTCCCGGCCCTGATCGCTTTACGGGACGCCGaggag                                                      >  8:190928/1‑132 (MQ=255)                                                                            cATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGCCGGCCCCCCAGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCCCTGCTCGCTTACCGGAACGCCGCTCGGGTc                                                   >  8:360341/1‑136 (MQ=255)                                                                            cATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCCAGCCGGCCCACCCCCCCCAGGTTAGGCGG‑CTTTCCGCCAACCGCCCCGCACGGAGAGCGTACCTGAACGCCGAGCGGGt                                                    >  3:235210/1‑135 (MQ=255)                                                                             aTGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGCCGGCCCCCCAGCACCAGATTCGCCGGCATTACCCCCACCG‑CCCGGCCCTGCTCGCTTACCTGGACGCCGCTcgtgtc                                                   >  1:70841/1‑130 (MQ=255)                                                                                             ccGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGCCGGCCCCCCAGCACCAGATTCTCCGG‑CTTAACGACAACGGCCCGGCCCGGCTCGCGTACCGGAACGCCGCTCTGGTCTACGCATTCACcaac                                    >  4:40208/1‑134 (MQ=255)                                                                                                   cccgccgGCGCCGCGCGGCGGGGCCGGGGGCCGCCACGGCCCACCAGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGttt                             <  4:6420/134‑1 (MQ=255)                                                                                                   cccgccgGCGCCGCGCGCCCGGGCCGCGGGCCGCCACGGCCCACCGGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGttt                             <  4:194503/134‑1 (MQ=255)                                                                                                   aCCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGCGCCGGCCCCCCAGCCCCAGATTCGCCGG‑GTTACCGACCACCGCCCGGCACTGCTCGCGTACCTGAACGCCGCTCTGGGGTACGAATTCACCAACAAGGttt                             >  4:287988/1‑135 (MQ=255)                                                                                                   aCCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGCCGGCCCCCCAGCACCCGATTAGCCGG‑ATTTCCGCCCACCTCCCGGCCCTGATCGCTTACCTGAACGCCGATCTGGTCGACGAATTAACCAACAAgggt                              >  4:176201/1‑132 (MQ=255)                                                                                                   aCCGCCTGCTCAGCGAGGCCCGGCAGCGAGCCGAGACGGCACACCAGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCACTGATCGCTTACCTGACCGCCGATCTGGTCTACGAATTCACCACCAAGGttt                             >  6:11490/1‑135 (MQ=255)                                                                                                          gCTCAGCGAGGCCCGGCAGCGGGCCGCGACGGCCCCCCAGCCCCAGATTAGCCGG‑ATTAACGCCAACGGCCCGGCCCTGATCGCGTACCTGAACGCCGCTCGGGTCTACGCAATCACCAACAAGGTTTACGAaga                      >  1:291441/1‑135 (MQ=255)                                                                                                                        ggCTGCGGCCCGACACGGCCCACCAGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGcc       <  8:330963/136‑1 (MQ=255)                                                                                                                              gagCCGAGACGGCCCACCAGCACCAGATTCGCCTC‑ATTCCCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACCCGGTTTACGAAGAGTGGTATTGCTGGCCGCATg  <  5:11490/133‑1 (MQ=255)                                                                                                                              cggCCCACACGGCCGACCATCACCAGATTCGCCTG‑ATTACCGACACCGTCCCGGCACTGATCGCTTACCTGAACTCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTTTTATTGCTGGCCGCATg  <  6:231864/134‑1 (MQ=255)                                                                                                                                    |                                                                                                                                  CCCATACGCCGCGTTCCAGCATTCAGTGGTGCTGGCCGAGCAGGTACGCGAACGCACCGATGCGTTGAACCAGGCCATGGCCGAGCTCAAGGCCGTCAACCGCCTGCTCAGCGAGGCCCGGCAGCGGGCCGAGACGGCCCACCAGCACCAGATTCGCCTG‑ATTACCGACAACGTCCCGGCACTGATCGCTTACCTGAACGCCGATCTGGTCTACGAATTCACCAACAAGGTTTACGAAGAGTGGTATTGCTGGCCGCATG  >  ACB122/3044495‑3044754

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