Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,167,515	C→G	7.5%	intergenic (‑40/+122)	creA ← / ← katG	creatinase/catalase‑peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,167,515	0	C	G	7.5%	165.5 / 4.5	53	intergenic (‑40/+122)	creA/katG	creatinase/catalase‑peroxidase
Reads supporting (aligned to +/- strand):  ref base C (21/28);  new base G (2/2);  total (23/30)
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.21e-01

CATATAGGCCCGCAGCCGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGC  >  NC_002947/4167371‑4167661                                                                                                                                                 |                                                                                                                                                    cATATAGGCCCGCAGCCGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAggccggcc                                                                                                                                                 <  6:331234/149‑1 (MQ=255)          ccGCAGCCGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGTGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTgag                                                                                                                                          >  4:153529/1‑147 (MQ=255)          ccGCAGCCGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCGGCCGGCCTgaaca                                                                                                                                        <  3:52823/149‑4 (MQ=255)          ccGCAGCCGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTgag                                                                                                                                          >  4:188451/1‑147 (MQ=255)            cgggcctatatggccgagcagCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGaa                                                                                                                                        <  8:148327‑M2/126‑1 (MQ=255)               gCCGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCAtt                                                                                                                                   <  3:266005/149‑1 (MQ=255)                ccGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTg                                                                                                                                  <  5:232458/149‑1 (MQ=255)                   ggCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCaa                                                                                                                               <  3:312802/149‑1 (MQ=255)                   ggCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCaa                                                                                                                               <  7:217504/149‑1 (MQ=255)                          cgagcagCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCGGCCggcctgaacagc                                                                                                                                      <  8:269884‑M2/128‑6 (MQ=255)                          cgagcagCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCGGCCggcctgaacagc                                                                                                                                      >  7:269884‑M2/8‑130 (MQ=255)                           gagcagCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAAcc                                                                                                                        >  1:141778‑M2/7‑148 (MQ=255)                               agCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACcgacga                                                                                                                   <  5:255636‑M2/147‑1 (MQ=255)                                        cTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTgag                                                                                                                                          <  1:140507/117‑1 (MQ=255)                                        cTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTgag                                                                                                                                          >  2:140507/1‑117 (MQ=255)                                        cTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACttt                                                                                                          >  5:148175/1‑149 (MQ=255)                                         tGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAg                                                                                                                  >  7:49334/1‑140 (MQ=255)                                         tGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAg                                                                                                                  <  8:49334/140‑1 (MQ=255)                                          gAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTg                                                                                                        >  2:16233/1‑149 (MQ=255)                                           aGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGTGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTg                                                                                                        <  3:153529/148‑1 (MQ=255)                                                aaaaGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAAc                                                                                                                         >  2:290434/1‑126 (MQ=255)                                                aaaaGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAAc                                                                                                                         <  1:290434/126‑1 (MQ=255)                                                 aaaTGTTGGCTGCACCTTGTCGCCATTGCTGATTTTCAGGGTCTTGTGCATTTTCATGGTCGGGTTTCCTTGGCGTTGTTAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGtatctat                                                                                                 <  5:266997/149‑1 (MQ=255)                                                    ggTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACcga                                                                                                                      <  3:36446/125‑1 (MQ=255)                                                    ggTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACcga                                                                                                                      >  4:36446/1‑125 (MQ=255)                                                           tGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTgg                                                                                        >  3:297020/1‑148 (MQ=255)                                                           tGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTgg                                                                                        >  1:66253/1‑148 (MQ=255)                                                           tGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTgg                                                                                        >  3:89653/1‑148 (MQ=255)                                                             cACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTgg                                                                                        <  4:297020/146‑1 (MQ=255)                                                                         aTTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAAc                                                                         <  2:141778/149‑1 (MQ=255)                                                                              ggATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACcgac                                                                                                                     <  1:226645/100‑1 (MQ=255)                                                                              ggATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACcgac                                                                                                                     >  2:226645/1‑100 (MQ=255)                                                                                              gggCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGgcgc                                                       <  6:148175/146‑1 (MQ=255)                                                                                               ggCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAg                                                   >  8:93213/1‑149 (MQ=255)                                                                                               ggCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACATGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTgcg                                                   >  4:240766/1‑147 (MQ=255)                                                                                                     tGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAgcgcg‑c                                             <  3:188451/148‑1 (MQ=255)                                                                                                                ggTTCCTTGTCGTTGTGAGCACGCTGTTCAGGTCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAg                                 >  5:230767/1‑149 (MQ=255)                                                                                                                ggTGCCTTGTCGTTGTGAGCCCGCTGTTCAGGCCCGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAg                                 <  5:346189/149‑1 (MQ=255)                                                                                                                      ttGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAgc                                                  <  2:228698/127‑1 (MQ=255)                                                                                                                      ttGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCA‑CCCCCCCCGCCAg                                 <  2:337654/143‑1 (MQ=255)                                                                                                                      ttGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCA‑CCCCCCCCGCCAg                                 >  1:337654/1‑143 (MQ=255)                                                                                                                      ttGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTCAGCGCG‑CCCCCCCCGCCAGGCAACg                           >  8:27364/1‑149 (MQ=255)                                                                                                                           gTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGCCCCCCCCCGCCAGGCAACGCTAt                       >  5:11180/1‑149 (MQ=255)                                                                                                                            ttGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCCGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTAttt                     <  1:16233/149‑1 (MQ=255)                                                                                                                                 accACGCTGATCCGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGa                 <  7:298672/146‑1 (MQ=255)                                                                                                                                 aGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGa                 <  4:89653/148‑1 (MQ=255)                                                                                                                                 aGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGa                 <  2:66253/148‑1 (MQ=255)                                                                                                                                      gCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAg           <  3:240766/149‑1 (MQ=255)                                                                                                                                      gCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAg           <  7:93213/149‑1 (MQ=255)                                                                                                                                      gCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAg           >  8:234596/1‑149 (MQ=255)                                                                                                                                      gCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAg           >  6:89426/1‑149 (MQ=255)                                                                                                                                        tGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGGG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCg         <  7:27364/149‑1 (MQ=255)                                                                                                                                         tctCAGGCGGCCGGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCgg        >  8:212843/3‑149 (MQ=255)                                                                                                                                          ttCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCgg        <  7:234596/148‑1 (MQ=255)                                                                                                                                              ggccggccGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAg   >  7:161796/1‑149 (MQ=255)                                                                                                                                               gccggccgCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAg   >  2:87723/1‑148 (MQ=255)                                                                                                                                               gccggccgCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGc  >  5:316279/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CATATAGGCCCGCAGCCGGGCGTGGCGAGCAGCGTATTCCTGAGCGGAAAAGGTTGGCTGCACCTTGTCGCCATTGCGGATTTTCAGGGTCTTGGGCATTTGCATGGTCGGGGTTCCTTGTCGTTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGC  >  NC_002947/4167371‑4167661

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