Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,717,498 G→C G374A (GGC→GCC)  ydhK → FUSC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,717,4980GC92.3% 34.6 / ‑2.8 13G374A (GGC→GCC) ydhKFUSC family protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (5/7);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTA  >  NZ_CP009273/1717414‑1717573
                                                                                    |                                                                           
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                        <  1:54483/90‑1 (MQ=255)
ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                        <  2:513202/90‑1 (MQ=255)
        ggCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTAtt                                                                >  1:93794/1‑90 (MQ=255)
            tGGCTCGTCATACCGACAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCtt                                                            <  2:309227/90‑1 (MQ=255)
             ggCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTCGATGATGATTGCCGCATGGAGTATTGCTTc                                                           >  2:563858/1‑90 (MQ=255)
                      tACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGa                                                  <  2:82831/90‑1 (MQ=255)
                         cGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGc                                               >  2:209266/1‑90 (MQ=255)
                                                gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTAt                                                                 <  1:576769/49‑1 (MQ=255)
                                                gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTAt                                                                 >  2:576769/1‑49 (MQ=255)
                                                 cGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                       <  1:255558/90‑1 (MQ=255)
                                                  gggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgca                      >  2:141016/1‑90 (MQ=255)
                                                               tttGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCg         <  1:491688/90‑1 (MQ=255)
                                                                      gTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactcta  >  2:70835/1‑90 (MQ=255)
                                                                                    |                                                                           
GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTA  >  NZ_CP009273/1717414‑1717573

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