Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	886,189	A→G	G7G (GGT→GGC)	grxA ←	glutaredoxin 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	886,189	0	A	G	95.2%	58.6 / ‑3.2	21	G7G (GGT→GGC)	grxA	glutaredoxin 1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (11/9);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTG  >  NZ_CP009273/886108‑886265                                                                                  |                                                                             ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAAt                                                                      >  1:227951/1‑90 (MQ=255)   gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                    >  2:492909/1‑90 (MQ=255)   gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                    >  1:433517/1‑90 (MQ=255)       ctgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGt                                                                <  1:276038/90‑1 (MQ=255)           aaaTCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                        <  1:317802/78‑1 (MQ=255)           aaaTCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                        >  2:317802/1‑78 (MQ=255)           aaaTCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGc                                                            >  1:147034/1‑90 (MQ=255)                atcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                              >  2:558354/1‑83 (MQ=255)                atcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                              <  1:558354/83‑1 (MQ=255)                     gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                  <  1:223293/90‑1 (MQ=255)                             cTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                                          <  2:345223/90‑1 (MQ=255)                                    tctcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTAt                                   <  1:492909/90‑1 (MQ=255)                                       cAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                <  2:433517/90‑1 (MQ=255)                                           cAGATCTTTTGCACGCACACAGTAAGGGCAATCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTa                            >  2:444322/1‑90 (MQ=255)                                                cTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGAt                       <  2:530914/90‑1 (MQ=255)                                                  tttGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCa                     <  2:227951/90‑1 (MQ=255)                                                   ttGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCAc                    >  1:177598/1‑90 (MQ=255)                                                         gcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTaaa              >  1:262548/1‑90 (MQ=255)                                                          cacacaGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTaaaa             >  2:310847/1‑90 (MQ=255)                                                                    gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCt   >  1:152594/1‑90 (MQ=255)                                                                     ggCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTg  <  1:335501/90‑1 (MQ=255)                                                                                  |                                                                             CTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTG  >  NZ_CP009273/886108‑886265

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