Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,266,510	C→T	A243A (GCC→GCT)	alr →	alanine racemase, biosynthetic, PLP‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,266,510	0	C	T	100.0%	185.8 / NA	53	A243A (GCC→GCT)	alr	alanine racemase, biosynthetic, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/24);  total (29/24)

GTGGCCACAGTCGCATTTTGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCCGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAG  >  NC_000913/4266378‑4266643                                                                                                                                     |                                                                                                                                      gtgGCCACAGTCGCATTTTGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCtgtgcgt                                                                                                                                 <  1:697361/139‑1 (MQ=255)    gCCACAGTCGCATTTTGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAg                                                                                                                              >  1:1096002/1‑139 (MQ=255)         aGTCGCATTTTGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATaa                                                                                                                         >  2:209412/1‑139 (MQ=255)         aGTCGCATTTTGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATaa                                                                                                                         >  1:436500/1‑139 (MQ=255)              cATTTTGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCg                                                                                                                    >  1:390778/1‑139 (MQ=255)                 tttGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGgag                                                                                                                 <  2:1096002/139‑1 (MQ=255)                      cTGGGTGCGCCCGGGCATCA‑TCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTg                                                                                                           <  1:485579/139‑1 (MQ=255)                         ggTGCGCCCGGGCATCATTCTTTATGGTGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATa                                                                                                                          <  1:1194311/122‑1 (MQ=255)                         ggTGCGCCCGGGCATCATTCTTTATGGTGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATa                                                                                                                          >  2:1194311/1‑122 (MQ=255)                          gTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTg                                                                                                        >  2:1061502/1‑139 (MQ=255)                               cccGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTAt                                                                                                   >  2:254426/1‑139 (MQ=255)                                 cGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTAtgg                                                                                                 <  2:436500/139‑1 (MQ=255)                                  gggCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTAtggt                                                                                                <  2:416863/139‑1 (MQ=255)                                                  ggCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCtgtgcgt                                                                                                                                 >  1:1084220/1‑90 (MQ=255)                                                  ggCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCtgtgcgt                                                                                                                                 <  2:1084220/90‑1 (MQ=255)                                                    cGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGaa                                                                                            >  1:288625/1‑125 (MQ=255)                                                    cGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGaa                                                                                            <  2:288625/125‑1 (MQ=255)                                                    cGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGaa                                                                              >  2:1367016/1‑139 (MQ=255)                                                     gTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGga                                                                                                                  <  1:356759/102‑1 (MQ=255)                                                     gTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGga                                                                                                                  >  2:356759/1‑102 (MQ=255)                                                      tctcGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGt                                                                                     <  2:784082/130‑1 (MQ=255)                                                      tctcGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGt                                                                                     >  1:784082/1‑130 (MQ=255)                                                      tctcGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACg                                                                            <  1:1377149/139‑1 (MQ=255)                                                      actcGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCCGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGTACCTGGGTAAGCGAACg                                                                            >  1:928435/2‑139 (MQ=255)                                                       ctcGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGt                                                                           >  1:266908/1‑139 (MQ=255)                                                       ctcGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGt                                                                           >  2:1281620/1‑139 (MQ=255)                                                         cgccgcTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGa                                                                         >  2:1171324/1‑139 (MQ=255)                                                                  aaGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCg                                                                                <  2:947679/123‑1 (MQ=255)                                                                  aaGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCg                                                                                >  1:947679/1‑123 (MQ=255)                                                                         cTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTa                                                         >  2:247430/1‑139 (MQ=255)                                                                           ccACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGt                                                       <  1:1367016/139‑1 (MQ=255)                                                                               cGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTcgcg                                                   <  1:555385/139‑1 (MQ=255)                                                                               cGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTcgcg                                                   <  1:1367384/139‑1 (MQ=255)                                                                                ggTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGa                                                  <  1:209412/139‑1 (MQ=255)                                                                                   gCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGga                                                                                                                  <  2:438578/72‑1 (MQ=255)                                                                                   gCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGga                                                                                                                  >  1:438578/1‑72 (MQ=255)                                                                                    ccGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATAc                                                                      <  1:207663/115‑1 (MQ=255)                                                                                    ccGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATAc                                                                      >  2:207663/1‑115 (MQ=255)                                                                                    ccGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATggg                                              >  2:739297/1‑139 (MQ=255)                                                                                      gATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCt                                            <  1:170719/139‑1 (MQ=255)                                                                                      gATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCt                                            <  2:266908/139‑1 (MQ=255)                                                                                           tGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGc                                       <  2:506819/139‑1 (MQ=255)                                                                                                           aTGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGga                                                                                                                  <  1:1426391/48‑1 (MQ=255)                                                                                                           aTGTCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGga                                                                                                                  >  2:1426391/1‑48 (MQ=255)                                                                                                              tCACTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGcgcgccg                    >  2:658662/1‑139 (MQ=255)                                                                                                                aCTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCTAACGTGATACCCGTCTTGGCGTAGTCGCGATg                                                <  2:525984/109‑1 (MQ=255)                                                                                                                aCTAACCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATg                                                >  1:525984/1‑109 (MQ=255)                                                                                                                   aaCCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGcgccgcgccg               >  1:1033329/1‑139 (MQ=255)                                                                                                                   aaCCTCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGcgccgcgccg               >  2:86785/1‑139 (MQ=255)                                                                                                                       tCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTAtggt                                                                                                <  2:218677/54‑1 (MQ=255)                                                                                                                       tCCAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTAtggt                                                                                                >  1:218677/1‑54 (MQ=255)                                                                                                                        ccAGCCTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCgg          <  1:1281620/139‑1 (MQ=255)                                                                                                                            ccTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGtt                                                                                                     >  2:1422557/1‑44 (MQ=255)                                                                                                                            ccTGATTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGtt                                                                                                     <  1:1422557/44‑1 (MQ=255)                                                                                                                                aTTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAg  >  2:1168180/1‑139 (MQ=255)                                                                                                                                aTTGCTGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAg  >  1:303155/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                      GTGGCCACAGTCGCATTTTGACTGGGTGCGCCCGGGCATCATTCTTTATGGCGTCTCGCCGCTGGAAGATCGCTCCACCGGTGCCGATTTTGGCTGTCAGCCAGTGATGTCACTAACCTCCAGCCTGATTGCCGTGCGTGAGCATAAAGCCGGAGAGCCTGTTGGTTATGGTGGAACCTGGGTAAGCGAACGTGATACCCGTCTTGGCGTAGTCGCGATGGGCTATGGCGATGGTTATCCGCGCGCCGCGCCGTCCGGTACGCCAG  >  NC_000913/4266378‑4266643

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