Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	961,535	C→T	A112V (GCA→GTA)	cmk →	cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	961,535	0	C	T	98.1%	162.0 / NA	52	A112V (GCA→GTA)	cmk	cytidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/23);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCG  >  NC_000913/961400‑961670                                                                                                                                        |                                                                                                                                        tGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTatta                                                                                                                                     >  2:967695/1‑140 (MQ=255) tGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTatta                                                                                                                                     >  2:214702/1‑140 (MQ=255)       tACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCCCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGt                                                                                                                               >  2:1021882/1‑140 (MQ=255)              tgtgTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGACTGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAAc                                                                                                                        >  1:207603/1‑140 (MQ=255)                  tCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAAcgcgc                                                                                                                    >  1:865829/1‑140 (MQ=255)                        aaTGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTcg                                                                                                              <  2:53102/140‑1 (MQ=255)                         aTGGCAATCTGGAAGTGATCCTCGAAGGGGGAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTcgc                                                                                                             <  2:518300/140‑1 (MQ=255)                          tGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAATCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTcgcg                                                                                                            >  2:610218/1‑140 (MQ=255)                           ggCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGa                                                                                                           <  1:610218/140‑1 (MQ=255)                           ggCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGa                                                                                                           <  2:568502/140‑1 (MQ=255)                           ggCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGa                                                                                                           <  2:983771/140‑1 (MQ=255)                               aTCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTa                                                                                                       >  1:998773/1‑140 (MQ=255)                               aTCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTa                                                                                                       <  1:214702/140‑1 (MQ=255)                                    gAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGGTCGTGAAGCATTATTGGGTCGCCAAGGCGCGTTTCGCGAATTACCAgg                                                                                                  >  2:389244/1‑140 (MQ=255)                                     aaGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATACGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGt                                                                                                 >  1:650393/1‑140 (MQ=255)                                      aGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTc                                                                                                <  2:213365/140‑1 (MQ=255)                                      aGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTc                                                                                                >  1:841008/1‑140 (MQ=255)                                         gATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGa                                                                                             <  2:1011947/140‑1 (MQ=255)                                          aTCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGCAGTATTATAGCGGCGCTAACGCGCGGTTCGGAAATTACCGGGTCTGAt                                                                                            >  2:496276/1‑140 (MQ=255)                                            ccTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTg                                                                                          >  2:617822/1‑140 (MQ=255)                                            ccTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCACCATCGCGCGTTTCGCGTATTACCAGGTCTGATTg                                                                                          >  2:390519/1‑140 (MQ=255)                                                gAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGa                                                                                      <  2:852076/140‑1 (MQ=255)                                                    gggAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGc                                                                                  <  2:226899/140‑1 (MQ=255)                                                       aaGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGTGTTTCGCGAATAACGAGGTCTAATTGCCGATGGCcgc                                                                               >  1:260091/1‑140 (MQ=255)                                                       aaGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCcgc                                                                               <  2:865829/140‑1 (MQ=255)                                                        aGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCcgcg                                                                              >  2:157525/1‑140 (MQ=255)                                                           tGTCAGCGGCGAAATCCGTACTTAGGAAGTGGCGAATACAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACa                                                                           <  1:496276/140‑1 (MQ=255)                                                            gTCAGCGGCGAAATTGGTACTCAGGTAGTAGCGAATGGAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACAt                                                                          <  1:390519/140‑1 (MQ=255)                                                            gTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACAt                                                                          >  1:778297/1‑140 (MQ=255)                                                              cAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGGCTGATTGCCGATGGCCGCGACATgg                                                                        >  2:956059/1‑140 (MQ=255)                                                               aGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATggg                                                                       <  2:778297/140‑1 (MQ=255)                                                                 cggcgAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGaa                                                                     >  1:657017/1‑140 (MQ=255)                                                                 cggcgAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGaa                                                                     <  1:169704/140‑1 (MQ=255)                                                                 cggcgAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCACCAACGCACGTTTCGCTAATGACCAGGCCTGATAGCCGTTGGCCGCGACATGGGaa                                                                     >  2:282647/1‑140 (MQ=255)                                                                 cggcgAAATTCGTACTCAGGAACTGGCGAATGCAGCTCCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGaa                                                                     <  1:463618/140‑1 (MQ=255)                                                                   gcgAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACg                                                                   <  2:95079/140‑1 (MQ=255)                                                                    cgAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACgg                                                                  >  1:3314/1‑140 (MQ=255)                                                                      aaaTTCGTACTCAGGGAGTGGCCAACGCATCTTCACATGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACggtg                                                                <  1:282647/140‑1 (MQ=255)                                                                      aaaTTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACggtg                                                                >  1:958905/1‑140 (MQ=255)                                                                         ttCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTa                                                             <  2:914885/140‑1 (MQ=255)                                                                                                    cTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGAAACGGTGATATACCCTGATGCCCCAGTGAAAATTTTc                                  >  1:310589/1‑140 (MQ=255)                                                                                                          aaGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGAc                            >  1:129238/1‑140 (MQ=255)                                                                                                          aaGTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGAc                            <  2:998773/140‑1 (MQ=255)                                                                                                            gTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGc                          <  2:650393/140‑1 (MQ=255)                                                                                                            gTCGCGGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGc                          >  2:1006178/1‑140 (MQ=255)                                                                                                              cgcgGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTGTTCCCTGCTGCACCCGTGAAAGTTTTCCTTGACGcct                        >  1:492132/1‑140 (MQ=255)                                                                                                              cgcgGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGcct                        >  2:634984/1‑140 (MQ=255)                                                                                                              cgcgGCATTCCCACGCGTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTAATGCACCAGTGAATATTTTCCTTGACGcct                        >  2:499807/1‑140 (MQ=255)                                                                                                                           cgcgTTCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACg           >  2:925166/1‑140 (MQ=255)                                                                                                                               ttCGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTgcg       >  1:131948/1‑140 (MQ=255)                                                                                                                                tCGTGAAGTACTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCCGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGGCTTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTgcgc      >  2:903350/1‑140 (MQ=255)                                                                                                                                 cGTGAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGTTGCACCAGGGAAAATTTTCCTTGACACCTCCACGGAAGAACGTGCGCa     >  2:468155/1‑140 (MQ=255)                                                                                                                                  gTGAAGTAATATTGCGTCGCCAACGCGCGTTCCTCAAAACACCAGGTCAGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCAt    <  2:492132/140‑1 (MQ=255)                                                                                                                                    gAAGTATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATcg  <  1:832871/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        TGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGATCCTCGAAGGGGAAGATGTCAGCGGCGAAATTCGTACTCAGGAAGTGGCGAATGCAGCTTCACAAGTCGCGGCATTCCCACGCGTTCGTGAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCG  >  NC_000913/961400‑961670

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