Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,508,711	T→G	L644R (CTG→CGG)	ydcP →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,508,711	0	T	G	100.0%	167.5 / NA	47	L644R (CTG→CGG)	ydcP	putative peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/26);  total (21/26)

AATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTT  >  NC_000913/1508565‑1508857                                                                                                                                                   |                                                                                                                                                   aaTGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAgcggc                                                                                                                                                  <  1:738137/149‑1 (MQ=255) aaTGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGcggctg                                                                                                                                                <  1:645865/151‑1 (MQ=255)    gCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGaa                                                                                                                                              <  2:700789/150‑1 (MQ=255)             tAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCg                                                                                                                                    >  1:36466/1‑151 (MQ=255)              aaCGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCg                                                                                                                                    <  2:840007/150‑1 (MQ=255)                 ggCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGc                                                                                                                                <  2:684834/151‑1 (MQ=255)                     aTGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGaaa                                                                                                                            <  1:424439/151‑1 (MQ=255)                       gAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGaaaa                                                                                                                           <  2:292329/150‑1 (MQ=255)                       gAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAg                                                                                                                          <  2:524662/151‑1 (MQ=255)                         aGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAggg                                                                                                                        <  2:782580/151‑1 (MQ=255)                            aTTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTaa                                                                                                                     <  1:272123/151‑1 (MQ=255)                            aTTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTa                                                                                                                      <  2:22294/150‑1 (MQ=255)                             ttAACCCTAACGGTTGATTGCCGCCAATGCGAGATGCACGTGATTGGCATTATCAAAAATCACATACTGAGAATGCCGTTACCGGGAAGCGTAGTGTCATCCGGAAGGCCGGATGAGCGGCTGAAAACATTATCGTTGCGAAAAgag                                                                                                                        >  1:849637/1‑145 (MQ=255)                              tAACGATAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTaa                                                                                                                     <  1:446945/149‑1 (MQ=255)                                  gCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGcc                                                                                                               >  1:505059/1‑151 (MQ=255)                                       aGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGttt                                                                                                          >  1:898379/1‑151 (MQ=255)                                       aGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGttt                                                                                                          <  1:94713/151‑1 (MQ=255)                                        gTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGtttt                                                                                                         <  2:788684/151‑1 (MQ=255)                                          ttGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCt                                                                                                       >  1:173309/1‑151 (MQ=255)                                           tGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCt                                                                                                       <  2:237008/150‑1 (MQ=255)                                            gATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTgg                                                                                                     <  2:830300/151‑1 (MQ=255)                                                    cccATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTcac                                                                                             >  1:88897/1‑151 (MQ=255)                                                       aTGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAAc                                                                                          >  1:348436/1‑151 (MQ=255)                                                         gCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACtt                                                                                        <  1:636858/151‑1 (MQ=255)                                                           gagaTGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTAt                                                                                      <  1:323526/151‑1 (MQ=255)                                                               tGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGa                                                                                   >  1:68650/1‑150 (MQ=255)                                                                gCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAAt                                                                                 <  1:215550/151‑1 (MQ=255)                                                                          ggCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGat                                                                        >  1:553487/1‑150 (MQ=255)                                                                           gCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGatat                                                                      >  2:463182/1‑151 (MQ=255)                                                                             aaaatcaaaaatcaCATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATtg                                                                    <  1:412184/151‑1 (MQ=255)                                                                                       atcaCATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGt                                                          >  2:811640/1‑151 (MQ=255)                                                                                         cacaTACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTaa                                                        <  1:517494/151‑1 (MQ=255)                                                                                          acaTACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTa                                                         >  2:26866/1‑149 (MQ=255)                                                                                          acaTACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAg                                                       <  2:883039/151‑1 (MQ=255)                                                                                            aTACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGc                                                     <  1:75289/151‑1 (MQ=255)                                                                                                 gAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGAGGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGGAAGGCGTTAc                                                >  1:15801/1‑151 (MQ=255)                                                                                                   aaaTGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAcgc                                              >  1:810649/1‑151 (MQ=255)                                                                                                   aaaTGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAcgc                                              >  1:347832/1‑151 (MQ=255)                                                                                                     aTGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTAcgccg                                            >  2:765581/1‑151 (MQ=255)                                                                                                        ccGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGCTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCa                                          >  2:637591/1‑150 (MQ=255)                                                                                                            tACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCgg                                     >  1:264211/1‑151 (MQ=255)                                                                                                             aCCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGc                                    >  2:164226/1‑151 (MQ=255)                                                                                                              ccGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGc                                    >  1:845500/1‑150 (MQ=255)                                                                                                                            ggCATCCGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAAc                     <  1:277828/151‑1 (MQ=255)                                                                                                                                 ccGTAAGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGt                <  2:360131/151‑1 (MQ=255)                                                                                                                                     aaGTCCGGATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCg            <  2:689643/151‑1 (MQ=255)                                                                                                                                           ggATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCGGGGTAAGGCGATACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTg      >  1:640639/1‑151 (MQ=255)                                                                                                                                           ggATGAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTg      <  1:539923/151‑1 (MQ=255)                                                                                                                                               gAGCGGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGttt  <  2:578780/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCGCCCATGCGAGATGCACGTCATTGGCAAAATCAAAAATCACATACTGAAAATGCCGTTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTT  >  NC_000913/1508565‑1508857

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