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.90.0% 29.6 / ‑3.7 10coding (165/459 nt)ybbJNfeD family protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (8/1);  total (8/2)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAT  >  NZ_CP009273/510064‑510224
                                                                                       |                                                                         
gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                         >  1:171074/1‑88 (MQ=255)
            gcgcACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGc                                                            >  2:509523/1‑90 (MQ=255)
                                     caccacAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCgg                                   >  2:184301/1‑90 (MQ=255)
                                        cacaGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTaa                                >  2:37299/1‑90 (MQ=255)
                                          caGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATc                              >  2:752983/1‑90 (MQ=255)
                                           aGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCa                             >  2:257311/1‑90 (MQ=255)
                                               aGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGc                         >  2:632196/1‑90 (MQ=255)
                                                  cggcgAGCAGCGTCAGGATGGCAAACATTCCCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGc                      <  1:632196/90‑1 (MQ=255)
                                                          agcGTCAGGAGGGCAAACATTCCCCCTTCCACCCCCACCCCAGCGGCACCAGCCAGACCCCCAGGCCGGTAATCACCGCTGCCAcgccgc              <  1:394545/90‑1 (MQ=255)
                                                                  gATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCacaac      >  2:314000/1‑90 (MQ=255)
                                                                      gCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAt  >  1:20591/1‑90 (MQ=255)
                                                                                       |                                                                         
GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAT  >  NZ_CP009273/510064‑510224

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