Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,062,210	T→C	G390G (GGT→GGC)	agp →	bifunctional glucose‑1‑phosphatase/inositol phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,062,210	0	T	C	95.7%	65.8 / ‑5.0	23	G390G (GGT→GGC)	agp	bifunctional glucose‑1‑phosphatase/inositol phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (14/8);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTATCAGAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAG  >  NZ_CP009273/1062126‑1062285                                                                                     |                                                                            gtATCAGAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGccc                                                                        >  1:356859/1‑90 (MQ=255)      agagTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGCGTGACGCTGGAATTAAGCGGCTGCCCGATAg                                                                   <  1:594109/90‑1 (MQ=255)       gagTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGa                                                                  >  1:695532/1‑90 (MQ=255)       gagTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGa                                                                  >  1:639685/1‑90 (MQ=255)         gTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACg                                                                >  1:114376/1‑90 (MQ=255)         gTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACg                                                                >  1:410439/1‑90 (MQ=255)           gCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCt                                                              >  1:2349/1‑90 (MQ=255)           gCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCt                                                              >  2:440594/1‑90 (MQ=255)              gAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGAt                                                           <  2:410439/90‑1 (MQ=255)                  aGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGtt                                                       >  1:476252/1‑90 (MQ=255)                             gCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATg                                            >  2:85041/1‑90 (MQ=255)                                 aTGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                        <  2:509503/90‑1 (MQ=255)                                 aTGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                        <  2:127871/90‑1 (MQ=255)                                 aTGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                        >  2:607453/1‑90 (MQ=255)                                 aTGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATa                                        <  2:476252/90‑1 (MQ=255)                                     gTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGtt                                    >  2:380596/1‑90 (MQ=255)                                             tGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCg                            <  2:14510/90‑1 (MQ=255)                                                    aCCTGCGCAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTtgaa                     >  2:129023/1‑90 (MQ=255)                                                        gcgcAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTtgaatgaa                 <  2:61851/90‑1 (MQ=255)                                                          gcAGCGTGTGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGc               <  2:94556/90‑1 (MQ=255)                                                               gtgtGACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGa          >  1:271509/1‑90 (MQ=255)                                                                  tgACGCTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAAt       <  2:356859/90‑1 (MQ=255)                                                                       cTGGAATTAAGCGGCTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAg  >  2:698763/1‑90 (MQ=255)                                                                                     |                                                                            GTATCAGAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACAG  >  NZ_CP009273/1062126‑1062285

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