Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,958,917	(T)8→7	coding (74/1266 nt)	rhlB ←	ATP‑dependent RNA helicase RhlB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,958,910	0	T	.	100.0%	87.2 / NA	21	coding (81/1266 nt)	rhlB	ATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (8/13);  total (8/13)

ACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCAT  >  NZ_CP009273/3958823‑3958990                                                                                        |                                                                                 aCCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCttt                                                                                <  1:256198/90‑1 (MQ=255)            ccTGCCCGGCTCCGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGc                                                                    <  1:287837/90‑1 (MQ=255)            ccTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGc                                                                    <  2:34045/90‑1 (MQ=255)              tGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCtt                                                                  >  2:161452/1‑90 (MQ=255)                  cGGTTACGTCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTAc                                                              <  2:436592/90‑1 (MQ=255)                     cTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAAc                                                           <  2:118809/90‑1 (MQ=255)                                  cGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                              >  1:502810/1‑90 (MQ=255)                                       gCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                                         <  2:260167/90‑1 (MQ=255)                                                    ggCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttc                            <  2:103815/90‑1 (MQ=255)                                                         gTGCCTGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttc                       >  1:457235/1‑90 (MQ=255)                                                              tGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa                  >  1:184210/1‑90 (MQ=255)                                                              tGAATGGGCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa                  >  1:78372/1‑90 (MQ=255)                                                                   gggCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtgt             >  1:96851/1‑90 (MQ=255)                                                                     gCGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATgtgt           >  2:288336/1‑90 (MQ=255)                                                                      cGTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtt          <  2:502810/90‑1 (MQ=255)                                                                       gTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGttt         >  1:181747/1‑90 (MQ=255)                                                                       gTACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGttt         <  2:7955/90‑1 (MQ=255)                                                                        tACAGTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGtttt        <  2:464761/90‑1 (MQ=255)                                                                            gTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTc    <  1:31602/90‑1 (MQ=255)                                                                            gTTATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTc    <  2:446808/90‑1 (MQ=255)                                                                              tATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCAt  <  2:240463/90‑1 (MQ=255)                                                                              tATGAAACCC‑TTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCAt  <  1:211050/90‑1 (MQ=255)                                                                                        |                                                                                 ACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATGTGTTTTGCTCAT  >  NZ_CP009273/3958823‑3958990

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