Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 247,461 G→A *250* (TAG→TAA yafL → putative lipoprotein and C40 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913247,4610GA100.0% 84.5 / NA 27*250* (TAG→TAAyafLputative lipoprotein and C40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/13);  total (14/13)

TGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGAT  >  NC_000913/247325‑247575
                                                                                                                                        |                                                                                                                  
tgtATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCCGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGa                                                                                                                  >  2:70663/1‑139 (MQ=255)
                         cGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCg                                                                                        >  2:101551/1‑139 (MQ=255)
                                         ggCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGt                                                                        <  1:346825/139‑1 (MQ=255)
                                          gCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACg                                                                             >  1:459130/1‑133 (MQ=255)
                                          gCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACg                                                                             <  2:459130/133‑1 (MQ=255)
                                                            gCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgcg                                                     >  1:283587/1‑139 (MQ=255)
                                                             ccGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgcgc                                                    <  2:283587/139‑1 (MQ=255)
                                                                  tAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTAcgccg                                                                                              >  2:7333/1‑92 (MQ=255)
                                                                  tAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTAcgccg                                                                                              <  1:7333/92‑1 (MQ=255)
                                                                  tAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCgg                                                                                       <  2:766520/99‑1 (MQ=255)
                                                                  tAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCgg                                                                                       >  1:766520/1‑99 (MQ=255)
                                                                       gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                      <  2:918056/127‑1 (MQ=255)
                                                                       gAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                      >  1:918056/1‑127 (MQ=255)
                                                                          ccTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTGCGCAGGTGCCTGATGCGACGCTGgc                                                      >  1:636598/1‑124 (MQ=255)
                                                                          ccTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGgc                                                      <  2:636598/124‑1 (MQ=255)
                                                                                     ggACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCcgc                            >  2:305468/1‑139 (MQ=255)
                                                                                      gACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCcgcg                           >  1:631657/1‑139 (MQ=255)
                                                                                       aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTg                                                                                 >  1:881778/1‑84 (MQ=38)
                                                                                       aCCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTg                                                                                 <  2:881778/84‑1 (MQ=38)
                                                                                                tGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCg                 <  1:101551/139‑1 (MQ=255)
                                                                                                      cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTgg                                                       <  2:595528/95‑1 (MQ=255)
                                                                                                      cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGc           >  2:647863/1‑139 (MQ=255)
                                                                                                      cgcgcAGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCCGATGCGACGCTgg                                                       >  1:595528/1‑95 (MQ=37)
                                                                                                           aGGATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTcc      <  1:544041/139‑1 (MQ=255)
                                                                                                            ggATTTTGACGGAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTccc     <  2:53274/139‑1 (MQ=255)
                                                                                                                       gAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGat  <  1:635820/131‑1 (MQ=255)
                                                                                                                       gAAGAGACGATTTTGTAAGACGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCAGCGAATATTTGCGAGCCGCTTTCCCGat  >  2:635820/1‑131 (MQ=255)
                                                                                                                                        |                                                                                                                  
TGTATTTGGGCGATGGGCAATTTATCGAGTCGCCACGTACCGGCGAAACCATTCGGATAAGCCGATTAGCCGAACCTTTCTGGCAGGACCATTTTTTGGGCGCGCGCAGGATTTTGACGGAAGAGACGATTTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACGAGCCCGCGAATATTTGCGAGCCGCTTTCCCGAT  >  NC_000913/247325‑247575

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: