Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,594,128 Δ1 bp coding (326/957 nt) cR ← Recombinational DNA repair protein (RecF pathway)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,594,1280C.100.0% 62.6 / NA 15coding (326/957 nt)cRRecombinational DNA repair protein (RecF pathway)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (9/6);  total (9/6)

GCTGGGCGGACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCTTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGCGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGTC  >  AM260479/2594027‑2594252
                                                                                                     |                                                                                                                            
gcTGGGCGGACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGccc                                                                                              >  1:822157/1‑133 (MQ=255)
        gACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATc                                                                                     >  2:1012569/1‑134 (MQ=255)
        gACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATc                                                                                     >  2:935735/1‑134 (MQ=255)
                  cGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCt                                                                           >  1:175958/1‑134 (MQ=255)
                  cGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCt                                                                           >  1:224132/1‑134 (MQ=255)
                   gACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCtt                                                                          <  2:175958/134‑1 (MQ=255)
                                             cGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 <  1:1012569/133‑1 (MQ=255)
                                             cGACCAGTGCCTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 <  2:224132/133‑1 (MQ=255)
                                                                            tGCCTTCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAg                 <  2:517106/134‑1 (MQ=255)
                                                                                ttCACAGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCt             >  1:842596/1‑134 (MQ=255)
                                                                                   acaGTGGCGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGTGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTc           <  2:842596/133‑1 (MQ=255)
                                                                                          cGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 >  1:537092/1‑88 (MQ=255)
                                                                                          cGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTg                                                 <  2:537092/88‑1 (MQ=255)
                                                                                          cGTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGt   >  1:946889/1‑134 (MQ=255)
                                                                                           gTTCAGAACGGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGTc  >  1:437780/1‑134 (MQ=255)
                                                                                                     |                                                                                                                            
GCTGGGCGGACTTGGGCCCGACGCCGGGCAGCACCCGCAGCGCTTCGACCAGTGCTTGCAGCGAGGTCGGGGACGATGCCTTCACAGTGGCGTTCAGAACGCGCAGCTTGAAGCCCGGGGGCAGCGGCAGGCCCGAGGTCATCGAGCCCATCTTTTCCTGCGAGGTGGCTTCGGCCTTGCGCACGGCGTCATTGAAGGCGGCGGCGACCAGGTCTTCCAGCAGGTC  >  AM260479/2594027‑2594252

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