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% 88.4 / NA 26S218P (TCC→CCC) ydcTputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/14);  total (12/14)

AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACA  >  NC_000913/1513337‑1513601
                                                                                                                                   |                                                                                                                                     
aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                                <  1:106806/139‑1 (MQ=255)
aGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcg                                                                                                                                <  2:36011/139‑1 (MQ=255)
  cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                              <  1:341759/139‑1 (MQ=255)
  cTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCcgcgcg                                                                                                                              >  2:169981/1‑139 (MQ=255)
      aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                            <  2:576003/137‑1 (MQ=255)
      aaCTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGAt                                                                                                                            >  1:576003/1‑137 (MQ=255)
                                        ttttATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCCATGTCCGATCGTGTGGCGGTTTTCAACAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGGCGTTTGTTGCCggg                                                                                        <  1:468461/139‑1 (MQ=255)
                                                                gggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATcccc                                                                                                                                    <  2:678539/71‑1 (MQ=255)
                                                                gggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATcccc                                                                                                                                    >  1:678539/1‑71 (MQ=255)
                                                                gggCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGttttt                                                                >  1:642818/1‑139 (MQ=255)
                                                                      aGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                             <  2:319285/136‑1 (MQ=255)
                                                                      aGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGAt                                                             >  1:319285/1‑136 (MQ=255)
                                                                               gATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTGTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCa                                                 <  1:44393/139‑1 (MQ=255)
                                                                                  gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                        >  2:313825/1‑97 (MQ=255)
                                                                                  gTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCggg                                                                                        <  1:313825/97‑1 (MQ=255)
                                                                                     cGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGaaa                                           >  2:134365/1‑139 (MQ=255)
                                                                                       aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCa                                                                                                         >  2:81802/1‑75 (MQ=255)
                                                                                       aTCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCa                                                                                                         <  1:81802/75‑1 (MQ=255)
                                                                                           tgtgGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTtgt                                     <  1:367938/139‑1 (MQ=255)
                                                                                                             tGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCtat                                                                                                                      <  1:587377/40‑1 (MQ=255)
                                                                                                             tGGACGCATTGAGCAGGTCGATCCCCCGCGCGATCTCtat                                                                                                                      >  2:587377/1‑40 (MQ=255)
                                                                                                                  gCATTGAGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCt              <  2:642818/139‑1 (MQ=255)
                                                                                                                        aGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGa                                                              >  1:288275/1‑85 (MQ=255)
                                                                                                                        aGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGa                                                              <  2:288275/85‑1 (MQ=255)
                                                                                                                        aGCAGGTCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGAcc        <  2:92860/139‑1 (MQ=255)
                                                                                                                              tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACa  >  2:581854/1‑139 (MQ=255)
                                                                                                                              tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACa  >  2:423872/1‑139 (MQ=255)
                                                                                                                              tCGATCCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACa  >  1:535263/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
AGCTGGAACTGAAAAAACTGCAACAGTCTCTCGGTATCACTTTTATCTTCGTTACCCACGATCAGGGCGAAGCGTTATCGATGTCCGATCGTGTGGCGGTTTTCAATAATGGACGCATTGAGCAGGTCGATTCCCCGCGCGATCTCTATATGCGCCCGCGCACGCCGTTTGTTGCCGGGTTCGTTGGTACATCGAATGTTTTTGATGGACTGATGGCAGAGAAACTTTGTGGCATGACGGGAAGCTTCGCCCTGCGACCGGAACA  >  NC_000913/1513337‑1513601

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