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	.	95.0%	72.9 / ‑2.8	20	coding (165/459 nt)	ybbJ	NfeD family protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (9/10);  total (9/11)
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

GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAA  >  NZ_CP009273/510064‑510226                                                                                        |                                                                            gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                           >  1:1321611/1‑88 (MQ=255) gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                           >  1:235100/1‑88 (MQ=255)                        gcgACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCAccagcca                                                   >  1:114635/1‑90 (MQ=255)                            cAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGacc                                               >  2:263955/1‑90 (MQ=255)                               ccATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGaccacc                                            <  1:20742/90‑1 (MQ=255)                                   ttCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCCCCAGGc                                        <  1:232467/90‑1 (MQ=255)                                      caccacAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCgg                                     <  2:198321/90‑1 (MQ=255)                                             gCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAc                              >  1:763064/1‑90 (MQ=255)                                              ccAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                             <  2:114635/90‑1 (MQ=255)                                                aGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGaccacc                                            >  2:628568/1‑73 (MQ=255)                                                aGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGaccacc                                            <  1:628568/73‑1 (MQ=255)                                                     gcgAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCa                      >  1:1290008/1‑90 (MQ=255)                                                       gAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcg                    >  2:1117431/1‑90 (MQ=255)                                                        agcagcGTCAGGATGGCAAACATTCCCCCTTGCC‑CTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgc                   <  1:1312301/90‑1 (MQ=255)                                                         gcagcGTCAGGATGGCAACCATTCCCCCTTG‑CACCCCCAACCCAGCGGCCCCAGCCAGACCCCCAGGCCGGTAATCACCGCTGCCAcgcc                  <  1:924643/90‑1 (MQ=255)                                                                 aGGATGGCAAACATTCCCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCcaca          <  2:32617/90‑1 (MQ=255)                                                                 aGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCcaca          >  2:1246284/1‑90 (MQ=255)                                                                  ggATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgcc                  >  1:327629/1‑81 (MQ=255)                                                                  ggATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgcc                  <  2:327629/81‑1 (MQ=255)                                                                  ggATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCacaa         <  1:654045/90‑1 (MQ=255)                                                                      ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa     >  1:298553/1‑90 (MQ=255)                                                                         aaaCATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATaa  <  2:235100/90‑1 (MQ=255)                                                                                        |                                                                            GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAA  >  NZ_CP009273/510064‑510226

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