Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 339,802 G→A E4K (GAA→AAA)  yahH → putative uncharacterized protein YahH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,8020GA100.0% 80.6 / NA 25E4K (GAA→AAA) yahHputative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/10);  total (15/10)

GCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGG  >  NC_000913/339724‑339878
                                                                              |                                                                            
gCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAaaaa                                                                            >  1:66013/1‑81 (MQ=255)
 cGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAaaaa                                                                            <  1:65980/80‑1 (MQ=255)
  gAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATc                                                                          >  2:66020/1‑81 (MQ=255)
   aaCGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCt                                                                         >  2:66032/1‑81 (MQ=255)
   aaCGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCt                                                                         <  2:65981/81‑1 (MQ=255)
      gTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTg                                                                        >  1:66028/1‑79 (MQ=255)
             cGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGa                                                               >  2:66033/1‑81 (MQ=255)
             cGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGa                                                               >  1:66019/1‑81 (MQ=255)
              gTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGa                                                               >  1:66026/1‑80 (MQ=255)
                   aCCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAg                                                           <  2:65982/79‑1 (MQ=255)
                   aCCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGc                                                          <  1:65983/80‑1 (MQ=255)
                       aCATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAgc                                                     >  1:66027/1‑81 (MQ=255)
                         aTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAgcg                                                    >  1:66024/1‑80 (MQ=255)
                              cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCa                                               <  1:65984/80‑1 (MQ=17)
                              cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAg                                              >  1:66040/1‑81 (MQ=17)
                                               gAACGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGc                             >  2:66036/1‑81 (MQ=255)
                                                aaCGCCTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGc                             <  1:65986/80‑1 (MQ=255)
                                                     cTTATCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATg                       >  2:66042/1‑81 (MQ=255)
                                                         tCCGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATgcggc                   <  2:65987/81‑1 (MQ=255)
                                                           cGGCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGt                 <  2:65988/81‑1 (MQ=255)
                                                             gCCTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTg                >  2:66039/1‑80 (MQ=255)
                                                               cTACGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAAc             <  2:65989/81‑1 (MQ=255)
                                                                  cGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGcc          >  1:66031/1‑81 (MQ=255)
                                                                  cGAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGcc          <  2:65990/81‑1 (MQ=255)
                                                                   gAATGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCt         >  2:66041/1‑81 (MQ=255)
                                                                    aaTGGCGCAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCtt        <  1:65991/81‑1 (MQ=37)
                                                                          gcAAAAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCgg  >  2:66054/1‑81 (MQ=17)
                                                                              |                                                                            
GCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGG  >  NC_000913/339724‑339878

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