Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,417,788	(C)6→7	coding (69/120 nt)	BW25113_RS26130 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,417,782	1	.	C	94.4%	46.9 / ‑2.6	18	coding (75/120 nt)	BW25113_RS26130	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (8/9);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGG  >  NZ_CP009273/1417721‑1417864                                                               |                                                                                   cAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                         <  1:340625/90‑1 (MQ=255)  aaaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑CCCCCCGACTGGCAGACACCGCTGAACgg                                                       >  2:82555/1‑90 (MQ=35)        gTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatta                                                  >  2:72426/1‑90 (MQ=255)          cTGTCAGTCAGTGCGTGAATCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                         >  2:263201/1‑81 (MQ=11)          cTGTCAGTCAGTGCGTGAATCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                         <  1:263201/81‑1 (MQ=11)          cTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGac                                                                   >  1:457551/1‑71 (MQ=255)          cTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGac                                                                   <  2:457551/71‑1 (MQ=255)                  cagtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                        >  2:368160/1‑90 (MQ=255)                      gCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCggg                                    >  2:152487/1‑90 (MQ=255)                        gTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGaa                                  <  1:564050/90‑1 (MQ=14)                              ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                            >  2:427718/1‑90 (MQ=255)                               caccacCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTgg                           >  2:541882/1‑90 (MQ=255)                                         tCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCa                 >  1:424103/1‑90 (MQ=255)                                            ggCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGc              <  1:351254/90‑1 (MQ=255)                                             gCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGcc             <  1:72426/90‑1 (MQ=255)                                              cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc            <  1:541882/90‑1 (MQ=255)                                                   aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACtt       <  2:424103/90‑1 (MQ=255)                                                        gCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAgg  <  2:535658/90‑1 (MQ=255)                                                               |                                                                                   CAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGG  >  NZ_CP009273/1417721‑1417864

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