Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,228,195 G→C G132A (GGC→GCC)  ygjJ → protein YgjJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,228,1950GC100.0% 83.7 / NA 24G132A (GGC→GCC) ygjJprotein YgjJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (1/23);  total (1/23)

TGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGGCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACTGCGCTGAAGGTACAGAA  >  NC_000913/3228059‑3228327
                                                                                                                                        |                                                                                                                                    
tGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCAc                                                                                                                                   <  1:320974/140‑1 (MQ=255)
  gTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCg                                                                                                                                 <  1:171829/140‑1 (MQ=255)
  gTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCg                                                                                                                                 <  1:245607/140‑1 (MQ=255)
        ccTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATTGTACCAAAGATGCCAccccaaaca                                                                                                                           >  2:121988/1‑133 (MQ=255)
                aCGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGt                                                                                                                   <  1:697142/140‑1 (MQ=255)
                             gACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGAtt                                                                                                      <  2:625601/140‑1 (MQ=255)
                                       gCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATTGTTCCCAAGATGGCCCCCAAAACATGCTCAGTCTGCGTCAGGGTTTCTCCCTggg                                                                                            >  1:493082/1‑139 (MQ=255)
                                        ccGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGGTATTTTTCCAAAAAAGGCACCCAAAAAAAACTCAGTCTCCGGCCGGAATGCTccccgggt                                                                                           >  2:563211/1‑133 (MQ=255)
                                            ccACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGAcc                                                                                       <  1:646863/140‑1 (MQ=255)
                                                      aaTGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAAtggtg                                                                             <  1:644661/140‑1 (MQ=255)
                                                         gCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTgggg                                                                          <  1:182062/140‑1 (MQ=255)
                                                          ccTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGc                                                                         <  1:182792/140‑1 (MQ=255)
                                                               gagacAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGAttt                                                                    <  1:275443/136‑1 (MQ=255)
                                                                    aGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCAccc                                                               <  1:706277/140‑1 (MQ=255)
                                                                      aCCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCAccccc                                                             <  1:679228/140‑1 (MQ=255)
                                                                                gggTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCt                                                   <  2:193492/140‑1 (MQ=255)
                                                                                   tggtTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACg                                                <  1:630321/140‑1 (MQ=255)
                                                                                    ggtTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGt                                               <  1:258017/140‑1 (MQ=255)
                                                                                      tttCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCa                                             <  2:628899/140‑1 (MQ=255)
                                                                                                    ttAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCcgc                               <  2:84417/140‑1 (MQ=255)
                                                                                                    ttAACCGTTTCTATGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCcgc                               <  1:115176/140‑1 (MQ=255)
                                                                                                                aTGGCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACt                   <  1:735376/140‑1 (MQ=255)
                                                                                                                  ggCGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACTgc                 <  1:62047/140‑1 (MQ=255)
                                                                                                                    cGGTTATCGTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACTgcgc               <  1:121988/140‑1 (MQ=255)
                                                                                                                            gTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACTGCGCTGAAGGTa       <  2:3988/140‑1 (MQ=255)
                                                                                                                            gTACCAATGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACTGCGCTGAAGGTa       <  2:493082/140‑1 (MQ=255)
                                                                                                                                 aaTGATGCCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACTGCGCTGAAGGTACAGaa  >  2:341162/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATGCCTGGGTGCCAGACCATAGCTGGGGTGGTTTCGAGGGCGGAATTAACCGTTTCTATGGCGGTTATCGTACCAATGATGGCACCGAAATCATGCTCAGTCTGCGTCAGGATTCCTCGCTGGATGACCTGCAATGGTGGGGCGATTTCACCCCCGATCTGGGCTACGTCATTCCCAATACCCGCGACATTATGACTGCGCTGAAGGTACAGAA  >  NC_000913/3228059‑3228327

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