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	95.5%	145.4 / 0.0	44	coding (278/711 nt)	BW25113_RS10440	GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (22/20);  minor base . (1/1);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAG  >  NZ_CP009273/2064172‑2064341                                                                                        |                                                                                   tACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGc                                                                                  >  2:299047/1‑90 (MQ=255)     gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                              >  1:18593/1‑90 (MQ=255)     gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                              >  1:800179/1‑90 (MQ=255)       ggCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTtg                                                                            >  1:272701/1‑90 (MQ=255)         cTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTtgtg                                                                          >  2:474039/1‑90 (MQ=255)           gATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCGGGGGCTGTTTgtggt                                                                        >  1:320255/1‑90 (MQ=255)            aTTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTgtggtg                                                                       >  2:899317/1‑90 (MQ=255)               aaaGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACg                                                                    >  1:174816/1‑90 (MQ=255)                     gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                              <  2:546035/90‑1 (MQ=255)                     gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                              >  2:25022/1‑90 (MQ=255)                     gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                              >  1:722003/1‑90 (MQ=255)                                 ctgtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAg                                                  <  1:792213/90‑1 (MQ=255)                                   gtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGcc                                                >  1:393642/1‑90 (MQ=255)                                    tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                               <  2:393642/90‑1 (MQ=255)                                    tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                               <  1:25022/90‑1 (MQ=255)                                       gtgGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGc                                            >  1:689263/1‑90 (MQ=255)                                       gtgGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGGGGTGGCGCGGGCCGACAAAACGGAGCCCTGc                                            >  2:770040/1‑90 (MQ=255)                                         ggATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCa                                          <  1:116313/90‑1 (MQ=255)                                          gATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCAt                                         <  1:900337/90‑1 (MQ=255)                                           aTGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATg                                        >  2:729506/1‑90 (MQ=255)                                             gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCTGG‑GCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     >  2:474701/1‑90 (MQ=255)                                             gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                      >  1:521331/1‑90 (MQ=255)                                             gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                      >  2:386807/1‑90 (MQ=255)                                             gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                      >  1:324769/1‑90 (MQ=255)                                              tgtcttttTGCCGACCCACCCTGCAGTGCGGACACCAGCAGGGGCTGTTTGTGGTCACGCAGGCCGCCAAAACGGAGCCCTGCCATGAAt                                     <  2:704917/86‑1 (MQ=255)                                              aGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     <  2:429363/90‑1 (MQ=255)                                              aGTATTTCTGGCGAACCATCAGGCAGTGAGGCCATCCGCAGGGGTTGTTTGTGGTGACCCAGGCCGACAAACCGGAGCCCTGCCATGAAt                                     <  2:113389/90‑1 (MQ=255)                                               gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATg                                    >  1:743353/1‑90 (MQ=255)                                                         cGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGcc                          <  2:722003/90‑1 (MQ=255)                                                            cacaTCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     <  2:687420/76‑1 (MQ=255)                                                            cacaTCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     >  1:687420/1‑76 (MQ=255)                                                                    gCAGTGCGGACATCAGCTGG‑GCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:474701/90‑1 (MQ=255)                                                                     cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  2:800179/90‑1 (MQ=255)                                                                     cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:885886/90‑1 (MQ=255)                                                                     cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  2:521331/90‑1 (MQ=255)                                                                     cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:907203/90‑1 (MQ=255)                                                                     cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:729506/90‑1 (MQ=255)                                                                        tGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGGCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtct           >  1:541037/1‑90 (MQ=255)                                                                            gACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTg       >  1:75877/1‑90 (MQ=255)                                                                               aTCAGCGGGGGCTGTTTGTGGGGCCGCAGGCCGACAAAACGGGGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcac    <  1:440292/90‑1 (MQ=255)                                                                                tCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcaca   <  2:781671/90‑1 (MQ=255)                                                                                 cagcagGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcacag  >  1:213163/1‑90 (MQ=255)                                                                                 cagcagGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcacag  <  1:116479/90‑1 (MQ=255)                                                                                     agGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  2:787254/74‑1 (MQ=255)                                                                                        |                                                                                   TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAG  >  NZ_CP009273/2064172‑2064341

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