Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,179 A→G N429S (AAC→AGC)  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,1790AG100.0% 175.6 / NA 46N429S (AAC→AGC) cstAcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/26);  total (20/26)

GTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAACCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGT  >  NC_000913/631031‑631325
                                                                                                                                                    |                                                                                                                                                  
gTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGcc                                                                                                                                                  <  2:1023385/151‑1 (MQ=255)
   tCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGcccg                                                                                                                                                >  2:975044/1‑148 (MQ=255)
     cGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGGGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGAtt                                                                                                                                             >  2:465323/1‑151 (MQ=255)
          gTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCa                                                                                                                                         <  2:1044742/150‑1 (MQ=255)
            atcatcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGc                                                                                                                                       <  1:647424/150‑1 (MQ=255)
              catcGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGaa                                                                                                                                    <  2:846906/151‑1 (MQ=255)
                 cGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTg                                                                                                                                 <  1:940249/151‑1 (MQ=255)
                     ccGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTg                                                                                                                             >  1:424030/1‑151 (MQ=255)
                     ccGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTg                                                                                                                             >  2:47956/1‑151 (MQ=255)
                        ggCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAAc                                                                                                                          >  2:354242/1‑151 (MQ=255)
                             gtATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTc                                                                                                                      <  1:459258/150‑1 (MQ=255)
                                  ttGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCa                                                                                                                 <  1:536409/150‑1 (MQ=255)
                                     ccATGAACAGCCCGATGGCGGGGCTGGCTCCGGCAGGGACGGCGGCTGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTc                                                                                                             >  2:634346/1‑151 (MQ=255)
                                       aTGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTccc                                                                                                           <  2:240575/151‑1 (MQ=255)
                                       aTGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTccc                                                                                                           >  2:770423/1‑151 (MQ=255)
                                          aaCAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTg                                                                                                        <  2:94098/151‑1 (MQ=255)
                                                     ggCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTgcg                                                                                             <  2:450897/151‑1 (MQ=255)
                                                       cGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGcc                                                                                           <  1:228601/151‑1 (MQ=255)
                                                             tGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGc                                                                                      >  2:281503/1‑150 (MQ=255)
                                                               gCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTgg                                                                                   <  1:149618/151‑1 (MQ=255)
                                                                    gtcAGGGACGGCGGATGGGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTg                                                                              >  2:50351/3‑151 (MQ=255)
                                                                     gCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTgg                                                                             >  1:26423/1‑151 (MQ=255)
                                                                             ggcggATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGcc                                                                     >  1:991174/1‑151 (MQ=255)
                                                                                ggATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTAc                                                                  <  1:371019/151‑1 (MQ=255)
                                                                                           cTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACgg                                                       >  1:28531/1‑151 (MQ=255)
                                                                                             tCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACgg                                                       <  1:93156/149‑1 (MQ=255)
                                                                                                gCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGgcgcg                                                   <  1:546828/150‑1 (MQ=255)
                                                                                                gCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGgcgc                                                    <  1:570387/149‑1 (MQ=255)
                                                                                                            gtgAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGGGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCat                                        >  2:119865/1‑149 (MQ=255)
                                                                                                                  aGCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATgatggatg                                >  2:129482/1‑151 (MQ=255)
                                                                                                                   gCTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGAtgt                               >  1:507354/1‑151 (MQ=255)
                                                                                                                    cTGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGAtgtg                              >  1:634698/1‑151 (MQ=255)
                                                                                                                     tGGGGCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCCACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTgg                             >  2:522726/1‑151 (MQ=255)
                                                                                                                       gggCTTTAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCg                           <  1:761826/151‑1 (MQ=255)
                                                                                                                             tAGCATTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTg                      <  1:929401/150‑1 (MQ=255)
                                                                                                                                 aTTACACCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATc                 <  2:371838/151‑1 (MQ=255)
                                                                                                                                    acacCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCAc               >  2:134887/1‑150 (MQ=255)
                                                                                                                                    acacCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACt              <  1:4614/151‑1 (MQ=255)
                                                                                                                                    acacCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACt              <  1:1027757/151‑1 (MQ=255)
                                                                                                                                     cacCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACtt             <  2:992386/151‑1 (MQ=255)
                                                                                                                                      acCAGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGGATCACtt             >  2:488163/1‑150 (MQ=255)
                                                                                                                                         aGATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATTGATGTGGCGTTCTGGTATCACTTCGcc         <  1:212023/151‑1 (MQ=255)
                                                                                                                                          gATACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCa        <  2:242754/151‑1 (MQ=255)
                                                                                                                                           aTACGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCAt       <  2:927804/151‑1 (MQ=255)
                                                                                                                                              cGCTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGCCGCTGGCGTTGGGGATGGCCTACATTCTGCCCGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCAttt     <  2:534789/150‑1 (MQ=255)
                                                                                                                                                cTAAGCCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt  >  1:93410/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGGGACGGCGGATGTGGTCGCTTCTGCCGCGCAGGTGGTGAGTAGCTGGGGCTTTAGCATTACACCAGATACGCTAAACCAGATTGCCAGCGAAGTGGGTGAACAGTCGATCATTTCCCGTGCGGGCGGTGCGCCGACGCTGGCGGTGGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGT  >  NC_000913/631031‑631325

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