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% 232.3 / NA 64F640S (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 (32/32);  total (32/32)

GCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAAC  >  NC_000913/631679‑631946
                                                                                                                                     |                                                                                                                                      
gCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCcaaca                                                                                                                                   <  1:154023/139‑1 (MQ=255)
      tGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTc                                                                                                                             <  2:162400/139‑1 (MQ=255)
              aGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGcc                                                                                                                     >  1:7831/1‑139 (MQ=255)
              aGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGcc                                                                                                                     >  1:122335/1‑139 (MQ=255)
                  tttAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                                 >  1:537643/1‑139 (MQ=255)
                  tttAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCGAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                                 <  2:691176/139‑1 (MQ=255)
                         cGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  1:486870/139‑1 (MQ=255)
                                 aaaGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACACCCGTCTGGATGCCGGGTTAACCATCTTCTTTa                                                                                                  >  2:374878/1‑139 (MQ=255)
                                     tCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGc                                                                                                                      <  2:995402/115‑1 (MQ=255)
                                     tCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGc                                                                                                                      >  1:995402/1‑115 (MQ=255)
                                       ggCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtggtgg                                                                                            >  1:598943/1‑139 (MQ=255)
                                                 ccATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTggtggt                                                                                  >  2:502349/1‑139 (MQ=255)
                                                          ataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCaacaac                                                                                                                                 <  2:805514/83‑1 (MQ=255)
                                                          ataaGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCaacaac                                                                                                                                 >  1:805514/1‑83 (MQ=255)
                                                              gTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTtc                                                                     <  2:7831/139‑1 (MQ=255)
                                                               ttCCAGGCAATGATCGACAACGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCt                                                                                                                            <  2:410620/83‑1 (MQ=255)
                                                               ttCCAGGCAATGATCGACAACGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCt                                                                                                                            >  1:410620/1‑83 (MQ=255)
                                                                    ggCAATGATCGACAGTGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACTACCGTCt                                                                                                                            <  2:937333/78‑1 (MQ=255)
                                                                    ggCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCt                                                                                                                            >  1:937333/1‑78 (MQ=255)
                                                                    ggCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  2:820366/96‑1 (MQ=255)
                                                                    ggCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          >  1:820366/1‑96 (MQ=255)
                                                                    ggCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGAtt                                                               >  2:1059893/1‑139 (MQ=255)
                                                                      cAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTaa                                                             <  1:1059893/139‑1 (MQ=255)
                                                                                   cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCg                                                <  1:783847/139‑1 (MQ=255)
                                                                                   cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCg                                                >  2:783847/1‑139 (MQ=255)
                                                                                   cGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGAGGCCGGGTTAACCATCTTCTTTATGGTGGTCGGGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCg                                                <  2:398681/139‑1 (MQ=255)
                                                                                    ggCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATg                                                                                                                       >  1:425830/1‑67 (MQ=255)
                                                                                    ggCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATg                                                                                                                       <  2:425830/67‑1 (MQ=255)
                                                                                     gCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGc                                              >  1:520062/1‑139 (MQ=255)
                                                                                      cAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                                  <  2:710245/134‑1 (MQ=255)
                                                                                      cAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                                  >  1:710245/1‑134 (MQ=255)
                                                                                            tCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGaaa                                       <  2:520062/139‑1 (MQ=255)
                                                                                               gTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtgg                                                                                               <  1:1309782/80‑1 (MQ=255)
                                                                                               gTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtgg                                                                                               >  2:1309782/1‑80 (MQ=255)
                                                                                                tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAAcc                                                                                                                                >  1:1522/1‑46 (MQ=255)
                                                                                                tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAAcc                                                                                                                                <  2:1522/46‑1 (MQ=255)
                                                                                                tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                                 <  2:994570/61‑1 (MQ=255)
                                                                                                tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGt                                                                                                                 >  1:994570/1‑61 (MQ=255)
                                                                                                tCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtggt                                                                                     >  1:905527/1‑89 (MQ=255)
                                                                                                tCGCAGTATACAGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtggt                                                                                     <  2:905527/89‑1 (MQ=255)
                                                                                                  gCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCg                                 <  2:537643/139‑1 (MQ=255)
                                                                                                     gTATACTGAGTCACAGCTGGTGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  1:673335/63‑1 (MQ=255)
                                                                                                     gTATACTGAGTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          >  2:673335/1‑63 (MQ=255)
                                                                                                             tgTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtg                                                                                       <  2:784814/73‑1 (MQ=255)
                                                                                                             tgTCACAGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtg                                                                                       >  1:784814/2‑74 (MQ=255)
                                                                                                              gTCACAGCTAGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtggt                                                                                     <  1:435343/75‑1 (MQ=255)
                                                                                                              gTCACAGCTAGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCgtggt                                                                                     >  2:435343/1‑75 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          >  2:709825/1‑50 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCAt                                                                                                          <  1:709825/50‑1 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTggtggt                                                                                  <  2:218082/74‑1 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTggtggt                                                                                  >  1:218082/1‑74 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                                  <  1:1053298/106‑1 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                                  >  2:1053298/1‑106 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCg                                                <  2:153458/108‑1 (MQ=255)
                                                                                                                  caGCTGGCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCg                                                >  1:153458/1‑108 (MQ=255)
                                                                                                                       ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAt                                                                                                 <  1:765319/54‑1 (MQ=255)
                                                                                                                       ggCGCAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAt                                                                                                 >  2:765319/1‑54 (MQ=255)
                                                                                                                        gcgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtgg                                                                                               >  1:300074/1‑55 (MQ=255)
                                                                                                                        gcgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtgg                                                                                               <  2:300074/55‑1 (MQ=255)
                                                                                                                         cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGAt                                                                <  1:1150134/85‑1 (MQ=255)
                                                                                                                         cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGAt                                                                >  2:1150134/1‑85 (MQ=255)
                                                                                                                         cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                                  <  1:404966/99‑1 (MQ=255)
                                                                                                                         cgcAACTGGTGTCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTg                                                  >  2:404966/1‑99 (MQ=255)
                                                                                                                                 gtgtCCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAAc  >  2:262657/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                      
GCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAAC  >  NC_000913/631679‑631946

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