Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,955,978 T→G Q49P (CAG→CCG)  ppdC ← prepilin‑type N‑terminal cleavage/methylation domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,955,9780TG100.0% 58.7 / NA 19Q49P (CAG→CCG) ppdCprepilin‑type N‑terminal cleavage/methylation domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (7/12);  total (7/12)

GCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCA  >  NZ_CP009273/2955892‑2956053
                                                                                      |                                                                           
gCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGtt                                                                          <  1:327380/90‑1 (MQ=255)
gCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGtt                                                                          >  2:218349/1‑90 (MQ=255)
       cTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTg                                                                   >  2:451850/1‑90 (MQ=255)
          cATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGtt                                                                >  2:519402/1‑90 (MQ=255)
              cGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTctg                                                            <  1:519402/90‑1 (MQ=255)
              cGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTctg                                                            <  2:442997/90‑1 (MQ=255)
              cGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTctg                                                            <  2:358514/90‑1 (MQ=255)
               ggTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTctgc                                                           >  1:99173/1‑90 (MQ=255)
                      cTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAg                                                    <  2:16493/90‑1 (MQ=255)
                        gCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGAc                                                  <  2:548032/90‑1 (MQ=255)
                        gCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGAc                                                  <  2:278083/90‑1 (MQ=255)
                        gCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTACTGGCAAGAc                                                  >  2:501747/1‑90 (MQ=255)
                              tGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGAACTGGTTTCTGCTGGCAAGACTGTTCa                                            >  2:373514/1‑90 (MQ=255)
                                              cgcgcgCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGaa                            <  2:25604/90‑1 (MQ=255)
                                                      gTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATa                    <  2:364466/90‑1 (MQ=255)
                                                               gctgcCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGa           <  2:135410/90‑1 (MQ=255)
                                                                ctgcCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATa                    >  2:522263/1‑80 (MQ=255)
                                                                ctgcCAGCCATGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATa                    <  1:522263/80‑1 (MQ=255)
                                                                        caTGCCGCCAGAGCGGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAAtcatca  <  2:159934/90‑1 (MQ=255)
                                                                                      |                                                                           
GCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTTCGCTGGAAACCCGATAACGCAGTGACAATCATCA  >  NZ_CP009273/2955892‑2956053

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