Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,950,998 T→C V56A (GTT→GCT)  yecD → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,9980TC100.0% 187.6 / NA 53V56A (GTT→GCT) yecDputative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/26);  total (27/26)

GATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAA  >  NC_000913/1950870‑1951134
                                                                                                                                |                                                                                                                                        
gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTg                                                                                                                               >  2:349249/1‑140 (MQ=255)
    gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGa                                                                                                                           <  1:785410/140‑1 (MQ=255)
    gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGa                                                                                                                           <  2:1227065/140‑1 (MQ=255)
     aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGAt                                                                                                                          <  1:603394/140‑1 (MQ=255)
     aTTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGAt                                                                                                                          <  2:188458/140‑1 (MQ=255)
            agaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTg                                                                                                                   <  2:326869/140‑1 (MQ=255)
             gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGTCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGa                                                                                                                  <  1:924277/140‑1 (MQ=255)
             gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTTTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGa                                                                                                                  >  1:1010190/1‑140 (MQ=255)
             gaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGa                                                                                                                  >  1:226319/1‑140 (MQ=255)
                ggCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGc                                                                                                               >  2:30372/1‑140 (MQ=255)
                          cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTTTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAg                                                                                                     <  2:1010190/140‑1 (MQ=255)
                          cTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAg                                                                                                     <  2:226319/140‑1 (MQ=255)
                               gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGt                                                                                                >  2:530377/1‑140 (MQ=255)
                               gCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGt                                                                                                >  1:867555/1‑140 (MQ=255)
                                    tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                                                                                           >  1:1213946/1‑140 (MQ=255)
                                        ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCC‑‑GTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCt                                                                                     >  2:525123/1‑140 (MQ=255)
                                             tACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTATGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAc                                                                                  <  1:972527/140‑1 (MQ=255)
                                             tACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAc                                                                                  <  1:832907/140‑1 (MQ=255)
                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   >  2:796537/1‑136 (MQ=255)
                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   >  2:21201/1‑136 (MQ=255)
                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   <  1:21201/136‑1 (MQ=255)
                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   >  2:608673/1‑136 (MQ=255)
                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   <  1:608673/136‑1 (MQ=255)
                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   <  1:796537/136‑1 (MQ=255)
                                                tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCgg                                                                               >  2:1130494/1‑140 (MQ=255)
                                                tGCCGATGAGGGAGTTAATCGCGCCGGGGAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCgg                                                                               >  2:1089174/1‑140 (MQ=255)
                                                 gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   >  1:33123/1‑135 (MQ=255)
                                                 gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa                                                                                   <  2:33123/135‑1 (MQ=255)
                                                 gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc                                                                              >  1:848683/1‑140 (MQ=255)
                                                 gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc                                                                              >  2:599026/1‑140 (MQ=255)
                                                  ccGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCt                                                                             >  2:957161/1‑140 (MQ=255)
                                                   cGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTa                                                                            >  1:26706/1‑140 (MQ=255)
                                                   cGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTa                                                                            <  1:30372/140‑1 (MQ=255)
                                                           tAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCC‑‑GTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGcc                                                                  <  1:525123/140‑1 (MQ=255)
                                                            aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGc                                                                   >  1:1142458/1‑140 (MQ=255)
                                                                   cgcgCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGaaaa                                                            >  1:266503/1‑140 (MQ=255)
                                                                    gcgcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAAt                                                           >  2:302520/1‑140 (MQ=255)
                                                                     cgcCGGGAAGCTGGCGGTGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAAtt                                                          >  1:241525/1‑140 (MQ=255)
                                                                       ccGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtgg                                                        >  2:168739/1‑140 (MQ=255)
                                                                       ccGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtgg                                                        >  1:463785/1‑140 (MQ=255)
                                                                               cTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACa                                                <  1:349249/140‑1 (MQ=255)
                                                                                tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACAt                                               <  1:1130494/140‑1 (MQ=255)
                                                                                 ggcggcGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATc                                              <  2:26706/140‑1 (MQ=255)
                                                                                    ggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctg                                           <  1:599026/140‑1 (MQ=255)
                                                                                    ggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctg                                           <  2:867555/140‑1 (MQ=255)
                                                                                      cgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCctgct                                         <  1:302520/140‑1 (MQ=255)
                                                                                        aaaTTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtgg                                                        <  1:566982/123‑1 (MQ=255)
                                                                                        aaaTTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATtgg                                                        >  2:566982/1‑123 (MQ=255)
                                                                                            ttCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTa                                   >  1:1259593/1‑140 (MQ=255)
                                                                                             tCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAg                                  <  2:1142458/140‑1 (MQ=255)
                                                                                                  ccAGCGGTCAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCa                             <  2:848683/140‑1 (MQ=255)
                                                                                                          cAGCCCGTGTTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAg                     <  1:530377/140‑1 (MQ=255)
                                                                                                                gtgtTTCTGGTGCGCGCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGatat               >  1:1057952/1‑140 (MQ=255)
                                                                                                                           gcgcgcTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAAtcatca    <  2:463785/140‑1 (MQ=255)
                                                                                                                           gcgcgcTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAAtcatca    >  2:1111544/1‑140 (MQ=255)
                                                                                                                            cgcgCTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCaa   >  2:287722/1‑140 (MQ=255)
                                                                                                                             gcgcTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCaaa  <  2:1213946/140‑1 (MQ=255)
                                                                                                                                |                                                                                                                                        
GATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAA  >  NC_000913/1950870‑1951134

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