Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	568,250	Δ1 bp	7.8%	coding (1246/2835 nt)	uvrA ←	uvrABC excinuclease protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	568,250	0	C	.	7.8%	160.0 / 7.4	51	coding (1246/2835 nt)	uvrA	uvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base . (2/2);  total (28/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 = 7.97e-01

GTTCGCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGA  >  NC_002947/568103‑568396                                                                                                                                                    |                                                                                                                                                    gTTCGCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTcg                                                                                                                                                    >  1:251680/1‑149 (MQ=255)     gCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCgcggc                                                                                                                                                >  2:206540/1‑149 (MQ=255)       aGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCCGCCAGTAACGGTCAGATGACCAAAGTAGTTGATGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCgcggcg                                                                                                                                               >  3:134569/1‑148 (MQ=255)            tCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgc                                                                                                                                         <  4:102592/148‑1 (MQ=255)              cTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgcgc                                                                                                                                       <  2:76349/148‑1 (MQ=255)                 gAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGt                                                                                                                                     <  2:251680/147‑1 (MQ=255)                  aGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCT‑GCGCCGCGAGGcgcgcca                                                                                                                                   <  4:244322/148‑4 (MQ=255)                      tttttGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCgg                                                                                                                              >  1:125865/1‑149 (MQ=255)                       ttttGCCGCGATTTCGCCACGCCTGCCAGTAAGGTTCAGTTGACCAATGTAGTTGCTGGCTTCACGGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCTCCTCGCGGCGC‑AGGCGCGTACCCAgg                                                                                                                              <  4:134569/148‑1 (MQ=255)                                  tttCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCAc                                                                                                                  >  5:163150/1‑149 (MQ=255)                                     cGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgcg                                                                                                                                        <  6:157548/124‑1 (MQ=255)                                     cGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGgcg                                                                                                                                        >  5:157548/1‑124 (MQ=255)                                      gCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGc                                                                                                                             <  1:97400/134‑1 (MQ=255)                                      gCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGc                                                                                                                             >  2:97400/1‑134 (MQ=255)                                       ccACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGTGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTg                                                                                                             >  7:120355/1‑149 (MQ=255)                                               gCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGTGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCaa                                                                                                     <  8:120355/149‑1 (MQ=255)                                                  aGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAga                                                                                                   >  1:95990/1‑147 (MQ=255)                                                            aGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGT‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAg                                                                                        <  1:206540/149‑1 (MQ=255)                                                            aGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGTCCAg                                                                                        >  6:256106/1‑149 (MQ=255)                                                                 aCCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGttct                                                                                    >  4:99332/1‑148 (MQ=255)                                                                  ccAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGttctt                                                                                   >  3:194427/1‑148 (MQ=255)                                                                   cAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTcgc                                                                                                                                                   >  5:15879/1‑84 (MQ=255)                                                                   cAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTcgc                                                                                                                                                   <  6:15879/84‑1 (MQ=255)                                                                   cAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTca                                                                                 >  8:37595/1‑149 (MQ=255)                                                                   cAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTca                                                                                 <  2:28351/149‑1 (MQ=255)                                                                    aaaGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTcac                                                                                >  1:36470/1‑149 (MQ=255)                                                                       gtagtTGCTGGCTTCACCGATCGGCAAATTGGGTACCGCCGGCAGGGTTTTCTCGCCCACCCACAGATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTcacgca                                                                             <  1:177786/149‑1 (MQ=255)                                                                                      aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGc                                                                                                        >  1:33579/1‑107 (MQ=255)                                                                                      aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGc                                                                                                        <  2:33579/107‑1 (MQ=255)                                                                                      aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCgg                                                              >  3:121712/1‑149 (MQ=255)                                                                                      aCCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGGGGCCGATTCgg                                                              >  8:208831/1‑149 (MQ=255)                                                                                              gCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                      >  4:144269/1‑149 (MQ=255)                                                                                              gCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                      >  6:183523/1‑149 (MQ=255)                                                                                               cAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCgg                                                      >  6:12141/1‑148 (MQ=255)                                                                                                      gTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAgcggcg                                              <  5:12141/149‑1 (MQ=255)                                                                                                       ttACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGgcgc                                             <  1:72680/149‑1 (MQ=255)                                                                                                           cgccgGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCaa                                         <  4:194427/149‑1 (MQ=255)                                                                                                              cgGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGt                                      <  7:37595/149‑1 (MQ=255)                                                                                                                    ggTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGca                                >  5:252406/1‑149 (MQ=255)                                                                                                                              cccacccaCACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTc                      >  4:24436/1‑149 (MQ=255)                                                                                                                              cccacccaCACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTc                      <  1:72658/149‑1 (MQ=255)                                                                                                                              cccacccaCACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCAAGCACGTTGGCCGATTCGGTTTCGCTGTCTCGGCGCTCAAGCTACGGCACACTGCCTTc                      >  8:114303/1‑149 (MQ=255)                                                                                                                                   ccacACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACg                 >  2:115421/1‑149 (MQ=255)                                                                                                                                    cacacaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACgg                <  3:144269/149‑1 (MQ=255)                                                                                                                                     acacaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCTCGGTGGCCGATTCGGTTTCGCGGTAGCGGAGCTCAAGGTTCGGCACAATGCCTTCGAACggg               >  6:188890/1‑149 (MQ=255)                                                                                                                                     acacaTGGCGCGCCTCGCGGCGC‑AGGAGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACggg               >  1:11399/1‑149 (MQ=255)                                                                                                                                      cacaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGt              <  3:99332/149‑1 (MQ=255)                                                                                                                                       acaTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTg             <  2:36470/149‑1 (MQ=255)                                                                                                                                        gggtacgcgcCT‑GCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGAGTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGc            >  7:147705/6‑149 (MQ=255)                                                                                                                                         ggtacgcgcCT‑GCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAAc                  <  7:125732/138‑1 (MQ=255)                                                                                                                                         ggtacgcgcCT‑GCGCCGCGAGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAAc                  >  8:125732/5‑142 (MQ=255)                                                                                                                                         aTGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGc            >  1:169521/1‑148 (MQ=255)                                                                                                                                                 ccTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTg   <  5:256106/149‑1 (MQ=255)                                                                                                                                                  cTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGa  <  6:163150/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GTTCGCAGATTTCCTTGAGGATTTTTGCCGCGATTTCGCCACGCCGGCCAGTAAGGGTCAGTTGACCAAAGTAGTTGCTGGCTTCACCGATCGGCAAATTGGTTACCGCCGGCAGGGTTTTCTCGCCCACCCACACATGGCGCGCCTCGCGGCGC‑AGGCGCGTACCCCGGCAATCCGGGCACGGTTGCGTGCCAAGGAATTTGGCCAGTTCTTCACGCACGGTGGCCGATTCGGTTTCGCGGTAGCGGCGCTCAAGGTTCGGCACAATGCCTTCGAACGGGTGCGATCGCTTGA  >  NC_002947/568103‑568396

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: ‑