Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,513,468 T→C S218P (TCC→CCC)  ydcT → putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,513,4680TC100.0% 117.1 / NA 34S218P (TCC→CCC) ydcTputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/16);  total (18/16)

AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACA  >  NC_000913/1513337‑1513601
                                                                                                                                   |                                                                                                                                     
aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGGTCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                                <  1:2789/139‑1 (MQ=255)
  cTGGAACTGAAAAAACTGCCACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCTTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                              <  1:233880/139‑1 (MQ=255)
                                        ttttATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                        <  1:375720/139‑1 (MQ=255)
                                          ttATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                      >  2:745696/1‑139 (MQ=255)
                                                                      aGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                            <  2:85231/73‑1 (MQ=255)
                                                                      aGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                            >  1:85231/1‑73 (MQ=255)
                                                                       gCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGAc                                                         <  2:643782/139‑1 (MQ=255)
                                                                                 tGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCaga                                               >  1:635262/1‑139 (MQ=255)
                                                                                    ccGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                      >  2:100445/1‑97 (MQ=255)
                                                                                    ccGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGtt                                                                                      <  1:100445/97‑1 (MQ=255)
                                                                                       aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                              <  1:111321/54‑1 (MQ=255)
                                                                                       aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                              >  2:111321/1‑54 (MQ=255)
                                                                                       aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                                <  2:805024/52‑1 (MQ=255)
                                                                                       aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                                >  1:805024/1‑52 (MQ=255)
                                                                                           tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGGGCGATCTCTATATGCGCCCGCGCa                                                                                                         <  1:512759/71‑1 (MQ=255)
                                                                                           tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATTTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                             <  1:625185/115‑1 (MQ=255)
                                                                                           tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATTTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                             >  2:625185/1‑115 (MQ=255)
                                                                                           tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCa                                                                                                         >  2:512759/1‑71 (MQ=255)
                                                                                           tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTtgt                                     >  1:785211/1‑139 (MQ=255)
                                                                                                  gTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                              >  1:366212/1‑43 (MQ=255)
                                                                                                  gTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                              <  2:366212/43‑1 (MQ=255)
                                                                                                    tttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACg                            >  2:839602/1‑139 (MQ=255)
                                                                                                     ttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtct                                                                                                                          >  2:14010/1‑44 (MQ=255)
                                                                                                     ttCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAtct                                                                                                                          <  1:14010/44‑1 (MQ=255)
                                                                                                             tGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCgcgc                                                                                                          >  1:313609/1‑52 (MQ=255)
                                                                                                             tGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCgcgc                                                                                                          <  2:313609/52‑1 (MQ=255)
                                                                                                               gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                        >  2:286315/1‑68 (MQ=255)
                                                                                                               gACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                        <  1:286315/68‑1 (MQ=255)
                                                                                                                    aTTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGc            >  1:314265/1‑139 (MQ=255)
                                                                                                                            ggTCAATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGaa    >  2:576684/1‑139 (MQ=255)
                                                                                                                             gTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAAc   <  2:314265/139‑1 (MQ=255)
                                                                                                                              tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCTAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGAcc        <  2:762030/133‑1 (MQ=255)
                                                                                                                              tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGAcc        >  1:762030/1‑133 (MQ=255)
                                                                                                                              tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACa  >  2:780329/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACA  >  NC_000913/1513337‑1513601

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