Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,812 T→C F640S (TTC→TCC)  cstA → carbon starvation protein involved in peptide utilization, APC peptide transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,8120TC100.0% 113.5 / NA 32F640S (TTC→TCC) cstAcarbon starvation protein involved in peptide utilization, APC peptide transporter family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

ACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACG  >  NC_000913/631676‑631936
                                                                                                                                        |                                                                                                                            
aCCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCca                                                                                                                            >  1:431025/1‑139 (MQ=255)
  cGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCcaac                                                                                                                          <  1:899203/139‑1 (MQ=255)
               aaaGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATg                                                                                                             <  2:334187/139‑1 (MQ=255)
                            cGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACGACCGTCTGGATGCCGGGTTAACCAt                                                                                                <  2:575314/139‑1 (MQ=255)
                                          ggCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtggtgg                                                                                  >  2:356987/1‑139 (MQ=255)
                                                cTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtgg                                                                            <  1:411362/139‑1 (MQ=255)
                                                           taataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                 <  1:356987/139‑1 (MQ=255)
                                                           taataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                 <  1:758167/139‑1 (MQ=255)
                                                                       ggCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGAtt                                                     >  2:105936/1‑139 (MQ=255)
                                                                                     gCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTTGATGCCGGGTTAACCATCTTCTTTATGGTGGGCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGc                                       <  2:458393/139‑1 (MQ=255)
                                                                                      cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCg                                      <  2:380004/139‑1 (MQ=255)
                                                                                          aaTATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtggtgg                                                                                  >  1:25208/1‑91 (MQ=255)
                                                                                          aaTATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtggtgg                                                                                  <  2:25208/91‑1 (MQ=255)
                                                                                                      cAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGa                      >  1:853360/1‑139 (MQ=255)
                                                                                                              tGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACg              >  1:242359/1‑139 (MQ=255)
                                                                                                                 gTCACAGCTGGCGCAACTGGTGTCCAACATCCGTCTGGATGCCGGGTTAACCAt                                                                                                <  2:631161/54‑1 (MQ=38)
                                                                                                                 gTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTCGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCg           >  2:550635/1‑139 (MQ=255)
                                                                                                                 gTCACAGCTGGCGCAACTGGGGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                >  1:631161/1‑54 (MQ=39)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                <  2:710534/52‑1 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                >  1:710534/1‑52 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtg                                                                                      <  1:6585/62‑1 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtg                                                                                      >  2:6585/1‑62 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATggtggt                                                                                 >  1:445240/1‑67 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATggtggt                                                                                 <  2:445240/67‑1 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaa         >  2:354268/1‑139 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                >  1:455327/1‑50 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                <  2:455327/50‑1 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATggtggt                                                                                 >  2:69901/1‑65 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATggtggt                                                                                 <  1:69901/65‑1 (MQ=255)
                                                                                                                          ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                        <  2:407410/101‑1 (MQ=255)
                                                                                                                          ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                        >  1:407410/1‑101 (MQ=255)
                                                                                                                          ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTATTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACg  <  1:670390/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                            
ACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACG  >  NC_000913/631676‑631936

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