Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,169,591	(C)8→9	intergenic (‑150/+12)	creA ← / ← katG	creatinase/catalase‑peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,169,583	1	.	C	100.0%	111.6 / NA	31	intergenic (‑142/+20)	creA/katG	creatinase/catalase‑peroxidase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/17);  total (14/17)

TTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCACTGCTGGCATACACCTCGCT  >  NC_002947/4169460‑4169704                                                                                                                             |                                                                                                                           tggtgAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGAATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑‑cccccccccgcc                                                                                                               <  1:382886/134‑5 (MQ=255)  tgtgAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑‑cccccccccgcca                                                                                                              <  2:467184/136‑6 (MQ=255)  tgtgAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑‑cccccccccgcca                                                                                                              <  1:366399/136‑6 (MQ=255)  tgtcAGCACGCTGTGCAGGCCGGCCGCCTGAGAAGCATTGCAACACAACCGACGAGCACACTTTTGTAACTATCCCACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑‑cccccccccgcca                                                                                                              <  2:70220/136‑6 (MQ=255)    ggAGCACGCGGTTCAGGCCGGCCGCCGGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAgg                                                                                                             <  2:94207/135‑1 (MQ=255)       gCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGcca                                                                                                          >  1:444796/1‑134 (MQ=255)       ccACGCTTTTCCGCCCGGCCGCCTGATAAGCTTTCCAACAGACCCGACGAGGACACTTTTTTATCTATCCGACCTGCAATTGGGGGGATAACCACTGCCCACGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCaa                                                                                                          <  2:1291282/135‑1 (MQ=255)                     ccggccgCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATc                                                                                            >  2:325547/1‑136 (MQ=255)                            ccTGATATGCATTGCAACAGAACGGACGAGGACAGTTTTTTATATATCCGACTTGGATTTGGCGGTATAACCACTGGCCAAGGGGGCGCTGCGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCgg                                                                                     <  2:487332/136‑1 (MQ=255)                            ccTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCa                                                                                           >  1:1026788/1‑130 (MQ=255)                            ccTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCa                                                                                           <  2:1026788/130‑1 (MQ=255)                            ccTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCgg                                                                                     >  2:111280/1‑136 (MQ=255)                                         gCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTc                                                                                   >  2:872667/1‑125 (MQ=255)                                         gCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTc                                                                                   <  1:872667/125‑1 (MQ=255)                                          cAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGAc                                                                       >  2:678809/1‑136 (MQ=255)                                                        aGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGCCCCCCCCCCGCCCGGCAACGCTTTTTCAGATCACAGCGGTCCAGCTCCCAGACATAATGCCAcgcg                                                          >  1:1037552/1‑136 (MQ=255)                                                                        tatcCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGAc                                         <  2:1208033/136‑1 (MQ=255)                                                                        aatcCGACCTGGAATTGGCGGGAAAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTACTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGAc                                         <  2:1102366/135‑1 (MQ=255)                                                                            cGACCTGGAATAGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTAAGATAAAAGCCGTCCAGCACCATGACCTTCTGCCAGGCGCCCACAAATTCTAGGACaaaa                                     >  1:779509/1‑135 (MQ=255)                                                                                tggtatttGGCGGTTTAACCCGTGCCCACGGGGGCGCTCCGGGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCAGCCAGGCGGCCACAAAGTCTCGGACAAACTTGt                                 <  1:79587/130‑1 (MQ=255)                                                                                      ttGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGc                           <  2:727323/136‑1 (MQ=255)                                                                                           ggtaTAACCACTGGCCTTGGGGGCGCTCTGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCTAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCa                      <  1:1003095/133‑1 (MQ=255)                                                                                           gggATAACCACTGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCa                      <  1:111280/136‑1 (MQ=255)                                                                                           gggATAACCACGGGCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCa                      <  1:1366573/136‑1 (MQ=255)                                                                                           gggAAAACCACTGGCCACGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCa                      <  1:1271593/136‑1 (MQ=255)                                                                                           ggcaaaaCCCCTGCCCAGGGGGGCGGGGCCCGCCC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCa                      <  1:1088071/131‑1 (MQ=255)                                                                                             caTAACCACTGCCCAAGGGGGCGCGGCGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCAct                    <  2:1305692/135‑1 (MQ=255)                                                                                                       ggCCAAGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCACTGCTGGCATac          <  2:1281327/136‑1 (MQ=255)                                                                                                           aaGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCCCCGCCATCACTGCTGGCATACAcca      >  1:743036/1‑135 (MQ=255)                                                                                                           aaGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCACTGCTGGCATACAccg      >  1:24889/1‑135 (MQ=255)                                                                                                           aaGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCACTGCTGGCATACACCt      >  2:1122065/1‑136 (MQ=255)                                                                                                           aaGGGGGCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCAAGCCCTTCTGCCAGGCGGCCACAAAGTCTCGGCCAAACATGTCACCCCCATCACTGCTGGCATACAcca      >  1:784384/1‑135 (MQ=255)                                                                                                           aaGGGGGCGCTCAGCGCGC‑CCCCCCCCGCCAGGCAACGCAATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTAGTCACCGCCATCACAGCTGGCATACAcca      >  1:546942/1‑135 (MQ=255)                                                                                                              ggggCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCACTGCTGGCATACACCTCGc   <  2:747565/136‑1 (MQ=255)                                                                                                               gggCGCTGAGCGCGC‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCACTGCTGGCATACACCTCGcg  >  1:1364665/1‑135 (MQ=255)                                                                                                                             |                                                                                                                           TTGTGAGCACGCTGTTCAGGCCGGCCGCCTGAGAAGCATTGCAACAGAACCGACGAGGACACTTTTGTATCTATCCGACCTGGAATTGGCGGGATAACCACTGGCCAAGGGGGCGCTGAGCGCG‑‑CCCCCCCCGCCAGGCAACGCTATTTCAGATCAAAGCGGTCCAGCTCCATGACCTTCTGCCAGGCGGCCACAAAGTCTCGGACAAACTTGTCACCGCCATCACTGCTGGCATACACCTCGCT  >  NC_002947/4169460‑4169704

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