Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,532,942	C→T	100%	A225V (GCC→GTC)	PP_4020 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,532,942	0	C	T	100.0%	174.2 / NA	52	A225V (GCC→GTC)	PP_4020	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/25);  total (27/25)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCCGCGCTGTACGACGCCAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGCCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCCTCTGGATGCCGACTTC  >  NC_002947/4532796‑4533087                                                                                                                                                   |                                                                                                                                                   ggCCGCGCTGTACGACGCCAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGtcc                                                                                                                                                  >  5:76106/1‑149 (MQ=255)  gCCGCGCTGTACGACGCCAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCg                                                                                                                                                 >  4:209983/1‑149 (MQ=255)     gcgcTGTACGACGCCAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ac                                                                                                                                             >  4:183097/1‑149 (MQ=255)         tGTACGACGCCAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑accgac                                                                                                                                         >  7:26868/1‑149 (MQ=255)         tGTACGACGCCAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑accgac                                                                                                                                         >  6:191008/1‑149 (MQ=255)                  ccAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCg                                                                                                                                  >  5:59670/1‑147 (MQ=255)                  ccAAAGGCGCCCACCTCAACCCGGTGAAACTGGCACAGGGCGTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGaa                                                                                                                                <  8:266237/149‑1 (MQ=255)                   cAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAg                                                                                                                               >  1:213631/1‑149 (MQ=255)                   cAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAg                                                                                                                               >  3:28315/1‑149 (MQ=255)                     aaGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGtt                                                                                                                             <  3:183097/149‑1 (MQ=255)                      aGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTc                                                                                                                            <  8:255224/149‑1 (MQ=255)                            ccACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGc                                                                                                                       >  3:305092/1‑148 (MQ=255)                              aCCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAg                                                                                                                    <  6:76106/149‑1 (MQ=255)                                   aaCCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTa                                                                                                                 <  5:191008/147‑1 (MQ=255)                                      ccGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTg                                                                                                                                     >  3:103067/1‑124 (MQ=255)                                      ccGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTg                                                                                                                                     <  4:103067/124‑1 (MQ=255)                                        gCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGc                                                                                                          <  6:40089/149‑1 (MQ=255)                                        gCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGc                                                                                                          <  6:59670/149‑1 (MQ=255)                                            aaaCTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTgg                                                                                                       <  3:179791/148‑1 (MQ=255)                                            aaaCTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTg                                                                                                        <  5:92367/147‑1 (MQ=255)                                                 ggCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCa                                                                                                 >  6:254171/1‑149 (MQ=255)                                                    acaGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACg                                                                                              >  5:215120/1‑149 (MQ=255)                                                         gCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTa                                                                                          >  4:163845/1‑148 (MQ=255)                                                               ccAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGcc                                                                                            <  3:41375/140‑1 (MQ=255)                                                               ccAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGcc                                                                                            >  4:41375/1‑140 (MQ=255)                                                                cAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACatcgatc                                                                                  <  4:131759/149‑1 (MQ=255)                                                                  gTACGATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCTGCTAGGAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCTTAGTTCGCTGCCAGTTATTGGTGCTGGCCTGCAACGc                                                                                             <  1:212540/136‑1 (MQ=255)                                                                  gTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCACGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGATCCCTGCCAGGTATTGCTGCTGGCCTGCAAAGc                                                                                             >  2:212540/1‑136 (MQ=255)                                                                      cATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCtt                                                                             >  2:14690/1‑148 (MQ=255)                                                                      cATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGACTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCtt                                                                             >  2:278223/1‑148 (MQ=255)                                                                               tGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGcc                                                                   >  6:93420/1‑149 (MQ=255)                                                                                               tACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGtt                                                   <  3:163845/149‑1 (MQ=255)                                                                                                 cGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGAACGCCTTGACCGCGGCCTGTCGCGCCGCCTGtt                                                   <  8:26868/147‑1 (MQ=255)                                                                                                  gAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGTGACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGc                                                 >  1:86523/1‑149 (MQ=255)                                                                                                  gAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGcc                                                >  6:334884/1‑149 (MQ=255)                                                                                                   aaCAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGccc                                               >  8:5140/1‑149 (MQ=255)                                                                                                          ccAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGtt                                        <  2:267050/149‑1 (MQ=255)                                                                                                           cAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTc                                       <  7:5140/149‑1 (MQ=255)                                                                                                                         ccAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGc                           <  4:28315/147‑1 (MQ=255)                                                                                                                         ccAGCAGGCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCa                         <  5:93420/149‑1 (MQ=255)                                                                                                                                gCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGcc                  >  6:72552/1‑149 (MQ=255)                                                                                                                                gCTCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGcc                  <  2:86523/149‑1 (MQ=255)                                                                                                                                  tCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGc                                                                    >  3:220681/1‑97 (MQ=255)                                                                                                                                  tCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGc                                                                    <  4:220681/97‑1 (MQ=255)                                                                                                                                  tCGCGACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCctc                >  5:210779/1‑149 (MQ=255)                                                                                                                                    gcgACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCctct               <  5:254171/148‑1 (MQ=255)                                                                                                                                    gcgACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCctct               <  2:213631/148‑1 (MQ=255)                                                                                                                                    gcgACGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCCTCTg              >  4:33424/1‑149 (MQ=255)                                                                                                                                       aCGGCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCCTCTGGAt           >  2:10724/1‑149 (MQ=255)                                                                                                                                         ggCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAg                                >  2:41996/1‑126 (MQ=255)                                                                                                                                         ggCTATGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAg                                <  1:41996/126‑1 (MQ=255)                                                                                                                                             aTGTCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCCTCTGGATGCCGAc     >  1:83702/1‑149 (MQ=255)                                                                                                                                                tCGTCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCCTCTGGATGCCGACTTc  <  1:10724/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GGCCGCGCTGTACGACGCCAAAGGTGCCCACCTCAACCCGCTGAAACTGGCACAGGGCCTGGCCAGTACCATCGAGGCTGCCGGCGGGCAAATCTACGAACAGAGCCAGGTGCTCAGCTACCAGCAGGCTCGCGACGGCTATGTCGCCCGT‑ACCGACCGTGGCGAAGTTCGCTGCCAGGTATTGGTGCTGGCCTGCAACGCCTACATCGATCGCCTTGACCGCGGCCTGTCGCGCCGCCTGTTGCCCGTCGGTTCCTACCAGGTCGCCACTGCGCCTCTGGATGCCGACTTC  >  NC_002947/4532796‑4533087

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