Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,431,880	T→C	8.9%	intergenic (‑9/‑41)	PP_1253 ← / → xenA	ArsR family transcriptional regulator/protein XenA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,431,880	0	T	C	8.9%	149.5 / 4.6	57	intergenic (‑9/‑41)	PP_1253/xenA	ArsR family transcriptional regulator/protein XenA
Reads supporting (aligned to +/- strand):  ref base T (21/30);  new base C (3/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCAA  >  NC_002947/1431736‑1432024                                                                                                                                                 |                                                                                                                                                   ccAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGATCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGa                                                                                                                                                <  8:175074/149‑1 (MQ=255)  cAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGAt                                                                                                                                               <  5:302776/149‑1 (MQ=255)  cAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGAt                                                                                                                                               <  8:268875/149‑1 (MQ=255)  cAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGAt                                                                                                                                               <  6:95382/149‑1 (MQ=255)      tGGCGGAAGATTTCGAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGACGGCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGaa                                                                                                                                           <  2:345502/149‑1 (MQ=255)           gACGATTTCCAGGCGCATGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAacaaca                                                                                                                                      <  7:104975/149‑1 (MQ=255)           gACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATCATAATCGaactatg                                                                                                                                      <  2:189643/149‑5 (MQ=255)            aCGATTTCCAGGCGCACGGGGTCGCTGAGCTCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCTCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAacaaca                                                                                                                                      <  1:134023/148‑1 (MQ=255)            aCGATTTCCAGGCGCACGCGGTCGCAGAGCGCATAGAGCAGGCGTTCAAGAATCAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCAAGGTTGACAAGAAAACGCGGGTTTCATTAATTGCCATAGTTCGATAAAGATCGaacaaaa                                                                                                                                      >  2:134023/1‑146 (MQ=255)             cGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATCATAATCGaactat                                                                                                                                       >  3:249214/1‑143 (MQ=255)             cGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATCATAATCGaactat                                                                                                                                       <  4:249214/146‑4 (MQ=255)                 ttCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTAtt                                                                                                                                 >  7:300413/1‑148 (MQ=255)                           aCGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAAccc                                                                                                                    <  5:75936/149‑1 (MQ=255)                                gTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAAACCGGAGt                                                                                                               >  6:221530/1‑149 (MQ=255)                                     tGAGCGCATAGAGCAGGCGTTCAAGAATGGGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTcc                                                                                                           >  2:262634/1‑148 (MQ=255)                                     tGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTcc                                                                                                           >  4:45051/1‑148 (MQ=255)                                      gAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGAttc                                                                                                                                                   >  6:252613/1‑108 (MQ=255)                                        gcgcATAGAGCAGGCGTTCAATAATGATCTCGTCGGGGTGTGGATGTTTATTGGCTCGCATTGTGGACATGATAACGGGGGTTTCATTAATTTCCATTGTTCGATTAGGATCGTACAACATTATTCAA‑A‑TGAACCCGGAGTTTCCCATg                                                                                                       <  6:173179/149‑1 (MQ=255)                                          gcATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTccc                                                                                                          >  5:212720/1‑144 (MQ=255)                                          gcATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTccc                                                                                                          <  6:212720/144‑1 (MQ=255)                                          gcATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTc                                                                                                     <  2:131340/149‑1 (MQ=255)                                            aTAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATTGAACAACAGTATTCAA‑A‑TGAACCCg                                                                                                                   <  5:168920/133‑1 (MQ=255)                                            aTAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATTGAACAACAGTATTCAA‑A‑TGAACCCg                                                                                                                   >  6:168920/1‑133 (MQ=255)                                               gagCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCa                                                                                                 <  7:267786/148‑1 (MQ=255)                                                       gTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAACA‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGa                                                                                         <  1:51600/149‑1 (MQ=255)                                                        ttCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAAc                                                                                       >  4:68860/1‑149 (MQ=255)                                                        ttCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAAc                                                                                       <  5:302549/149‑1 (MQ=255)                                                             gAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTa                                                                                   >  1:266119/1‑148 (MQ=255)                                                               aTGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTacac                                                                                <  1:89002/149‑1 (MQ=255)                                                                          ggggTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑ACTGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGAc                                                                      <  1:262634/149‑1 (MQ=255)                                                                          ggggTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACCCCCTCAAATACg                                                                     >  4:201824/1‑149 (MQ=255)                                                                                 ggATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCt                                                              >  5:148918/1‑149 (MQ=255)                                                                                   aTGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCTAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGc                                                            <  3:68860/149‑1 (MQ=255)                                                                                          taGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAAcc                                                      <  3:45051/148‑1 (MQ=255)                                                                                           aGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCg                                                     >  5:86519/1‑148 (MQ=255)                                                                                             gCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGAc                                                                      <  1:144646/129‑1 (MQ=255)                                                                                             gCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGAc                                                                      >  2:144646/1‑129 (MQ=255)                                                                                                     ggTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGc                                               >  8:196122/1‑144 (MQ=255)                                                                                                     ggTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGc                                               <  7:196122/144‑1 (MQ=255)                                                                                                       ttGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTccg                                        >  1:21601/1‑149 (MQ=255)                                                                                                            aTGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGcc                                              <  6:143175/138‑1 (MQ=255)                                                                                                            aTGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGcc                                              >  5:143175/1‑138 (MQ=255)                                                                                                              gATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGAtgt                                 >  5:270037/1‑149 (MQ=255)                                                                                                               aTAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGAtgt                                 <  5:221530/148‑1 (MQ=255)                                                                                                                tAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGc                               >  4:144742/1‑149 (MQ=255)                                                                                                                            cATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGc                    <  6:148918/148‑1 (MQ=255)                                                                                                                             aTTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                   <  8:300413/148‑1 (MQ=255)                                                                                                                             aTTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                   >  2:163698/1‑148 (MQ=255)                                                                                                                              ttAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCg                  >  2:282277/1‑148 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                   >  6:18329/1‑140 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                   <  5:18329/140‑1 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACgg           <  3:144742/148‑1 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACgg           <  6:270037/148‑1 (MQ=255)                                                                                                                                      ttttgTTCGATCATAATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACACTCAAATACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCa         >  7:220338/5‑149 (MQ=255)                                                                                                                                      tgtAGTTCGAACATAATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGATTCCCATGTCCGCTCTGTTCGAACCCTACACCCTGAAAGACGTCACCCTCCGTAACCGTATCGCCATTCCGCCCCTGTGCCAGTACATGGCCTAAGACTGCa         >  4:28354/3‑149 (MQ=255)                                                                                                                                             cGATTATGATCGCACAACAGTATTCAT‑A‑TGAACCCGGTGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGCCGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATg                      <  4:77222/129‑1 (MQ=255)                                                                                                                                             cGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGAGCCAGTACATg                      >  3:77222/1‑129 (MQ=255)                                                                                                                                              gATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCaa  <  6:86519/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CCAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCAA  >  NC_002947/1431736‑1432024

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