Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,473,944 (G)6→7 coding (550/1758 nt) ynbC → bifunctional alpha/beta hydrolase/class I SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,473,9381.G95.7% 75.6 / ‑3.2 23coding (544/1758 nt)ynbCbifunctional alpha/beta hydrolase/class I SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (13/9);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGGCGTTATGGCATCGTCTGCGTGGTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGT  >  NZ_CP009273/1473854‑1474016
                                                                                      |                                                                              
gcgcTGGCGTTATGGCATCGTCTGCGTGGTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACggggg                                                                           >  1:404381/1‑90 (MQ=255)
     ggCGTTATGGCATCGTCTGCGTGGTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACgggggggcga                                                                      >  2:681359/1‑86 (MQ=255)
                         gTGGTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTg                                                   <  1:681359/90‑1 (MQ=255)
                            gTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGa                                                          <  1:71977/80‑1 (MQ=255)
                            gTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGa                                                          >  2:71977/1‑80 (MQ=255)
                            gTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATc                                                >  2:301240/1‑90 (MQ=255)
                              cTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATcac                                              >  1:486542/1‑90 (MQ=255)
                                          aaTTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGc                                  >  1:503942/1‑90 (MQ=255)
                                          aaTTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGc                                  >  2:809155/1‑90 (MQ=255)
                                          aaTTCCTATGTGAAAGGACAGTTTTTGCCCCCCGATCGGGAATGGGGGGGCGAGCTTCAATAATGTTCCGCTGATGACACGGGCGATTGc                                  <  2:293627/90‑1 (MQ=255)
                                                      aaaGGACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCtt                      >  2:636054/1‑90 (MQ=255)
                                                         ggACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCt                   >  1:735091/1‑90 (MQ=255)
                                                         ggACGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCt                   <  2:363657/90‑1 (MQ=255)
                                                           aCGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATc                                                >  2:795696/1‑59 (MQ=255)
                                                           aCGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATc                                                <  1:795696/59‑1 (MQ=255)
                                                           aCGCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCg                 >  1:922169/1‑90 (MQ=255)
                                                             gCTATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGAt               >  2:166903/1‑90 (MQ=255)
                                                               tATTTGACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACCCAGACGATTGCCGTTAATATCTTGCTCGAtct             <  2:486542/90‑1 (MQ=255)
                                                                    gACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCg                                                      <  1:964694/44‑1 (MQ=38)
                                                                    gACCCACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCg                                                      >  2:964694/1‑44 (MQ=38)
                                                                      cccACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACaaa      <  1:994174/90‑1 (MQ=255)
                                                                      cccACGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACaaa      <  2:735091/90‑1 (MQ=255)
                                                                          cGATCGGCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGt  >  1:679274/1‑90 (MQ=255)
                                                                               ggCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTg                                   <  1:933394/52‑1 (MQ=255)
                                                                               ggCAACGGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTg                                   >  2:933394/1‑52 (MQ=255)
                                                                                      |                                                                              
GCGCTGGCGTTATGGCATCGTCTGCGTGGTCTGTTTTTTATTAATTCCTATGTGAAAGGACGCTATTTGACCCACGATCGGCAACGGGGGGCGAGTTTCAATAATGATCCGCTGATCACACGGGCGATTGCCGTTAATATCTTGCTCGATCTCTACAAAACGT  >  NZ_CP009273/1473854‑1474016

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