Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	280,792	C→A	16.0%	P83P (CCG→CCT)	fliY ←	cystine‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	280,792	0	C	A	16.0%	52.1 / 4.7	25	P83P (CCG→CCT)	fliY	cystine‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (9/12);  new base A (2/2);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

CAGTGCCTGGATGCCGGAGACAGTGTAGGGCTCGGAGAAGTCGTACTTTTTCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAGGCCGACGTTGATGACGCC  >  NC_002947/280654‑280929                                                                                                                                           |                                                                                                                                          cAGTGCCTGGATGCCGGGGACAGTGTAGGGCTCGGAGAAGTCGTAATCCGTCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCtt                                                                                                                                <  1:86550/150‑1 (MQ=255)            tGCCGGAGACAGTGTAGGGCTCGGAGAAGTCGTACTTTTTCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGAccc                                                                                                                     >  1:110371/1‑150 (MQ=255)            tGCCGGAGACAGTGTAGGGCTCGGAGAAGTCGTACGTTTTCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGTCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGAccc                                                                                                                     <  2:319777/150‑1 (MQ=255)               cGGAGACAGTGTAGGGCTCGGAGAAGTCGTACTTTTTCTTGCGCTCTTCGGAGATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCa                                                                                                                  <  1:110403/150‑1 (MQ=255)                      aGTGTCGGGCTCGCAGAAGTCGTACTTTTCCTTGCGCTCTTCCGAAATGGTCCCCTGATTGACCCCTCCATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCCCTTGGTCCGCCGGATCCTGGCCTCGCCCCACAGCCCCtt                                                                                                           <  2:231886/150‑1 (MQ=255)                                      aaGTCGTACTTTTTCTTGCGCTCTTCGGAACTGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGAc                                                                                           <  1:202161/150‑1 (MQ=255)                                             aCTTTTTCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGcc                                                                                                                            <  2:558276/110‑1 (MQ=255)                                             aCTTTTTCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGcc                                                                                                                            >  1:558276/1‑110 (MQ=255)                                               tttttCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACt                                                                                  >  2:493070/1‑150 (MQ=255)                                                         gcTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCgg                                                                        >  2:777014/1‑150 (MQ=255)                                                          ctctTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGt                                                                       <  2:355815/150‑1 (MQ=255)                                                          ctctTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGt                                                                       <  1:493070/150‑1 (MQ=255)                                                                            ctaaTTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCGCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTc                                                     >  2:150253/4‑150 (MQ=255)                                                                                  gACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGaa                                               >  1:593996/1‑150 (MQ=255)                                                                                    ccACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTAGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGc                                             >  2:745954/1‑150 (MQ=255)                                                                                    ccACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTAGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGc                                             >  2:746098/1‑150 (MQ=255)                                                                                            gccAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCTACGTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGc                                     <  2:110371/149‑1 (MQ=255)                                                                                               aGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCggg                                  >  2:642601/1‑150 (MQ=255)                                                                                                            ccAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTAACACCCATCTCCTTTGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGa                     <  1:777014/150‑1 (MQ=255)                                                                                                             cAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAg                    <  2:594053/150‑1 (MQ=255)                                                                                                              aGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAgg                   >  1:437346/1‑150 (MQ=255)                                                                                                                   ggcgAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAGGCCGAc              <  1:642601/150‑1 (MQ=255)                                                                                                                         gATGCCGTCCCACTTGGTAGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAGGCCGACGTtgat        <  1:746098/150‑1 (MQ=255)                                                                                                                         gATGCCGTCCCACTTGGTAGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAGGCCGACGTtgat        <  1:745954/150‑1 (MQ=255)                                                                                                                               gTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAGGCCGACGTTGATGACGcc  >  2:544511/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                          CAGTGCCTGGATGCCGGAGACAGTGTAGGGCTCGGAGAAGTCGTACTTTTTCTTGCGCTCTTCGGAAATGGTCACCTGATTGACCACTACATCCAGGCGTTTGGATTCCAGCGCGGCGAGGATGCCGTCCCACTTGGTCGGCTGGATCTTGGCCTTGACCCCCAGCTCCTTGGCCAGCAGCTCGGACAGCTCCACTTCGAAGCCGGTCAGTTTGCCGTTTTCGTCCTGGAAGCTGAATGGCGGGTAGGTGCCCTCGAGGCCGACGTTGATGACGCC  >  NC_002947/280654‑280929

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑