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% 25.6 / NA 8intergenic (+162/‑58)cafA2/moaA1Ribonuclease G and E/molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (5/3);  total (5/3)

CCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCT  >  AM260479/2787799‑2788009
                                                                                                                          |                                                                                         
ccACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTgg                                                                                <  1:324657/134‑1 (MQ=255)
         ctcCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCAcgccgc                                                                       >  1:580822/1‑134 (MQ=255)
             ttttttCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAgg                                                                   >  1:277011/1‑134 (MQ=255)
                                    aTGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGccca                                            <  2:144953/134‑1 (MQ=255)
                                                  tctcCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGa                              >  1:823923/1‑134 (MQ=255)
                                                                              cGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCt  >  2:550647/1‑134 (MQ=255)
                                                                                                      aGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGAcc   >  1:918871/1‑109 (MQ=255)
                                                                                                      aGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGAcc   <  2:918871/109‑1 (MQ=255)
                                                                                                                          |                                                                                         
CCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCT  >  AM260479/2787799‑2788009

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