Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	359,776	T→A	R96R (CGT→CGA)	cynS →	cyanate aminohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	359,776	0	T	A	98.3%	224.6 / ‑3.2	58	R96R (CGT→CGA)	cynS	cyanate aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (22/35);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

TGGGTCAGCAGGCGCTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGTTTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACGTGCGGTCATCA  >  NC_000913/359631‑359921                                                                                                                                                  |                                                                                                                                                  tGGGTCAGCAGGCGCTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCt                                                                                                                                               >  2:45102/1‑150 (MQ=255)  gggTCAGCAGGCGCTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCt                                                                                                                                               <  2:328078/149‑1 (MQ=255)  gggTCAGCAGGCGCTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTa                                                                                                                                              <  2:657504/150‑1 (MQ=255)  gggTCAGCAGGCGCTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTAt                                                                                                                                             <  2:15374/151‑1 (MQ=255)  gggTCAGCAGGCGCTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTAt                                                                                                                                             <  2:120923/151‑1 (MQ=255)             cgcTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCa                                                                                                                                  <  2:321469/151‑1 (MQ=255)              gcTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAGTGTTGc                                                                                                                                   >  2:88881/1‑149 (MQ=255)               cTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAgg                                                                                                                                >  2:252869/1‑151 (MQ=255)                ttCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAgg                                                                                                                                >  1:512546/1‑150 (MQ=255)                    tGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTAc                                                                                                                           >  2:283764/1‑151 (MQ=255)                            ccgcccgcCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAAc                                                                                                                   >  1:56991/1‑151 (MQ=255)                             cgcccgccTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTATTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGGCCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAAcc                                                                                                                  <  1:452716/151‑1 (MQ=255)                             cgcccgccTGGTCGGGGCGAAGCTGGAACTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAAcc                                                                                                                  <  1:477855/151‑1 (MQ=255)                                ccgccTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTg                                                                                                               >  1:473817/1‑151 (MQ=255)                                         cGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTg                                                                                                      <  2:526791/151‑1 (MQ=255)                                         cGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTg                                                                                                      >  1:452121/1‑151 (MQ=255)                                            ggCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGtt                                                                                                   <  1:271710/151‑1 (MQ=255)                                               gAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCAt                                                                                                <  2:328446/151‑1 (MQ=255)                                                 aGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATga                                                                                              <  1:29871/151‑1 (MQ=255)                                                   cTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATg                                                                                               >  2:377945/1‑148 (MQ=255)                                                             acgcAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCg                                                                                   >  2:139902/1‑150 (MQ=255)                                                              cgaAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGAt                                                                                 <  1:547083/151‑1 (MQ=255)                                                                 aGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGc                                                                              <  2:744149/151‑1 (MQ=255)                                                                   aCTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCat                                                                            <  2:657035/151‑1 (MQ=255)                                                                     tCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCatt                                                                           >  1:7446/1‑150 (MQ=255)                                                                       cATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCattatt                                                                        <  2:423605/151‑1 (MQ=255)                                                                        aTTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCattatt                                                                        >  2:45624/1‑150 (MQ=255)                                                                        aTTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCAttatta                                                                       <  1:572516/151‑1 (MQ=255)                                                                             aCTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAgcgc                                                                   <  1:330478/150‑1 (MQ=255)                                                                             aCTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAgcg                                                                    <  1:568445/149‑1 (MQ=255)                                                                                gTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGAtt                                                               <  2:251384/151‑1 (MQ=255)                                                                                   gCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTaa                                                             <  1:38878/150‑1 (MQ=255)                                                                                     aGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACtt                                                          <  2:3584/151‑1 (MQ=255)                                                                                             ccACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTc                                                   <  2:427467/150‑1 (MQ=255)                                                                                             ccACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTc                                                   <  2:7446/150‑1 (MQ=255)                                                                                             ccACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCg                                                  <  2:650707/151‑1 (MQ=255)                                                                                             ccACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCg                                                  >  1:544882/1‑151 (MQ=255)                                                                                              cACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGa                                                 <  1:356858/151‑1 (MQ=255)                                                                                                  gCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGtt                                             <  1:74722/151‑1 (MQ=255)                                                                                                      ggCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAg                                          <  2:107679/150‑1 (MQ=255)                                                                                                      ggCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGa                                         <  2:466190/151‑1 (MQ=255)                                                                                                       gCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGaa                                        <  2:87129/151‑1 (MQ=255)                                                                                                        cTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGaa                                        >  1:278964/1‑150 (MQ=255)                                                                                                            aTTGATGACCGTATTCCAACTGACCCGACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTgg                                   <  2:300352/151‑1 (MQ=255)                                                                                                              tgatgaCCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGc                                  <  1:38066/150‑1 (MQ=255)                                                                                                                atgcCCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCgg                                >  2:143753/1‑150 (MQ=255)                                                                                                                atgaCCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCgg                                >  2:452878/1‑150 (MQ=255)                                                                                                                  gaCCGTATTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGAcc                             >  2:678119/1‑151 (MQ=255)                                                                                                                        aTTCCAACTGACCCAACGGTGTATCGATTCTATGAACTGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGAGGTTAAGAAAGTGGCGGACCCg                           >  2:204949/1‑147 (MQ=255)                                                                                                                        aTTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGaa                        <  1:127290/150‑1 (MQ=255)                                                                                                                        aTTCCAACTGACCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGa                         >  1:599178/1‑149 (MQ=255)                                                                                                                                  aCCCAACGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACg             >  2:124332/1‑151 (MQ=255)                                                                                                                                      aaCGATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACGTGCg         >  2:474951/1‑151 (MQ=255)                                                                                                                                        cGATGGATCGATTCTATGAAATGGTGCAGGTGTGCGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGAAGGCATTATTAGCGCGATTAACGTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGAGGCGGACGAGCgg        >  2:674766/1‑150 (MQ=255)                                                                                                                                         gATGTATCGTTTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACGTGCGGtc      <  1:526364/151‑1 (MQ=255)                                                                                                                                         gATGTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACGTGCGGtc      <  2:229696/151‑1 (MQ=255)                                                                                                                                            gTATCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACGTGCGGtcat    <  1:120190/150‑1 (MQ=255)                                                                                                                                              aTCGATTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACGTGCGGtcatca  <  2:742562/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TGGGTCAGCAGGCGCTTCCTGCCGACGCCGCCCGCCTGGTCGGGGCGAAGCTGGATCTCGACGAAGACTCCATTCTACTGTTGCAGATGATTCCACTGCGTGGCTGCATTGATGACCGTATTCCAACTGACCCAACGATGTATCGTTTCTATGAAATGTTGCAGGTGTACGGTACAACCCTGAAAGCGTTGGTTCATGAGAAATTTGGCGATGGCATTATTAGCGCGATTAACTTCAAACTCGACGTTAAGAAAGTGGCGGACCCGGAAGGTGGCGAACGTGCGGTCATCA  >  NC_000913/359631‑359921

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