Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	539,504	G→A	7.3%	A120T (GCG→ACG)	allB →	allantoinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	539,504	0	G	A	7.3%	136.1 / 2.2	57	A120T (GCG→ACG)	allB	allantoinase
Reads supporting (aligned to +/- strand):  ref base G (21/31);  new base A (2/2);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TATGAAACCGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTATCGACAACGACTTCCGTGATGTAAACGACT  >  NC_000913/539363‑539645                                                                                                                                              |                                                                                                                                               taggAAACCGGGACTCGCGCAGCGGCAAAAGGGGGTAGCACCGCCATGATCGAAATGCCGCTCAACCAGCTGGCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGc                                                                                                                                               <  2:365136/140‑1 (MQ=255) tATGAAACCGGTGCTCGCGCGGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCGGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGc                                                                                                                                               <  1:859785/143‑1 (MQ=255) tATGAAACCGGTACACGCGCAGCGGCAAAAGGTGGGAGCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGc                                                                                                                                               <  1:672921/143‑1 (MQ=255)         cGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGCTGCTGCACAAc                                                                                                                                       >  1:697786/1‑143 (MQ=255)         cGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTAAATTGAACTGAAGTTCGATGCCGCTAACGGCAAGCTGACTATTGATGCGGcaccac                                                                                                                                       >  1:736801/1‑143 (MQ=255)              gtttcggaGCGGCAATAGGTGGTATCACCACCATGATCGACATGCCGCTCAACCAGCTGCCGGCAACGGTTGACCGCGCTTCAATTGGACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGt                                                                                                                                  <  2:111284/136‑1 (MQ=255)              cTCGCGCAGCGGCAAGAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATACGGCACAACCCGGt                                                                                                                                  >  2:956613/1‑143 (MQ=255)              cTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGt                                                                                                                                  <  1:1095663/143‑1 (MQ=255)              cTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATACGGCACAACTCGGt                                                                                                                                  >  1:1058099/1‑143 (MQ=255)                    cAGCGGCAAACGGTGGTATCACCACCATGATCGACATGCCGCTCAACTAGCTGCCTGCAACGGATGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTCAAGCCAAGCTGACTATTGAGGCTGCGCGACTCGGGCGCCtc                                                                                                                            >  2:840586/1‑142 (MQ=255)                       cGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGgtg                                                                                                                         >  2:52608/1‑143 (MQ=255)                         gCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTAGCCTGGTGTc                                                                                                                       >  1:965400/1‑143 (MQ=255)                           aaaaGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCAAAACTCGGTGGCCTGGTGTCtt                                                                                                                     >  1:1104670/1‑143 (MQ=255)                            aaaGGTGGTATCACCACCATGCTCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTa                                                                                                                    <  1:349244/143‑1 (MQ=255)                            aaaGGTGGTATCACCACCATGATCGAAATGCCTCTCAACCAGCGGCCTGCAACGGTTGACCGCGCGTCAATTGAACTGAAGTTCGATGCCTCTAAAGGCGTGCTGTCTAGTGATGCGGCACAACTCGGTGGCCGGGTGTCTTa                                                                                                                    >  2:482322/1‑143 (MQ=255)                                  ggtATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑At                                                                                                             >  2:180421/1‑143 (MQ=255)                                  ggtATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑At                                                                                                             >  1:444754/1‑143 (MQ=255)                                   gtATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATc                                                                                                            <  1:185798/143‑1 (MQ=255)                                     cacatcactaTGATAGAAAAGCCGCTCAACCAGCTCCCTGCAACGTTTGACCGCGCTTCAATTGACGTGAAGTCGGAGTCCGCTAACGGCAAGCTTACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGa                                                                                                          <  2:134212/138‑1 (MQ=255)                                           accaTGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCCATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACACCCCGATGGCCTGGGTTCTTACGAC‑ATCGACCGCCt                                                                                                    >  2:1048766/1‑143 (MQ=255)                                                 aTCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAACTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGa                                                                                              >  1:568695/1‑143 (MQ=255)                                                   cGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGGCTTACAAC‑ATCGACCGTCTGCATGAGc                                                                                            >  2:920736/1‑143 (MQ=255)                                                    gAAATGCCGCTCAAGCAGCGGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCt                                                                                           <  1:753782/143‑1 (MQ=255)                                                    gAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCGAATGAACTGAAGTTCGATGCCGCTAAAGGCAAGCGGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCt                                                                                           <  2:329895/143‑1 (MQ=255)                                                    gAAAGGCCGGTCACTCAGCGGCCTGCAACGGTCGACCGGGCTGCAATTGAAAGGAAGTTCGATGCCGCGAACGGCAAGCTGACTATTGATGCGGCACAACCCGGCGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCt                                                                                           <  1:766520/143‑1 (MQ=255)                                                      aaTGGCGCGCGACCGGCTGCCTGCAACGGTTGACCGCGCTTCAATTGCACTCAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTgg                                                                                         <  2:33677/143‑1 (MQ=255)                                                      aaTGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTCTTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACCAC‑ATCGACCGTCGACATGAGCTgg                                                                                         >  1:289787/1‑143 (MQ=255)                                                            gCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACACGGATGCCTGGTGGCCTAGATC‑ATCGACCGTCTGCATGAGCTGGAGGAAg                                                                                   >  2:848238/1‑143 (MQ=255)                                                                     gCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTg                                                                          <  2:568695/143‑1 (MQ=255)                                                                      ctgccggcAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCGAGCTGACTATTGAGGCGGCACAACTCGGTGGCCTGGTGTCTCACAAC‑CTCGACCGTCTGCGTGAGCTGGGTGAAGGGGGcgtcgt                                                                         >  2:460676/1‑143 (MQ=255)                                                                        gcctgcAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCg                                                                       <  1:720955/143‑1 (MQ=255)                                                                         gctgcAACGGTTGAGCGCGCTCCAACTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCgg                                                                      <  1:482322/142‑1 (MQ=255)                                                                            gcAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTAGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCtt                                                                   <  2:965400/143‑1 (MQ=255)                                                                                 ggTTGATCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATGGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAAt                                                              <  1:955722/143‑1 (MQ=255)                                                                                      aCCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCttt                                                         >  1:305366/1‑142 (MQ=255)                                                                                           tcTTCAATGATACTGAAGTGCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGTCCTGGTGTCTTACAGC‑ATCGACCTTCGGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGc                                                    <  1:183218/142‑1 (MQ=255)                                                                                           gcTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCCGGATCAAGTGGGCGCTGTCGGCTTCAAACGCTTCGTTGc                                                    >  2:481097/1‑143 (MQ=255)                                                                                            cTTCAATTGAAGTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCTGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCg                                                   <  1:442658/143‑1 (MQ=255)                                                                                            cTTCAATGGAACTGAAGTTCGATGCCGCTAAAGGCAAGCGGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCg                                                   <  1:429026/143‑1 (MQ=255)                                                                                                aaTTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCt                                               >  2:692713/1‑143 (MQ=255)                                                                                                    gAACGGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTgg                                           <  1:443370/143‑1 (MQ=255)                                                                                                     aaCTGAAGTTCGAGGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGc                                          <  1:833740/143‑1 (MQ=255)                                                                                                       gtGAAGTTCGATGCGGGTAAAGGCAAGATGACTATTGACGCGGCACAGCTCGGTGGCCTGGTGTCTGACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCGTCAAATGCGTCGTTGCGACCTGTGGCGa                                        <  2:1005380/142‑1 (MQ=255)                                                                                                        tGAAGTTCGATGCCGCTAAGGGCAAGCTGACTATTGGTACGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGAt                                       <  2:1058099/143‑1 (MQ=255)                                                                                                        tGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATACGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGAt                                       <  1:956613/143‑1 (MQ=255)                                                                                                           aGTTCGATGACGCTAAAGGCAAGCTGACTATTGATGCTGCATAACTCGGTGGCCTGGTGTCTT‑CATCTTTCTACCTGCTCCACGAGATGATTTATATTGTCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATcgc                                    <  1:262183/143‑1 (MQ=255)                                                                                                            gttgggtGCCGCCAAAGGTCAGCCGACTACTGATCCGGCGCAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGAGGAAGTGGGCATTTCCGTGTTTAAATGCTTCGTTGCGACCTGTGGCGATcgcg                                   <  2:1107817/137‑1 (MQ=255)                                                                                                             ttCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCgg                                  >  1:614509/1‑143 (MQ=255)                                                                                                             ttCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACAGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGCCCTGTCGCGATCGCgg                                  >  1:296868/1‑143 (MQ=255)                                                                                                             tgtgATGGCGCGAAAGGCGAGCTGCCGAGTGAGGCGGCACATCTCGGTGGCCGGGTGTCTTACAAC‑ATCGTCCGTCTGCATGTGCTGGATGAAGTGGGCGTTGTCGGCTTCAAGTGCTTCGTTGCGACCTGTGGCGATCGCgg                                  <  2:887960/140‑1 (MQ=255)                                                                                                                gagggctcTAGAGGCAAGCTGACTGTTGATGCGGCACAACTCGGTGGCCTTGTGGCTTAGAAC‑ATCGACCGTCTGTATGAGCTGGACGACGGGGGCGTCGTCGGCTTCAAATGCCTCGTGGCGACCTTTGGCGATCGCGGTAt                               <  2:254017/136‑1 (MQ=255)                                                                                                                 aTGCCGCTAAAGGCAAGCTGACTATTTACGCGGCACAACTCGGTCGCCTGGTGTCTTACAAC‑ATCGACCGCCTGCATGAGCTGGATCACGTGCCCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTATc                              <  2:835743/143‑1 (MQ=255)                                                                                                                             gCAAGCTGACTATGGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTATCGACAACGACTTc                  <  1:283272/143‑1 (MQ=255)                                                                                                                                aGCGGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTATCGACAACGACTTCCGt               <  1:166961/143‑1 (MQ=255)                                                                                                                                   tGACCATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTATCGACAACGACTTCCGTGAt            <  1:166519/143‑1 (MQ=255)                                                                                                                                       tATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTCTCGACAACGACTTCCGTGATGTaa        >  2:680312/1‑143 (MQ=255)                                                                                                                                             tGCGGCACAACTCGGGGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCAGGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTATCGACAACGACTTCCGTGATGTAAACGACt  <  1:236543/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               TATGAAACCGGTACTCGCGCAGCGGCAAAAGGTGGTATCACCACCATGATCGAAATGCCGCTCAACCAGCTGCCTGCAACGGTTGACCGCGCTTCAATTGAACTGAAGTTCGATGCCGCTAAAGGCAAGCTGACTATTGATGCGGCACAACTCGGTGGCCTGGTGTCTTACAAC‑ATCGACCGTCTGCATGAGCTGGATGAAGTGGGCGTTGTCGGCTTCAAATGCTTCGTTGCGACCTGTGGCGATCGCGGTATCGACAACGACTTCCGTGATGTAAACGACT  >  NC_000913/539363‑539645

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑