Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 510,152 Δ1 bp coding (164/459 nt) ybbJ ← NfeD family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273510,1510C.83.3% 11.1 / ‑3.4 6coding (165/459 nt)ybbJNfeD family protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (4/1);  total (4/2)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCA  >  NZ_CP009273/510064‑510229
                                                                                       |                                                                              
gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                              >  1:263006/1‑88 (MQ=255)
  gCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcactc                                                                            >  1:695250/1‑86 (MQ=255)
  gCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcactc                                                                            >  2:259269/1‑86 (MQ=255)
                                   tCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGcc                                          <  1:98679/90‑1 (MQ=255)
                                                          agcGTCAGGATGGCAAACATTCCCCCTTCCCCCCCCACCCCAGCGGCCCCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgccgc                   <  1:512188/90‑1 (MQ=255)
                                                           gcGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCt                  >  1:316887/1‑90 (MQ=255)
                                                           gcGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCt                  >  2:592430/1‑90 (MQ=255)
                                                                     ggCAAACATTACCCCTTGCACTCCCAACCCAGGGGCACCAGCCAGACCACCAAGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa        >  1:447343/1‑90 (MQ=255)
                                                                           cATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCa  >  1:40602/1‑90 (MQ=255)
                                                                                       |                                                                              
GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCA  >  NZ_CP009273/510064‑510229

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