Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_LN832404 4,023,947 G→A 33.2% intergenic (+246/+135) 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_LN8324044,023,9470GA33.2% 26.7 / 11.2 18intergenic (+246/+135)EK12AG100_RS19960/fadANAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base A (3/3);  total (10/8)
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

AGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACC  >  NZ_LN832404/4023665‑4024240
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                     
aGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCGGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGtt                                                                                                                                                                                                                                                                                         <  1:812330/298‑1 (MQ=18)
     gCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGttcccgc                                                                                                                                                                                                                                                                                    >  1:129678/1‑293 (MQ=18)
           ttGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGttccagccccatc                                                                                                                                                                                                                                                                              >  2:357560/1‑287 (MQ=18)
                 aTGCGTTTGCAGTTATCGGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                                                                                                                                                                                        <  1:762389/298‑1 (MQ=17)
                         gCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaa                                                                                                                                                                                                                                                                >  2:94154/1‑295 (MQ=17)
                                 cTGAGAAATAAAAAAACCCGCCCCTGTCAGGCGGGAAGAACGGAAACTAAACTGTTATGAAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCGTATCGTGACTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGttaaacccgctc                                                                                                                                                                                                                                                        <  1:1020169/298‑12 (MQ=14)
                                   gagaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGACCGGTGGCCCGGCTCAGGCGTTCACGCCGCACCCGGCACGTGGttaaacccgcccaa                                                                                                                                                                                                                                                      >  2:10801/1‑285 (MQ=17)
                                    agaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCACGTGGttaaacccgctccca                                                                                                                                                                                                                                                     >  2:82474/1‑284 (MQ=14)
                                                                                                                                                                                                                     tatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTc                                                                   <  1:95694/285‑1 (MQ=255)
                                                                                                                                                                                                                      atctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCACACTGAACCTCTTTGCGTTCCACCAGACTcc                                                                  >  2:873689/13‑297 (MQ=255)
                                                                                                                                                                                                                                  atgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCCAACTGAACGTCTTTGCGTCCCATCAGATTCCGCAGCGTGTTGc                                                      >  1:962196/1‑298 (MQ=255)
                                                                                                                                                                                                                                        gATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATAc                                                >  1:644495/1‑298 (MQ=255)
                                                                                                                                                                                                                                         aTGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATAcg                                               <  1:422409/298‑1 (MQ=255)
                                                                                                                                                                                                                                          tGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATAcgc                                              <  1:357560/298‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGa                                      >  1:506789/1‑298 (MQ=255)
                                                                                                                                                                                                                                                               cGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATCCGCGCCCCGGAACAACCCCCCGGa                         >  2:419310/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                gTGCCTACAACTAGGCCGAACCGTAGGCCGGATAAGGCGTTGACGCCGCAGCCGGCAAGTGGTGCATGATGCCTGATGGGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGAt                        <  2:388231/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                gTGCCTACAAATAGGCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGAt                        <  1:44498/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 cGAACCGTAGGCCGGATAAGGCGTTTACTCCGCATCCGGCAATTGGGGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTGCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATcgc        <  1:473386/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        gttgggCGGGGTAAGTCGTCTAGGTCGCTGCCGGGAGTGGCTGGATGGTGCGGGATGGGACGGTGGCTCGACTGATCAGGCCTACATTGGGGCCGGGTCGGTAGACCTGAGAAGGCGTCCCCGCCTCATCCGGCAAGTGGTTAAACCCGCGCAAACACCGGCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAGCGTCTGTGCGTTCCAACAGGTTCAGCAGCGTGGGGCGGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCAc   <  2:324485/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGACCGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATAGCCGCACCGGAACAACCCAGCGGCCGACCCAGCCCCATCGCGCcccc  >  2:740710/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                         gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCAcc  >  1:143014/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                     
AGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACC  >  NZ_LN832404/4023665‑4024240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: