Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_LN832404 4,023,962 T→C 39.1% intergenic (+261/+120) 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,9620TC39.1% 23.1 / 19.2 23intergenic (+261/+120)EK12AG100_RS19960/fadANAD(P)H‑flavin reductase/acetyl‑CoA C‑acyltransferase FadA
Reads supporting (aligned to +/- strand):  ref base T (11/3);  new base C (4/5);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATC  >  NZ_LN832404/4023676‑4024253
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                   
ttGGCGAAGCGTTTGCATTTATCTGAGAAATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCAcgc                                                                                                                                                                                                                                                                                                 <  1:366242/292‑1 (MQ=18)
  ggCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCAcgccg                                                                                                                                                                                                                                                                                               >  1:542120/1‑292 (MQ=18)
       tGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCACGCCGCATcc                                                                                                                                                                                                                                                                                          <  1:136422/292‑1 (MQ=18)
       tGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCACGCCGCATcc                                                                                                                                                                                                                                                                                          >  1:299922/1‑292 (MQ=18)
          tttgttaTTTAGCTTAGTGATAAAAAAACCCGCCCCTTGCAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGAGCTATGACGTAGCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCGTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCACGCCGCATCCGGc                                                                                                                                                                                                                                                                                       <  2:780562/286‑1 (MQ=14)
                aTTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATCAGGCGTTTACGCCGCCTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTACCAGGCCTAC‑CTTGGTGCCGGATCGGTAGACCGCCTAA‑GGCGATCACGCCGAATCCGGCaggtgg                                                                                                                                                                                                                                                                                 >  1:607211/1‑287 (MQ=17)
                 tttATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCACGCCGCATCCGGCAAGTGGt                                                                                                                                                                                                                                                                                >  1:543537/1‑292 (MQ=17)
                  ttATCGGAGATAGAAAAAGACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGGTATTCAGTGGCATTTAGATCTATTACGTAGCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACTCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCACGCCGCATCCGGCAAGTGGtt                                                                                                                                                                                                                                                                               <  2:441864/292‑2 (MQ=14)
                        gagaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAC‑ATTGGTGCCGGATCGGTAGACCGGATAA‑GGCGTTCACGCCGCATCCGGCAAGTGGttaaaccc                                                                                                                                                                                                                                                                         <  2:234276/292‑8 (MQ=14)
                                                                                                                                                                                                         gtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTCAACCCGCTCCAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCCGACCAAACTGAACGCCTTTGCGTTCCAt                                                                                       >  1:328008/15‑292 (MQ=255)
                                                                                                                                                                                                         gtatctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGGTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGCCTTTGCGTTCCAt                                                                                       >  1:456412/15‑292 (MQ=255)
                                                                                                                                                                                                             ctggcaaaccatgCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCCCAAACACCGTCGCAATACCCTGACCCACAACGAGACACACCCTCCGCAGCCCAAACCGACCGTCCTTGCGGTCCATCaaa                                                                                   >  2:323647/11‑290 (MQ=18)
                                                                                                                                                                                                                             gATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATCCACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGc                                                                   >  1:702934/1‑292 (MQ=255)
                                                                                                                                                                                                                                gctgcGCTGTGTCGTCTGTGCGTGCCTACAAATAGT‑CC‑GAACCGTGGGCCGGATAA‑GGCGGTTGCGCCGCATCCGGCAATTGGTGCATGATGCCTGAGGCGACGCTGGCGCGTCTTATCAGGCCTACATGGGGGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTTGTGCTGa                                                                <  2:456412/288‑1 (MQ=37)
                                                                                                                                                                                                                                 cGACGCTGTCGCGTCTTACCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCCACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCCTCCGGCAAGCGGTTAAACCCGCTCAAACACCCTCGCAATACCCCGACCCAGCCCGTTACACATCCCCGCCAGACCCAACTGAACGTCCTCCCCTTCCATCAGACTCAGCAGCGTGGGGCTGAt                                                               >  2:28836/1‑292 (MQ=25)
                                                                                                                                                                                                                                             gTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCCTGGTGCTGATACGCGCACCGGa                                                   >  1:161680/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                 ctggccggAACCGTAGG‑CAGATAAGGGCGTTTACGCCGGCGTCGGCAATTGGTGCATGGTGCCTGATGCGACGCGGGCGCGTCGTATCAGGCCTGCAGTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACTCCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAGCGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGAcc                                 <  2:607211/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCCATACCCTGACCCAGACCGATACACATCGTCTCCCGCCCAAACTGCACGTCTATGCGTTCCAGCAGATTACGCAGCGTGGCTCTTATACGCCCCCCGGAACACCACAGCGGATGACCCAg                              >  2:396279/1‑292 (MQ=37)
                                                                                                                                                                                                                                                                    ccc‑GAACCGGAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGAGCATGATGCCTGGTGCGACGCTGGCGCGTCTTAACAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAgc                             <  1:16518/291‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaCCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCACCGTGGTTCTGATACGCGCACCGGACCAACCCCGCGGATCACCCCACGCGa                         >  2:92861/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                            cGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCTGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCCGACCCATACCGATACACATCGACGCCAGACCAAACCGAACCTCCTTCCGCTCTATCAGATTCAGCAGCCTGGTACTGACACCACCCCCGGAACCACCCAGTCGATgcacggccgctatcg                      >  2:254162/1‑278 (MQ=25)
                                                                                                                                                                                                                                                                             gTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGCTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCCGATACGCGCACCGGAACAACCCAGCGGATGGCCCAGCGCGATccc                     >  1:368088/1‑290 (MQ=255)
                                                                                                                                                                                                                                                                                   cGGGTAA‑GGCGTTTAAGCCGGAGCCGGCAATTGGTGCATGATGGCTGATGCGGCGCTGGCGCGTCTTATCAGGCCTACATGGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACACCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCAcc               <  2:542120/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      aTAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACCGAACGTCTTTGCGTTCCATCAGATTCCGCAGCGTGGTGCTGATCCGCGCACCGGCACACACCAGCGGATGACCCCGCGCGATCGCGCCCCCGtt            >  2:216119/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                                           ggCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAgg        <  1:275213/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             cGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGtt      <  1:467044/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 tACGCCGGATCCGGGAATTGGTGCATGTTGCCTGATGGGACGCGGGCGCGTCTTATCAGGCGTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTGCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGAtc  <  2:702934/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                   
TTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGT‑CC‑GAACCGTAGGCCGGATAA‑GGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGCGGATGACCCAGCGCGATCGCGCCACCGTTGAGGTTGATC  >  NZ_LN832404/4023676‑4024253

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