Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,643,899	G→A	G428G (GGC→GGT)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,899	0	G	A	100.0%	222.6 / NA	68	G428G (GGC→GGT)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/34);  total (34/34)

ATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCA  >  NC_000913/3643789‑3644031                                                                                                               |                                                                                                                                     aTCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCccgc                                                                                                          >  2:170554/1‑139 (MQ=255)     cTTCGTCGTGAGTTAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCcgcg                                                                                                      <  1:1129220/139‑1 (MQ=255)        cgtcgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                              <  1:39659/112‑1 (MQ=255)        cgtcgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                              >  2:39659/1‑112 (MQ=255)        cgtcgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCtt                                                                                                   >  2:1054081/1‑139 (MQ=255)           cgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTcatc                                                                                                                        >  1:1620844/1‑115 (MQ=255)           cgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTcatc                                                                                                                        <  2:1620844/115‑1 (MQ=255)             tgagtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                              >  2:1064443/1‑139 (MQ=255)               agtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            >  2:1182417/1‑139 (MQ=255)               agtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            >  2:49033/1‑139 (MQ=255)                gtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCAc                                                                                           >  1:1040209/1‑139 (MQ=255)                       gcgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACACGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                    <  1:1064443/139‑1 (MQ=255)                       gcgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                    <  2:1040209/139‑1 (MQ=255)                        cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            >  2:387060/1‑130 (MQ=255)                        cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            <  1:387060/130‑1 (MQ=255)                             gTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATc                                                                              <  1:490577/139‑1 (MQ=255)                                          cGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATc                                                                              <  1:1090002/126‑1 (MQ=255)                                          cGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATc                                                                              >  2:1090002/1‑126 (MQ=255)                                            ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                               >  2:481983/1‑123 (MQ=255)                                            ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                               <  1:481983/123‑1 (MQ=255)                                            ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTAcc                                                               >  2:532017/1‑139 (MQ=255)                                             gCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcg                                                              >  2:781733/1‑139 (MQ=255)                                                      gTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCt                                                                                                                                   <  1:792338/61‑1 (MQ=255)                                                      gTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCt                                                                                                                                   >  2:792338/1‑61 (MQ=255)                                                         aCAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            <  2:77108/97‑1 (MQ=255)                                                         aCAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            >  1:77108/1‑97 (MQ=255)                                                          cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                 <  2:1227086/139‑1 (MQ=255)                                                          cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                 <  1:1314413/139‑1 (MQ=255)                                                          cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                 <  1:532017/139‑1 (MQ=255)                                                           aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                              >  1:49532/1‑93 (MQ=255)                                                           aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                              <  2:49532/93‑1 (MQ=255)                                                           aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCg                                                                             >  2:1322743/1‑110 (MQ=255)                                                           aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCg                                                                             <  1:1322743/110‑1 (MQ=255)                                                           aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAtt                                                >  2:1664546/1‑139 (MQ=255)                                                               cAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGTCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCg                                            <  1:1182417/139‑1 (MQ=255)                                                                   tACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCAt                                        <  1:781733/139‑1 (MQ=255)                                                                     cgcgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                               >  1:750905/1‑98 (MQ=255)                                                                     cgcgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                               <  2:750905/98‑1 (MQ=255)                                                                        gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                          <  1:1209978/84‑1 (MQ=255)                                                                        gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                          >  2:1209978/1‑84 (MQ=255)                                                                        gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCa                                                                                      >  1:1215075/1‑88 (MQ=255)                                                                        gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCa                                                                                      <  2:1215075/88‑1 (MQ=255)                                                                        gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                       >  1:1139817/1‑135 (MQ=255)                                                                        gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                       <  2:1139817/135‑1 (MQ=255)                                                                         aCCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCccg                                                                                                           >  1:627143/1‑66 (MQ=255)                                                                         aCCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCccg                                                                                                           <  2:627143/66‑1 (MQ=255)                                                                          ccGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTc                                 >  1:1013677/1‑139 (MQ=255)                                                                             gTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCATTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCa                                                                                      >  2:1117258/1‑83 (MQ=255)                                                                             gTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGTTTGTTTTCACGGGCa                                                                                      <  1:1117258/83‑1 (MQ=255)                                                                                   gCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACg                        <  1:1054081/139‑1 (MQ=255)                                                                                   gCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACg                        <  1:1664546/139‑1 (MQ=255)                                                                                        gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                          <  2:1546191/68‑1 (MQ=255)                                                                                        gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                          >  1:1546191/1‑68 (MQ=255)                                                                                        gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcgcg                                                            >  1:345266/1‑98 (MQ=255)                                                                                        gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcgcg                                                            <  2:345266/98‑1 (MQ=255)                                                                                            cATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGtt                                                       <  1:1405595/99‑1 (MQ=255)                                                                                            cATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGtt                                                       >  2:1405595/1‑99 (MQ=255)                                                                                              tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            >  1:1383187/1‑60 (MQ=255)                                                                                              tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                            <  2:1383187/60‑1 (MQ=255)                                                                                                 gCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGTGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGcc          >  1:1114866/1‑139 (MQ=255)                                                                                                   tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                               >  2:1645891/1‑68 (MQ=255)                                                                                                   tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                               <  1:1645891/68‑1 (MQ=255)                                                                                                   tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCAGACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAg        >  1:1079573/1‑139 (MQ=255)                                                                                                   tttACGCATCTGACCTACGCAGTCAGCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAg        <  2:1079573/139‑1 (MQ=255)                                                                                                       cGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                       >  1:1574428/1‑104 (MQ=255)                                                                                                       cGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                       <  2:1574428/104‑1 (MQ=255)                                                                                                       cGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACAt    >  2:1377421/1‑139 (MQ=255)                                                                                                         cATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCTCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCa  <  1:170554/139‑1 (MQ=255)                                                                                                               |                                                                                                                                     ATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCA  >  NC_000913/3643789‑3644031

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