Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	510,152	Δ1 bp	coding (164/459 nt)	ybbJ ←	NfeD family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	510,151	0	C	.	87.5%	20.8 / ‑2.6	8	coding (165/459 nt)	ybbJ	NfeD family protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (3/4);  total (3/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAT  >  NZ_CP009273/510068‑510224                                                                                    |                                                                          ttttGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTccc                                                                    <  2:68074/90‑1 (MQ=255)                     cgACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCAccagccag                                                <  2:209433/90‑1 (MQ=255)                           ccATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCC‑CTCCCAGCCCAGCGGCACCAGCCAGaccccc                                          <  1:7628/90‑1 (MQ=255)                             aTTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAg                                        <  1:312177/90‑1 (MQ=255)                                                                  ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa   >  1:201590/1‑90 (MQ=255)                                                                   gCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAt  >  2:286194/1‑90 (MQ=255)                                                                           tACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                           >  1:58522/1‑57 (MQ=255)                                                                           tACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                           <  2:58522/57‑1 (MQ=255)                                                                                    |                                                                          TTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAT  >  NZ_CP009273/510068‑510224

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