Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	118,912	A→G	A60A (GCA→GCG)	ampD →	1,6‑anhydro‑N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	118,912	0	A	G	100.0%	170.6 / NA	48	A60A (GCA→GCG)	ampD	1,6‑anhydro‑N‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/27);  total (21/27)

ACATTACGATTGCCGCCCGGATGACGAAACACCCACCCTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCACATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCG  >  NC_000913/118783‑119048                                                                                                                                  |                                                                                                                                         acaTTACGATTGCCGCCCGGATGACGAAACACCCACCCTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCctttct                                                                                                                                <  2:642439/140‑1 (MQ=255)      aCGATTGCCGCCCGGATGACGAAACACCCACCCTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCt                                                                                                                           <  1:684204/140‑1 (MQ=255)             ccgccCGGATGACGAAACACCCACCCTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCg                                                                                                                    >  1:361095/1‑140 (MQ=255)                                 ccacccTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCt                                                                                                <  1:225474/140‑1 (MQ=255)                                 ccacccTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCt                                                                                                <  1:463299/140‑1 (MQ=255)                                  cacccTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTc                                                                                               <  1:441635/140‑1 (MQ=255)                                    cccTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCAc                                                                                             >  2:456684/1‑140 (MQ=255)                                            gtggtgCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGAt                                                                                     <  1:456684/140‑1 (MQ=255)                                               gtgCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCg                                                                                  >  2:160261/1‑140 (MQ=255)                                                  cacaATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCg                                                                               >  2:200940/1‑140 (MQ=255)                                                   acaATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGt                                                                              >  1:244322/1‑140 (MQ=255)                                                   acaATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGt                                                                              <  2:361095/140‑1 (MQ=255)                                                   acaATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGt                                                                              <  1:379183/140‑1 (MQ=255)                                                     aaTATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGa                                                                            <  1:200940/140‑1 (MQ=255)                                                     aaTATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGa                                                                            <  2:664331/140‑1 (MQ=255)                                                       tatTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGt                                                                              <  1:398962/136‑1 (MQ=255)                                                       tatTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGt                                                                              >  2:398962/1‑136 (MQ=255)                                                          tAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTg                                                                       >  1:211149/1‑140 (MQ=255)                                                            gCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGaa                                                                     >  1:73207/1‑140 (MQ=255)                                                             ccTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGt                                                                              <  1:686666/130‑1 (MQ=255)                                                             ccTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGt                                                                              >  2:686666/1‑130 (MQ=255)                                                             ccTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGaaa                                                                    <  1:635656/140‑1 (MQ=255)                                                             ccTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGaaa                                                                    <  1:322905/140‑1 (MQ=255)                                                              cTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAAt                                                                   >  2:727573/1‑140 (MQ=255)                                                                 ccgccAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGt                                                                >  2:475674/1‑140 (MQ=255)                                                                 ccgccAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGt                                                                >  2:397398/1‑140 (MQ=255)                                                                 ccgccAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGt                                                                >  2:312116/1‑140 (MQ=255)                                                                 ccgccAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTCCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGt                                                                >  1:711144/1‑140 (MQ=255)                                                                 ccgccAGGCGAGTTTGGCGGTCCGTGGATCGACGCATCATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCAGTGATgg                                                                         >  2:317979/1‑131 (MQ=255)                                                                 ccgccAGGCGAGTTTGGCGGTCCGTGGATCGACGCATCATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCAGTGATgg                                                                         <  1:317979/131‑1 (MQ=255)                                                                    ccAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCa                                                             <  1:805751/140‑1 (MQ=255)                                                                                      ccGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATaa                                           >  2:755442/1‑140 (MQ=255)                                                                                         tGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACg                                        <  1:727573/140‑1 (MQ=255)                                                                                                       attCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCgg                          <  1:755442/140‑1 (MQ=255)                                                                                                        ttCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCggg                         <  1:475674/140‑1 (MQ=255)                                                                                                         tCACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGa                        <  2:244322/140‑1 (MQ=255)                                                                                                          cACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAg                       <  2:83363/140‑1 (MQ=255)                                                                                                          cACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAg                       <  2:73207/140‑1 (MQ=255)                                                                                                          cACTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAg                       <  2:491908/140‑1 (MQ=255)                                                                                                           aCTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGt                      <  2:472749/140‑1 (MQ=255)                                                                                                           aCTGGAACTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGt                      <  2:31679/140‑1 (MQ=255)                                                                                                                aaCTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGtctctc                 >  1:474364/1‑140 (MQ=255)                                                                                                                  cTATTGATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAg               >  1:42269/1‑140 (MQ=255)                                                                                                                       gATCCGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCa          >  2:529771/1‑140 (MQ=255)                                                                                                                          ccGCAGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAggg       >  2:252632/1‑140 (MQ=255)                                                                                                                              aGGCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGgcgc   <  1:160261/140‑1 (MQ=255)                                                                                                                               ggCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGcgcg  >  2:609687/1‑140 (MQ=255)                                                                                                                               ggCGCATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGcgcg  <  2:211149/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         ACATTACGATTGCCGCCCGGATGACGAAACACCCACCCTGCTGGTGGTGCACAATATTAGCCTGCCGCCAGGCGAGTTTGGCGGTCCGTGGATCGACGCATTATTCACTGGAACTATTGATCCGCAGGCACATCCTTTCTTTGCTGAGATCGCCCATTTGCGCGTCTCCGCTCACTGTTTGATTCGCCGTGATGGTGAAATAGTCCAGTATGTTCCTTTCGATAAACGTGCATGGCATGCGGGAGTCTCTCAGTATCAGGGGCGCG  >  NC_000913/118783‑119048

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: ‑