Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_LN832404	4,023,942	C→G	53.3%	intergenic (+241/+140)	EK12AG100_RS19960 → / ← fadA	NAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_LN832404	4,023,942	0	C	G	53.3%	4.4 / 23.6	15	intergenic (+241/+140)	EK12AG100_RS19960/fadA	NAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base G (5/3);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGAT  >  NZ_LN832404/4023670‑4024231                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                  gCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGttc                                                                                                                                                                                                                                                                               <  1:383580/294‑2 (MQ=18)  ccTGTGGGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGttca                                                                                                                                                                                                                                                                              <  1:549634/294‑3 (MQ=18)     gTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCT‑CCATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCAcgc                                                                                                                                                                                                                                                                           >  1:399649/1‑294 (MQ=18)             aTGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCATCTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                                                                                                   <  1:640829/294‑1 (MQ=17)                 gTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCaa                                                                                                                                                                                                                                                               >  1:115414/1‑294 (MQ=17)                         ttATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaa                                                                                                                                                                                                                                                       >  1:16710/1‑291 (MQ=17)                         ttATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaa                                                                                                                                                                                                                                                       >  1:91164/1‑291 (MQ=17)                                 gaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaaacccgctc                                                                                                                                                                                                                                               >  1:581241/1‑283 (MQ=14)                                                                                                                         aCCATGCCCGGGGCGCAGACGTTGCGTGTTATCGTGGCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGGCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGtt                                                                                                                                                      <  2:570800/294‑1 (MQ=25)                                                                                                                                                                                                                    ctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGAtt                                                           >  1:152755/11‑294 (MQ=255)                                                                                                                                                                                                                           acctcgCCGGATGCTCCGCTGTAGGGTCTTTCTGGGCCTACAAATATT‑CC‑GAAACTTGGTCCGGAGAAGGCATTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCagcagc                                                    <  2:516760/289‑1 (MQ=37)                                                                                                                                                                                                                                      tGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATAAACATCGTCGCCCGACCCAACTGAACGACTTTGCCTTCAACACGATTCACCAGCgtgggcctgcc                                         >  2:169521/1‑287 (MQ=25)                                                                                                                                                                                                                                       gcggcgCTGGCGTGTCTGATCGTGCCTACAAATAGT‑CC‑GAACGTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGGGCATGATGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATa                                        <  2:175119/294‑1 (MQ=255)                                                                                                                                                                                                                                         gACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGCCCAAACTGAACGTTCTTGCGTTCCATCAGACTCaacagccgggttctcatccc                                      >  2:71307/1‑280 (MQ=255)                                                                                                                                                                                                                                               ggcgcgTGTTATCGGGCCGAAAAAGAGT‑CC‑GAACCGGAGGCCGGATAAGGCGTTTACGCCGCAGCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCAcc                                <  2:115221/292‑1 (MQ=255)                                                                                                                                                                                                                                                        tATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCCCCGGACCAAcc                       >  1:200642/1‑294 (MQ=255)                                                                                                                                                                                                                                                             tgccgtcAAAGAGT‑CC‑GAGACGTTGGCCGGATAAGGCGGTCACGTCGGTTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCGACATGGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCg                  <  2:171210/288‑1 (MQ=38)                                                                                                                                                                                                                                                                         gT‑CC‑GAACCGCAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAgcg      >  1:50206/1‑294 (MQ=255)                                                                                                                                                                                                                                                                               gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGAt  <  1:128936/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                  GCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGAT  >  NZ_LN832404/4023670‑4024231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑