Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,787,923 +C intergenic (+165/‑55) cafA2 → / → moaA1 Ribonuclease G and E/molybdopterin biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,787,9201.C100.0% 29.3 / NA 9intergenic (+162/‑58)cafA2/moaA1Ribonuclease G and E/molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (4/5);  total (4/5)

ATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGA  >  AM260479/2787835‑2788032
                                                                                      |                                                                                                                
aTGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGccca                                                                   <  2:999156/134‑1 (MQ=255)
                             ggTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTc                                      >  2:260654/1‑134 (MQ=255)
                                    cACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGAtttt                               <  2:324546/134‑1 (MQ=255)
                                                              cagAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTg                        <  1:688345/115‑1 (MQ=255)
                                                              cagAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTg                        >  2:688345/1‑115 (MQ=255)
                                                                    gCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCcaca                                                             <  1:562854/72‑1 (MQ=255)
                                                                    gCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCcaca                                                             >  2:562854/1‑72 (MQ=255)
                                                                         cGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGa  <  1:507327/126‑1 (MQ=255)
                                                                         cGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGa  >  2:507327/1‑126 (MQ=255)
                                                                                      |                                                                                                                
ATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGA  >  AM260479/2787835‑2788032

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