Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,390,883 Δ1 bp coding (272/1470 nt) rng ← ribonuclease G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,390,8830C.90.9% 35.5 / ‑3.2 11coding (272/1470 nt)rngribonuclease G
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (7/3);  minor base G (0/1);  total (7/4)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCG  >  NZ_CP009273/3390805‑3390960
                                                                              |                                                                             
gcACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGGCACCGTGAAtt                                                                   >  2:604461/1‑90 (MQ=255)
                  ggttCTTTCCCCCCCTGCACCATCAGATCTGGCCCCTGACGACCCAGTCCCGAGATGTGGGCGCCGTGAATTGCTTTTTTTCTCCCCCCg                                                 <  2:286493/87‑1 (MQ=255)
                                  gCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGGCACCGTGAATTGCttt                                                              >  1:594005/1‑61 (MQ=255)
                                  gCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGGCACCGTGAATTGCttt                                                              <  2:594005/61‑1 (MQ=255)
                                  gCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGtgtg                                 >  2:593101/1‑90 (MQ=255)
                                   ccccATAAGTTTTTCCCCCTGCCGACCCAGTCCCGAGATGGCGGCCCGGTAATTTGTTTTTTTTCTTCCCCCGCCACACATTCGGTGTGc                                <  1:109612/88‑1 (MQ=255)
                                                tGCCCCTGACGAACCAGTTCCGAGATGTCGGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCgg                   >  2:216182/1‑90 (MQ=255)
                                                     cTGACGAACCAGTTCCGAGATGTCGGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCa              >  1:444298/1‑90 (MQ=255)
                                                     cTGACGAACCAGTTCCGAGATGTCGGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCa              >  1:486103/1‑90 (MQ=255)
                                                                gTTCCGAGTTGTCGGCCCCGTAATTTTTTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcgc   <  2:27988/90‑1 (MQ=255)
                                                                 ttCCGAGATGTCGGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcgcg  >  1:274744/1‑90 (MQ=255)
                                                                              |                                                                             
GCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCG  >  NZ_CP009273/3390805‑3390960

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