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.80.0% 19.5 / 1.1 10coding (165/459 nt)ybbJNfeD family protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (3/5);  total (3/7)
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.51e-01

TGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCA  >  NZ_CP009273/510063‑510229
                                                                                        |                                                                              
tgtgCTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGca                                                                               <  2:697415/90‑2 (MQ=255)
 gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                              >  1:651006/1‑88 (MQ=255)
 gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                              >  2:719951/1‑88 (MQ=255)
  tgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcact                                                                             >  2:775292/1‑87 (MQ=255)
                 ccccGCCGCGACAACCATTTCCCCCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCa                                                             <  2:382394/89‑1 (MQ=255)
                  cccGCCGCGACAACCATTTCCCCCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCAc                                                            <  2:106256/90‑1 (MQ=255)
                              aCCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGaccac                                                <  2:453834/90‑1 (MQ=255)
                               ccATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGaccacc                                               <  1:151752/90‑1 (MQ=255)
                                        ccacAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACAACCAGGCCggga                                      >  2:201317/1‑88 (MQ=255)
                                                  gcggcgAGCAGCGTCAGGATGGCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTg                            >  2:292118/1‑90 (MQ=255)
                                                                 aGGATGGCAAACATTCCCCCTTGCCCCCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCcaca             <  2:369803/90‑1 (MQ=255)
                                                                  ggATGGCAAACATTCCCCCTTCCCCCCCCACCCCAGCGGCCCCAGCCAGACCCCCAGGCCGGTAATCCCCGCTGCCACGCCGCTCCacaa            <  2:27171/90‑1 (MQ=255)
                                                                      ggCAAACATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa        <  1:60464/90‑1 (MQ=255)
                                                                            cATTACCCCTTGCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCa  >  2:653968/1‑90 (MQ=255)
                                                                                        |                                                                              
TGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCA  >  NZ_CP009273/510063‑510229

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