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% 215.2 / NA 61F640S (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 (30/31);  total (30/31)

ACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAAC  >  NC_000913/631676‑631946
                                                                                                                                        |                                                                                                                                      
aCCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCca                                                                                                                                      >  2:1824790/1‑139 (MQ=255)
     aGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTAGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCaacaac                                                                                                                                 <  2:958722/139‑1 (MQ=255)
     aGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCAGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCaacaac                                                                                                                                 <  1:659630/139‑1 (MQ=255)
                            cGGATGCGAAAGTCGGCTTTCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  1:1374938/139‑1 (MQ=255)
                             ggATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAtc                                                                                                         >  1:938574/1‑139 (MQ=255)
                                           gCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTc                                                                                                                             >  1:1118914/1‑105 (MQ=255)
                                           gCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTc                                                                                                                             <  2:1118914/105‑1 (MQ=255)
                                           gCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          >  1:1467294/1‑124 (MQ=255)
                                           gCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  2:1467294/124‑1 (MQ=255)
                                               ccTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtg                                                                                       <  1:1824790/139‑1 (MQ=255)
                                                        tGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTg                                                                              >  1:20291/1‑139 (MQ=255)
                                                           taataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           <  2:647422/139‑1 (MQ=255)
                                                           taataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           <  2:20291/139‑1 (MQ=255)
                                                           taataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           <  2:1074083/139‑1 (MQ=255)
                                                             ataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCg                                                                                                                    <  2:376340/96‑1 (MQ=255)
                                                             ataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCg                                                                                                                    >  1:376340/1‑96 (MQ=255)
                                                                  ttCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAt                                                                                                 <  1:1369085/110‑1 (MQ=255)
                                                                  ttCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAt                                                                                                 >  2:1369085/1‑110 (MQ=255)
                                                                    ccAGTCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           <  2:1164035/130‑1 (MQ=255)
                                                                    ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           <  2:885181/130‑1 (MQ=255)
                                                                    ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           >  1:1164035/1‑130 (MQ=255)
                                                                    ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           >  1:885181/1‑130 (MQ=255)
                                                                    ccAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCg                                                                  >  1:1053111/1‑139 (MQ=255)
                                                                            tGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGtt                                                                      <  1:159565/127‑1 (MQ=255)
                                                                            tGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGtt                                                                      >  2:159565/1‑127 (MQ=255)
                                                                            tGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGAc                                                          >  1:569435/1‑139 (MQ=255)
                                                                                    aGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCAcgac                                                  <  1:254799/139‑4 (MQ=255)
                                                                                      cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTaa                                                                                                              <  1:1172650/77‑1 (MQ=255)
                                                                                      cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTaa                                                                                                              >  2:1172650/1‑77 (MQ=255)
                                                                                           atatTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtg                                                                                       >  1:734563/1‑95 (MQ=255)
                                                                                           atatTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtg                                                                                       <  2:734563/95‑1 (MQ=255)
                                                                                                 cGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGa                                     >  1:113417/1‑139 (MQ=255)
                                                                                                      cAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGa                                >  2:620065/1‑139 (MQ=255)
                                                                                                           taCTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCa                           <  1:620065/139‑1 (MQ=255)
                                                                                                                aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          >  1:523656/1‑55 (MQ=255)
                                                                                                                aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  2:523656/55‑1 (MQ=255)
                                                                                                                aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGc                      >  1:210385/1‑139 (MQ=255)
                                                                                                                aGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGc                      <  2:569435/139‑1 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGc                                                                                                                      >  1:1602263/1‑40 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGc                                                                                                                      <  2:1602263/40‑1 (MQ=38)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  2:1475232/52‑1 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          >  1:1475232/1‑52 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaa                   >  2:1129814/1‑139 (MQ=255)
                                                                                                                   cacaGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaa                   >  1:1185185/1‑139 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGcc                                                                                                                     <  1:179610/39‑1 (MQ=37)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGcc                                                                                                                     >  2:179610/1‑39 (MQ=37)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          >  2:1452672/1‑50 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  1:1452672/50‑1 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAt                                                                                                 >  2:500803/1‑59 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAt                                                                                                 <  1:500803/59‑1 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           <  1:1270436/81‑1 (MQ=255)
                                                                                                                     caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTACTTTATGGTGGTCGTGGTGGTTCTGGCa                                                                           >  2:1270436/1‑81 (MQ=255)
                                                                                                                      aGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaaga                <  2:375770/139‑1 (MQ=255)
                                                                                                                      aGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaaga                >  1:1743566/1‑139 (MQ=255)
                                                                                                                      aGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaaga                >  1:375770/1‑139 (MQ=255)
                                                                                                                        cTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCgaaagaaa              >  1:886484/1‑139 (MQ=255)
                                                                                                                          ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACg            <  2:886484/139‑1 (MQ=255)
                                                                                                                            cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGtt                                                                      <  2:140875/79‑1 (MQ=255)
                                                                                                                            cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGtt                                                                      >  1:140875/1‑79 (MQ=255)
                                                                                                                                    gtgtCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAAc  <  1:1129814/139‑1 (MQ=255)
                                                                                                                                    gtgtCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAAc  <  2:1185185/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAAC  >  NC_000913/631676‑631946

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