Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,805,092	Δ1 bp	7.3%	coding (67/618 nt)	PP_3363 ←	thymidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,805,092	0	C	.	7.3%	175.1 / 7.7	55	coding (67/618 nt)	PP_3363	thymidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base . (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

TTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGATTCTCGACGGGAGTA  >  NC_002947/3804946‑3805235                                                                                                                                                    |                                                                                                                                                ttCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCt                                                                                                                                               <  2:287224/149‑1 (MQ=255)    gCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATAACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTc                                                                                                                                              >  3:287442/1‑147 (MQ=255)    gCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTGCGGGAAAGCTGGGGGACCAGGGCCATGGGCTCGCCCCGGTGAGGAGCATTCTTTCTTTCGCAACGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGGGACGGCCTCca                                                                                                                                            <  6:291998/149‑1 (MQ=255)     aggaCCTTGAGTTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCcag                                                                                                                                           >  2:70627/4‑149 (MQ=255)     aggaCCTTGAGTTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCcag                                                                                                                                           >  4:61718/4‑149 (MQ=255)          cTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCg                                                                                                                                      >  5:352425/1‑149 (MQ=255)          cTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCg                                                                                                                                      >  6:342118/1‑149 (MQ=255)             gAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTg                                                                                                                                    <  5:342118/148‑1 (MQ=255)                   aaGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCtt                                                                                                                               <  6:141349/147‑1 (MQ=255)                     ggTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCg                                                                                                                           <  7:166834/149‑1 (MQ=255)                                gggCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTT‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTggg                                                                                                                <  1:70627/149‑1 (MQ=255)                                 ggCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGc                                                                                                               >  4:75693/1‑149 (MQ=255)                                    ttCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCa                                                                                                            >  7:29146/1‑149 (MQ=255)                                    ttCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGG‑CTCCAGCAGCGTGGTCTTGCCGGTACCCTTgg                                                                                                                 >  4:187402/1‑143 (MQ=255)                                    ttCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGG‑CTCCAGCAGCGTGGTCTTGCCGGTACCCTTgg                                                                                                                 >  4:181723/1‑143 (MQ=255)                                    ttCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGG‑CTCCAGCAGCGTGGTCTTGCCGGTACCCTTgg                                                                                                                 <  3:187402/143‑1 (MQ=255)                                    ttCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGG‑CTCCAGCAGCGTGGTCTTGCCGGTACCCTTgg                                                                                                                 <  3:181723/143‑1 (MQ=255)                                      cTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCAt                                                                                                           <  3:61718/148‑1 (MQ=255)                                            aGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTATCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAgg                                                                                                     >  2:122446/1‑148 (MQ=255)                                                  ttGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATaaa                                                                                              <  6:314358/149‑1 (MQ=255)                                                       cAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcg                                                                                         >  4:172732/1‑149 (MQ=255)                                                         gggCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcggg                                                                                       <  7:957/149‑1 (MQ=255)                                                            ccATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcgggcgg                                                                                    <  8:29146/149‑1 (MQ=255)                                                                  tctcGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCTGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGAAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATc                                                                              <  4:287442/149‑1 (MQ=255)                                                                      gCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTg                                                                           >  4:31642/1‑148 (MQ=255)                                                                      gCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGc                                                                          <  3:75693/149‑1 (MQ=255)                                                                         ccGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAAt                                                                       <  3:172732/149‑1 (MQ=255)                                                                               aGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGt                                                                 <  8:213395/149‑1 (MQ=255)                                                                                  aTCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTgg                                                              >  7:66891/1‑149 (MQ=255)                                                                                    cATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGt                                                                            <  4:43029/133‑1 (MQ=255)                                                                                    cATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGt                                                                            >  3:43029/1‑133 (MQ=255)                                                                                            ctttgGCACAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGGGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCCTCGTGCAATACCAGTTGGCTCAAgcgc                                                     <  8:150185/148‑1 (MQ=255)                                                                                            ctttcGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAgcgc                                                     <  1:122446/148‑1 (MQ=255)                                                                                               tcGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGCAGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTggg                                                                                                                <  1:257518/87‑1 (MQ=255)                                                                                               tcGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTggg                                                                                                                >  2:257518/1‑86 (MQ=255)                                                                                                cGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGa                                                 >  1:1354/1‑148 (MQ=255)                                                                                                      ggCGGATCACCTTGTTCCCGTCCAGCCTTAACGCTTG‑AGTGACGGCCTCCAGCAGCGGGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGAcc                                          <  6:352425/149‑1 (MQ=255)                                                                                                          gATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCAtt                                       >  8:48818/1‑148 (MQ=255)                                                                                                               ccTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCg                                  >  6:63621/1‑148 (MQ=255)                                                                                                                    ttCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTg                            >  3:19799/1‑149 (MQ=255)                                                                                                                    ttCCCGTCCAGCCGTAACGATTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTg                            >  4:266935/1‑149 (MQ=255)                                                                                                                       ccGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCt                         >  7:212899/1‑149 (MQ=255)                                                                                                                            cAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCaa                     <  4:19799/148‑1 (MQ=255)                                                                                                                               ccGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGa                 >  8:37730/1‑149 (MQ=255)                                                                                                                                cGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGAt                >  5:206507/1‑149 (MQ=255)                                                                                                                                 gTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGAtt               <  8:66891/149‑1 (MQ=255)                                                                                                                                    aCGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGATtct             <  7:37730/148‑1 (MQ=255)                                                                                                                                     cGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGATTCTCg           <  5:63621/149‑1 (MQ=255)                                                                                                                                      gCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGATTCTCGa          >  5:137722/1‑149 (MQ=255)                                                                                                                                            aGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCggg                                                  >  8:171099/1‑104 (MQ=255)                                                                                                                                            aGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCggg                                                  <  7:171099/104‑1 (MQ=255)                                                                                                                                             gTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGATTCTCGACGGGa     <  3:31642/148‑1 (MQ=255)                                                                                                                                                aCGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATAc                                                                     <  7:211480/81‑1 (MQ=255)                                                                                                                                                aCGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATAc                                                                     >  8:211480/1‑81 (MQ=255)                                                                                                                                                aCGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGATTCTCGACGGGAGTa  <  6:206507/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTG‑AGTGACGGCCTCCAGCAGCGTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGGTTCATCGTGCAATACCAGTTGGCTCAAGCGCGGGATAAGACCATTGATCGTAGCTGCCTTCAAGCGATTCTCGACGGGAGTA  >  NC_002947/3804946‑3805235

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