Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,230,399 C→T G504G (GGC→GGT ygjK → glycoside hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,230,3990CT100.0% 100.0 / NA 32G504G (GGC→GGTygjKglycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/19);  total (13/19)

GCCGGAATATGGCGCGACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGCCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATAAATATGTCGCTAATGGCGGCAAACGTAG  >  NC_000913/3230267‑3230530
                                                                                                                                    |                                                                                                                                   
gCCGGAATATGGGGCGACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACg                                                                                                                              <  2:192708/140‑1 (MQ=255)
gCCGGAATATGGCGCGACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACg                                                                                                                              <  2:803018/140‑1 (MQ=255)
 ccGGAATATGGCGCGACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGa                                                                                                                             <  1:720323/140‑1 (MQ=255)
 ccGGAATATGGCGCGACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGa                                                                                                                             <  1:696913/140‑1 (MQ=255)
 ccGGAATATGGCGCGACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGa                                                                                                                             <  1:674889/140‑1 (MQ=255)
              cgACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAAtt                                                                                                                >  1:78514/1‑140 (MQ=255)
              cgACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAAtt                                                                                                                >  2:779041/1‑140 (MQ=255)
              cgACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTCTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAAtt                                                                                                                <  2:593327/140‑1 (MQ=255)
               gACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTc                                                                                                               >  1:191874/1‑140 (MQ=255)
                      gACAAAGCCCACAACTCTGAGACCGTCGAAATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGcaca                                                                                                        <  1:361451/140‑1 (MQ=255)
                      gACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGcaca                                                                                                        <  1:103122/140‑1 (MQ=255)
                               cacaACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTtgctgc                                                                                               <  2:78514/140‑1 (MQ=255)
                                           agCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATc                                                                                   <  2:158465/140‑1 (MQ=255)
                                           agCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATc                                                                                   >  2:798421/1‑140 (MQ=255)
                                            gcggcgAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCg                                                                                  <  1:523616/140‑1 (MQ=255)
                                                              cGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATgccgcc                                                                <  2:509285/140‑1 (MQ=255)
                                                               ggTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTACGGCACAGATTGATGCCACGTGGGAGTCGGGACGTGATGATGccgccg                                                               >  2:404012/1‑140 (MQ=255)
                                                                   aaaaaaGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGAACAGGTTGATGCGTCGTGGGAATCGGGGCGCGATGATGCCGCCGTCtt                                                           >  2:270362/1‑140 (MQ=255)
                                                                                      agaCGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGaa                                        <  2:158431/140‑1 (MQ=255)
                                                                                        aCGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACa                                      >  2:478847/1‑140 (MQ=255)
                                                                                         cGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAg                                     <  1:478847/140‑1 (MQ=255)
                                                                                         cGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAg                                     <  1:779041/140‑1 (MQ=255)
                                                                                            aGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTg                                  <  1:798421/140‑1 (MQ=255)
                                                                                             gTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTgg                                 >  1:748442/1‑140 (MQ=255)
                                                                                             gTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCATGGGAATCGGGGCGTGATGAAGCCGCCGTCTTAGGATTTATCAACAAAGAACAGTTgg                                 >  1:258893/1‑140 (MQ=255)
                                                                                                tGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATa                              <  2:748442/140‑1 (MQ=255)
                                                                                                              aCGCCCGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACCAAGCACTGTTGGATAAATATGTCGCTAAt                >  1:292950/1‑140 (MQ=255)
                                                                                                                 cccGCGTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATAAATATGTCGCTAATggc             >  1:3168/1‑140 (MQ=255)
                                                                                                                    gcgTGGTGGAGAAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATAAATATGTCGCTAATggcggc          >  1:32961/1‑140 (MQ=255)
                                                                                                                     cgTGGTGGAGAAAGGTCAGTACAACCGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATAAATATGTCGCTAATGGCGGCa         <  1:270362/140‑1 (MQ=255)
                                                                                                                          tggAGAAAGGTCAGTCCGACAGGCTGGAATTTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATAAATATGTCGCTAATGGCGGCAAACGt    <  1:423018/140‑1 (MQ=255)
                                                                                                                            gagaAAGGTCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATAAATATGTCGCTAATGGCGGCAAACGTAg  >  2:430426/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GCCGGAATATGGCGCGACCCGCGACAAAGCCCACAACACTGAGAGCGGCGAGATGCTGTTTACGGTGAAAAAAGGCGATAAAGAAGAGACGCAGTCTGGGCTGAACAACTACGCCCGCGTGGTGGAGAAAGGCCAGTACGACAGTCTGGAAATTCCGGCACAGGTTGCTGCATCGTGGGAATCGGGGCGTGATGATGCCGCCGTCTTTGGGTTTATCGACAAAGAACAGTTGGATAAATATGTCGCTAATGGCGGCAAACGTAG  >  NC_000913/3230267‑3230530

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