Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	176,699	Δ1 bp	coding (976/1518 nt)	dgt →	dGTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	176,698	0	C	.	80.0%	20.1 / 4.7	10	coding (975/1518 nt)	dgt	dGTPase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base . (5/3);  total (7/3)
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.19e-01

CACGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGC  >  NZ_CP009273/176609‑176766                                                                                          |                                                                     cACGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGc                                                                      <  1:511201/90‑1 (MQ=255)   cGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGcgc                                                                    <  1:43916/90‑3 (MQ=255)           ggggCCAGCATGAGAAAGGCTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAG‑CGCAGTACGGa                                                           >  2:564751/1‑90 (MQ=255)                 aGCATGAGAAAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATc                                                      >  2:270729/1‑90 (MQ=255)                      gagaAAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAG‑CGCAGTACGGAAGATCAGtttt                                                >  1:132779/1‑90 (MQ=255)                               tCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAG‑CGCAGTACGGAAGATCAGTTTTTTATGTAtt                                       >  2:198714/1‑90 (MQ=255)                                            tggtggTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAG‑CGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAacac                          >  1:373027/1‑90 (MQ=255)                                              gtggtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAG‑CGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACAccc                        <  1:564751/90‑1 (MQ=255)                                                ggtTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCt                       >  1:27847/1‑90 (MQ=255)                                                          gCCTGGGAAAAATCACGCTCAAATAGTTTAAG‑CGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTgg            <  2:132779/90‑1 (MQ=255)                                                              gggAAAAATCACGCTCAAATAGTTTAAG‑CGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTAcc        >  2:233995/1‑90 (MQ=255)                                                                    aaTCACGCTCAAATAGTTTAAG‑CGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTAcgc  <  2:284362/90‑1 (MQ=255)                                                                                          |                                                                     CACGAAGCGTGGGGCCAGCATGAGAAAGGTTCGCTCTTTTCGCTGGTGGTTGAAAATGCCTGGGAAAAATCACGCTCAAATAGTTTAAGCCGCAGTACGGAAGATCAGTTTTTTATGTATTTACGGGTAAACACCCTAAATAAACTGGTACCCTACGC  >  NZ_CP009273/176609‑176766

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