Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,842,591 T→C L151L (TTA→TTG hycI ← protease involved in processing C‑terminal end of HycE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,842,5910TC100.0% 86.5 / NA 26L151L (TTA→TTGhycIprotease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/11);  total (15/11)

CGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCG  >  NC_000913/2842461‑2842726
                                                                                                                                  |                                                                                                                                       
cGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTacgc                                                                                                                                 >  2:435578/1‑139 (MQ=255)
    cGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAg                                                                                                                             <  1:518753/139‑1 (MQ=255)
                 cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATttcctt                                                                                                                >  2:162142/1‑139 (MQ=255)
                              acgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAg                                                                                                   <  2:316058/138‑1 (MQ=255)
                               cgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGt                                                                                                  <  2:74484/139‑1 (MQ=255)
                                                 aTTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCt                                                                                <  1:110229/139‑1 (MQ=255)
                                                                  gtgcACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCTCCAACTGCTCTATCCCCCAATTTCCTTCTCAGCCTCCCAGTCGTTGATAAGGGGTTCCACCGGCATCTTTATTCgc                                                               >  1:160523/4‑138 (MQ=255)
                                                                    tCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCt                                                             <  2:551835/139‑1 (MQ=255)
                                                                       ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCGAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTggg                                                          >  1:396321/1‑139 (MQ=255)
                                                                       ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGt                                                                                                  <  1:478047/99‑1 (MQ=255)
                                                                       ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGt                                                                                                  >  2:478047/1‑99 (MQ=255)
                                                                                      cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTgc                                                                                                                                   <  1:541380/51‑1 (MQ=255)
                                                                                      cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTgc                                                                                                                                   >  2:541380/1‑51 (MQ=255)
                                                                                      cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAgccgcc                                                                                                                        >  1:251460/1‑62 (MQ=255)
                                                                                      cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAgccgcc                                                                                                                        <  2:251460/62‑1 (MQ=255)
                                                                                           tCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCcg                                      >  1:391219/1‑139 (MQ=255)
                                                                                                aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAg                                                                                                           <  2:491929/65‑1 (MQ=255)
                                                                                                aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAg                                                                                                           >  1:491929/1‑65 (MQ=255)
                                                                                                aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                                 >  1:586034/1‑139 (MQ=255)
                                                                                                aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                                 >  1:564106/1‑139 (MQ=255)
                                                                                                aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGCGTCATCGGGTAGTAAAAGCCGACGat                                 >  2:469279/1‑139 (MQ=255)
                                                                                                     ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCTTCGGgtagta                                             <  2:141512/122‑1 (MQ=255)
                                                                                                     ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGgtagta                                             >  1:141512/1‑122 (MQ=255)
                                                                                                                atactCCTCCTCCACCGCCATCTGCGCGAAGCCGCCCTTTACTTCCCGGCCTTCCAGTAGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCg                 <  2:160523/138‑1 (MQ=255)
                                                                                                                             aCCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACtt    >  2:135942/1‑139 (MQ=255)
                                                                                                                               cGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCg  >  2:100965/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                       
CGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCACTTCG  >  NC_000913/2842461‑2842726

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