Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,014,062 A→G D55G (GAC→GGC)  ygfM → molybdopterin‑dependent oxidoreductase FAD‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,014,0620AG93.3% 43.1 / ‑3.7 15D55G (GAC→GGC) ygfMmolybdopterin‑dependent oxidoreductase FAD‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (9/5);  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 = 9.18e-01

GGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTA  >  NZ_CP009273/3013989‑3014139
                                                                         |                                                                             
ggCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGAt                                                               <  1:144867/90‑1 (MQ=255)
                  acacCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGAtt                                             >  1:17883/1‑90 (MQ=255)
                   cacCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTg                                            >  1:398782/1‑90 (MQ=255)
                       aaCCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGgcgc                                        >  1:208827/1‑90 (MQ=255)
                       aaCCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGgcgc                                        >  2:36945/1‑90 (MQ=255)
                            gTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGt                                   >  1:220558/1‑90 (MQ=255)
                                   tAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCAc                                                     <  1:31364/65‑1 (MQ=255)
                                   tAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCAc                                                     >  2:31364/1‑65 (MQ=255)
                                        aGATTGCCATTTCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGc                       <  1:493907/90‑1 (MQ=255)
                                                 tttCCTTACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTg              >  2:424123/1‑90 (MQ=255)
                                                  ttCCTTACAGGATCTGGAACTGGGCTGGGTTGACt                                                                    <  1:253471/35‑1 (MQ=37)
                                                  ttCCTTACAGGATCTGGAACTGGGCTGGGTTGACt                                                                    >  2:253471/1‑35 (MQ=37)
                                                       tACAGGATCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATgcgc        <  1:424123/90‑1 (MQ=255)
                                                             aTCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTa  >  2:199777/1‑90 (MQ=255)
                                                             aTCTGGAACTGGGCTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTa  >  2:364266/1‑90 (MQ=255)
                                                                         |                                                                             
GGCAGCAAACTCAATGCTACACCAACCCGTACCGATAAAAAGATTGCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTA  >  NZ_CP009273/3013989‑3014139

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