Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,064,259	T→G	V93G (GTG→GGG)	BW25113_RS10440 →	GTPase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,064,259	0	T	G	96.2%	86.2 / ‑3.9	26	V93G (GTG→GGG)	BW25113_RS10440	GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (11/14);  minor base A (0/1);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGC  >  NZ_CP009273/2064172‑2064343                                                                                        |                                                                                     tACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGc                                                                                    <  1:85796/90‑1 (MQ=255)                            gTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGAc                                                             >  1:320968/1‑86 (MQ=255)                            gTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGAc                                                             <  2:320968/86‑1 (MQ=255)                                   gtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGcc                                                  >  1:165936/1‑90 (MQ=255)                                    tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                                 <  2:469431/90‑1 (MQ=255)                                    tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                                 <  2:165936/90‑1 (MQ=255)                                    tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGAGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                                 <  2:322830/90‑1 (MQ=255)                                         ggATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCa                                            <  1:310465/90‑1 (MQ=255)                                         ggATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCa                                            <  2:478365/90‑1 (MQ=255)                                             gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                >  2:488329/1‑66 (MQ=255)                                             gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                                <  1:488329/66‑1 (MQ=255)                                             gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                        >  1:165960/1‑90 (MQ=255)                                               gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                                 >  2:573429/1‑79 (MQ=255)                                               gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                                 <  1:573429/79‑1 (MQ=255)                                                                     cAGTGCGGACATCAGCAGGGGTTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt                <  1:93337/90‑1 (MQ=255)                                                                     cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt                <  2:521213/90‑1 (MQ=255)                                                                        tGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACgg                                                      <  2:172844/49‑1 (MQ=255)                                                                        tGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACgg                                                      >  1:172844/1‑49 (MQ=255)                                                                        tGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtct             >  1:499586/1‑90 (MQ=255)                                                                         gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGcc                            >  2:360125/1‑74 (MQ=255)                                                                         gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCg                           >  1:550630/1‑75 (MQ=255)                                                                         gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCg                           <  2:550630/75‑1 (MQ=255)                                                                         gCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtctc            >  1:61032/1‑90 (MQ=255)                                                                         gCGGACATCAGCAGGGGCTGTTTGTGGGGACGCAGGCCGGCAAAACGGAGCCCTGCCATGAATGGGATATGGcc                            <  1:360125/74‑1 (MQ=255)                                                                                tCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcaca     <  1:280718/90‑1 (MQ=255)                                                                                   gcagGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcacaggc  >  2:314386/1‑90 (MQ=255)                                                                                        |                                                                                     TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGC  >  NZ_CP009273/2064172‑2064343

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