Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	544,754	C→T	R180H (CGT→CAT)	allC ←	allantoate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	544,754	0	C	T	100.0%	218.1 / NA	67	R180H (CGT→CAT)	allC	allantoate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/34);  total (33/34)

GTTTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCA  >  NC_000913/544618‑544887                                                                                                                                         |                                                                                                                                      gTTTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGAtgc                                                                                                                                     >  1:1712312/1‑139 (MQ=255)             gTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTggg                                                                                                                         <  1:1394005/139‑1 (MQ=255)                cAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCg                                                                                                                      <  2:735814/139‑1 (MQ=255)                     ttACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCgg                                                                                                                 <  2:597906/139‑1 (MQ=255)                           tataACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGt                                                                                                           <  2:1712312/139‑1 (MQ=255)                             taACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTGGTGGGGCGTTCGGAAGAGTaa                                                                                                         >  1:648872/1‑139 (MQ=255)                               acgacgCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                           <  2:76370/119‑1 (MQ=255)                               acgacgCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                           >  1:76370/1‑119 (MQ=255)                                  acgCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCg                                                                                                    <  1:229885/139‑1 (MQ=255)                                    gCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                  <  1:1518109/139‑1 (MQ=255)                                    gCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                  <  1:946120/139‑1 (MQ=255)                                       gCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGc                                                                                               >  1:1749884/1‑139 (MQ=255)                                                 tGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTGGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgca                                                                                     <  2:648872/139‑1 (MQ=255)                                                 tGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgca                                                                                     <  1:422919/139‑1 (MQ=255)                                                               cAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTAtt                                                                       >  1:725380/1‑139 (MQ=255)                                                                  ttgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGCGTTAGTGGGGCGTTCGGAAGAGTaaa                                                                                                        >  1:950303/1‑103 (MQ=255)                                                                  ttgattgCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTaaa                                                                                                        <  2:950303/103‑1 (MQ=255)                                                                          ccATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                            <  2:1749884/139‑1 (MQ=255)                                                                          ccATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                            <  2:1075845/139‑1 (MQ=255)                                                                              tACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATcaca                                                        >  2:1526367/1‑139 (MQ=255)                                                                                   tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGatat                                                   >  2:1394566/1‑139 (MQ=255)                                                                                   tCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGatat                                                   >  1:1340613/1‑139 (MQ=255)                                                                                          cacaGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAc                                            >  1:1256168/1‑139 (MQ=255)                                                                                             aGCCCTGTTGAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               <  2:1095856/101‑1 (MQ=255)                                                                                             aGCCCTGTTCAATATGCAGTTCAACAACGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               >  1:1095856/1‑101 (MQ=255)                                                                                             aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               >  2:450485/1‑101 (MQ=255)                                                                                             aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               <  1:450485/101‑1 (MQ=255)                                                                                             aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCa                                                             >  1:827227/1‑119 (MQ=255)                                                                                             aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCa                                                             <  2:827227/119‑1 (MQ=255)                                                                                             aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtc                                         >  2:638046/1‑139 (MQ=255)                                                                                             aGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATACGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtc                                         >  2:811793/1‑139 (MQ=255)                                                                                              gCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATca                                                          <  1:1356821/121‑1 (MQ=255)                                                                                              gCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATca                                                          >  2:1356821/1‑121 (MQ=255)                                                                                               cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                           <  1:1324732/55‑1 (MQ=255)                                                                                               cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTg                                                                                                                           >  2:1324732/1‑55 (MQ=255)                                                                                               cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTa                                                                                                          <  1:897742/72‑1 (MQ=255)                                                                                               cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTa                                                                                                          >  2:897742/1‑72 (MQ=255)                                                                                               cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtcgt                                       >  1:1387/1‑139 (MQ=255)                                                                                               cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtcgt                                       <  2:725380/139‑1 (MQ=255)                                                                                               cccTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCAcgtcgt                                       <  2:1387/139‑1 (MQ=255)                                                                                                    ttCAATATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGa                                  >  1:565970/1‑139 (MQ=255)                                                                                                       aaTATGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTc                               >  2:464563/1‑139 (MQ=255)                                                                                                         tatGCGGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                             >  1:409970/1‑139 (MQ=255)                                                                                                         tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                  <  1:68943/70‑1 (MQ=255)                                                                                                         tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCa                                                                                                  >  2:68943/1‑70 (MQ=255)                                                                                                         tatGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGc                             <  2:565970/139‑1 (MQ=255)                                                                                                           tGCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCa                           >  2:329715/1‑139 (MQ=255)                                                                                                            gCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                            >  2:1050183/1‑105 (MQ=255)                                                                                                            gCAGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                            <  1:1050183/105‑1 (MQ=255)                                                                                                              aGTTCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGcc                        >  1:1168624/1‑139 (MQ=255)                                                                                                                ttCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCa                      >  2:373372/1‑139 (MQ=255)                                                                                                                ttCAACAAAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCa                      >  1:1715290/1‑139 (MQ=255)                                                                                                                  caacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCaaa                    <  1:329715/139‑1 (MQ=255)                                                                                                                   aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAACACTATTTCCTTTGGCAt                                                            <  2:1147834/98‑1 (MQ=255)                                                                                                                   aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCAt                                                            >  1:1147834/1‑98 (MQ=255)                                                                                                                   aacaaAGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCaaaa                   >  2:609834/1‑139 (MQ=255)                                                                                                                        aGGCTTTAATATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATAtt              <  2:1256168/139‑1 (MQ=255)                                                                                                                        aGGCTTTAATATCCTGATGAGAAGTTAGTGGGGCGTTCGGAAGAGTAACTTCTCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTCGCATCACAGATCATCCGCACGTCGTCAGGATTCCCCAGCCCAAAAATAtt              >  2:51785/1‑139 (MQ=255)                                                                                                                          gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    >  2:508573/1‑35 (MQ=255)                                                                                                                          gCTTTAATATCCTGATGCGGAGTTAGTGGGGCGtt                                                                                                                    <  1:508573/35‑1 (MQ=255)                                                                                                                               aaTATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTg       <  1:464563/139‑1 (MQ=255)                                                                                                                               aaTATCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTg       <  2:1715290/139‑1 (MQ=255)                                                                                                                                 tatCCTGATGCGGAGTTAGGGGGGCGTCCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGcc     <  1:811793/139‑1 (MQ=255)                                                                                                                                   tCCTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGcccc   <  1:638046/139‑1 (MQ=255)                                                                                                                                    ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgcat                                                                                    <  2:965576/57‑1 (MQ=255)                                                                                                                                    ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTcatcgcat                                                                                    >  1:965576/1‑57 (MQ=255)                                                                                                                                    ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               <  2:1162579/62‑1 (MQ=255)                                                                                                                                    ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACa                                                                               >  1:1162579/1‑62 (MQ=255)                                                                                                                                    ccTGATGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTTCCCCa  <  1:1394566/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      GTTTGATTCGCCGTTCAGCGTTACCGTATAACGACGCTGCCCGACAATTGCATTCACCACGCCAATTGATTGCCCATTACTTTCCAGCACACAGCCCTGTTCAATATGCAGTTCAACAAAGGCTTTAATATCCTGACGCGGAGTTAGTGGGGCGTTCGGAAGAGTAAATCCGCAAGCCTTCATCGCATCGACAAAACTATTTCCTTTGGCATCACAGATATTCCGCACGTCGTCAGGATTCGCCAGCCCAAAAATATTTTTACTGCCCCA  >  NC_000913/544618‑544887

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