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	100.0%	75.2 / NA	25	G25G (GGA→GGG)	BW25113_RS26130	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (11/14);  total (11/14)

GGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGG  >  NZ_CP009273/1417701‑1417864                                                                                   |                                                                                   ggTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccg                                                                            >  2:1577981/1‑88 (MQ=18)   tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccgac                                                                          <  1:1169956/90‑5 (MQ=18)   tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑cccccccgac                                                                          <  2:625294/90‑5 (MQ=18)      tcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTgg                                                                        >  1:660371/1‑90 (MQ=255)                cacaTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                              <  2:1508164/90‑1 (MQ=255)                     cAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                         <  1:84580/90‑1 (MQ=255)                            gTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatta                                                  >  1:252331/1‑90 (MQ=255)                              cTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattatt                                                >  2:1459636/1‑90 (MQ=255)                              cTGTCAGTCAGCGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattatt                                                >  1:846189/1‑90 (MQ=11)                                gtcagtcaGTGCGTGAAGCCACCCCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTc                                              <  1:798797/90‑1 (MQ=11)                                 tcagtcagTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCa                                             >  1:432566/1‑90 (MQ=255)                                        gtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCg                                      <  2:903771/90‑1 (MQ=255)                                        gtGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCg                                      >  1:903771/1‑90 (MQ=255)                                              gAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCt                                                                                     <  2:832254/37‑1 (MQ=40)                                              gAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCt                                                                                     >  1:832254/1‑37 (MQ=40)                                                aGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGAc                              <  2:1412185/90‑1 (MQ=255)                                                 gCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCGCCCTCCGGGAACGACt                             <  1:291933/90‑1 (MQ=17)                                                  ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                            <  1:240815/90‑1 (MQ=255)                                                  ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                            >  1:123852/1‑90 (MQ=255)                                                  ccaccaCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                            >  2:240815/1‑90 (MQ=255)                                                       accGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                       <  1:430866/90‑1 (MQ=255)                                                       accGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAat                       <  1:1481719/90‑1 (MQ=255)                                                        ccGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAatg                      <  1:491698/90‑1 (MQ=255)                                                            cTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGc                  <  1:1459636/90‑1 (MQ=255)                                                              ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAg                >  1:445738/1‑90 (MQ=255)                                                               cGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACt                             >  2:1566017/1‑76 (MQ=255)                                                               cGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACt                             <  1:1566017/76‑1 (MQ=255)                                                                  cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGccc            <  2:432566/90‑1 (MQ=255)                                                                     ggATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAAc         >  1:1423222/1‑90 (MQ=255)                                                                            gCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAgg  <  1:107820/90‑1 (MQ=255)                                                                                   |                                                                                   GGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCT‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGG  >  NZ_CP009273/1417701‑1417864

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