Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,842,573 C→T *157* (TAG→TAA 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,5730CT100.0% 87.9 / NA 28*157* (TAG→TAAhycIprotease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/12);  total (16/12)

ATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGAT  >  NC_000913/2842446‑2842695
                                                                                                                               |                                                                                                                          
aTCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCttcttc                                                                                                                 >  1:251981/1‑139 (MQ=255)
           aTGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACt                                                                                                      >  2:421715/1‑139 (MQ=255)
              tCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTgcg                                                                                                   <  1:421715/139‑1 (MQ=255)
                  aCGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGaa                                                                                               <  2:251981/139‑1 (MQ=255)
                     aaCCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgcc                                                                                            >  1:416363/1‑139 (MQ=255)
                                cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                  >  1:519600/1‑106 (MQ=255)
                                cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                  <  2:519600/106‑1 (MQ=255)
                                cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttcctt                                                                                 >  2:675088/1‑139 (MQ=255)
                                         cgcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCtt                                                                        >  1:558018/1‑139 (MQ=255)
                                             gcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                  <  2:747402/93‑1 (MQ=255)
                                             gcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                  >  1:747402/1‑93 (MQ=255)
                                                       cAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                          <  2:97661/139‑1 (MQ=255)
                                                                    ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                             >  1:616/1‑139 (MQ=255)
                                                                                                 aCATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                        <  1:618221/131‑1 (MQ=255)
                                                                                                 aCATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                        >  2:618221/1‑131 (MQ=255)
                                                                                                     cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCt                              >  2:47781/1‑121 (MQ=255)
                                                                                                     cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCt                              <  1:47781/121‑1 (MQ=255)
                                                                                                     cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaaa            >  2:796978/1‑139 (MQ=255)
                                                                                                       cTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGc                               <  1:545014/118‑1 (MQ=255)
                                                                                                       cTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGc                               >  2:545014/1‑118 (MQ=255)
                                                                                                           ccTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCcga      <  2:616/139‑1 (MQ=255)
                                                                                                            cTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACg                                                       >  2:207237/1‑89 (MQ=255)
                                                                                                            cTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACg                                                       <  1:207237/89‑1 (MQ=255)
                                                                                                               aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCGGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat  >  2:107012/1‑139 (MQ=255)
                                                                                                               aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACg                                                       <  1:478905/86‑1 (MQ=255)
                                                                                                               aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGttt                                                   <  1:820959/90‑1 (MQ=255)
                                                                                                               aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGttt                                                   >  2:820959/1‑90 (MQ=255)
                                                                                                               aTTCCTTAATGAAAAATGACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACg                                                       >  2:478905/1‑86 (MQ=255)
                                                                                                                               |                                                                                                                          
ATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGAT  >  NC_000913/2842446‑2842695

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