Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	3,845,703:1	+G	100%	intergenic (‑70/+28)	PP_3394 ← / ← PP_3395	putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	3,845,703	1	.	G	100.0%	146.2 / NA	39	intergenic (‑70/+28)	PP_3394/PP_3395	putative 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/Transcriptional regulator, LysR family
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (27/12);  total (27/12)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCT  >  AE015451/3845477‑3845906                                                                                                                                                                                                                                    |                                                                                                                                                                                                            gcaAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATc                                                                                                                                                                                                     >  1:688482/1‑236 (MQ=255)                             ttCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCcg                                                                                                                                                                          <  1:732693/235‑1 (MQ=255)                                            gCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTGAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGc                                                                                                                                                             >  2:407191/1‑233 (MQ=255)                                            gCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGccc                                                                                                                                                                           >  2:375027/1‑219 (MQ=255)                                            gCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                                                                                           >  1:386810/1‑235 (MQ=255)                                             ccGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATcc                                                                                                                                                                                                    >  1:65828/1‑193 (MQ=255)                                              cGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                                                                                               <  1:643927/197‑1 (MQ=255)                                                                                ggCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgag                                                                                                                       <  2:452568/235‑1 (MQ=255)                                                                                        tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTACCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                                                                                               >  1:210953/1‑235 (MQ=255)                                                                                          ggCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTg                                                                                                                                                        >  2:508641/1‑192 (MQ=255)                                                                                                  cACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCAGCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGAGCAGAAAAGACGTGTTTCGCGGGCTCATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggcgg                                                                                                                                       >  2:721900/1‑201 (MQ=255)                                                                                                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCa                                                                                                                                                                                         <  2:155857/140‑1 (MQ=255)                                                                                                             tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                                                                                              >  2:160453/1‑183 (MQ=255)                                                                                                                gCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                                                                                       >  1:537326/1‑235 (MQ=255)                                                                                                                  aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGcc                                                                                                                                                                                          >  1:590926/1‑134 (MQ=255)                                                                                                                          gcgcAGGCCGACTTCATTGATTGTTATTGTTGTCACCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                                                                                            <  1:613471/220‑1 (MQ=255)                                                                                                                            gcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                                                                                               >  1:316982/1‑135 (MQ=255)                                                                                                                            gcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                                                                                               <  2:316981/135‑1 (MQ=255)                                                                                                                                gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCAcc                                                                      <  1:440089/236‑1 (MQ=255)                                                                                                                                   gACTTCAGTGATTGTTATTGTGGTCATCGCGGTCGCGCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGcc                                                                                                                                                                      <  2:16452/137‑1 (MQ=255)                                                                                                                                     cTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                                                                                      >  2:416329/1‑183 (MQ=255)                                                                                                                                     cTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                                                                                      <  1:416330/183‑1 (MQ=255)                                                                                                                                           tgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                                                                                               >  2:373448/1‑184 (MQ=255)                                                                                                                                                  tattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcag                                                     >  1:321337/1‑235 (MQ=255)                                                                                                                                                          tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                                                                                             >  1:700345/1‑107 (MQ=255)                                                                                                                                                          tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCCACAGGTGGTCaa                                                                                                                                   >  2:174994/1‑149 (MQ=255)                                                                                                                                                              cgcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCAt                                                                                        >  1:142140/1‑188 (MQ=255)                                                                                                                                                               gcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCa                                                                                                                                                      >  2:195971/1‑125 (MQ=255)                                                                                                                                                               gcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCa                                             >  2:672202/1‑230 (MQ=255)                                                                                                                                                                  ctctctcTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGc                                    <  2:703355/235‑1 (MQ=255)                                                                                                                                                                    ctctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAg                                  >  1:729910/1‑236 (MQ=255)                                                                                                                                                                                  cGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCCATGACCCGCAGGTCGAACAAcgccg                    >  1:646732/1‑236 (MQ=255)                                                                                                                                                                                  cGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAAcgcc                     >  1:470033/1‑235 (MQ=255)                                                                                                                                                                                    ttGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCa                  <  1:672203/236‑1 (MQ=255)                                                                                                                                                                                        gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCCGTCAGCGCAATGACCCGCAGGTCGAACAAcgccg                    >  1:30177/1‑230 (MQ=255)                                                                                                                                                                                        gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCGATGACCCGCAGGTCGAACAAcgc                      >  1:692074/1‑228 (MQ=255)                                                                                                                                                                                        gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAg               >  1:646080/1‑235 (MQ=255)                                                                                                                                                                                                       ggATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgcg                                                                 <  2:757616/170‑1 (MQ=255)                                                                                                                                                                                                                     aaaGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCct  >  1:613000/1‑219 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                            GCAAGCGTGCCGGCACGAACGACGCGACTTCCATATGGCGCACGCCGGCGCCGTAAGCGGCGTCGATCCAGCGGCGTTTGGCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGACCCGCAGGTCGAACAACGCCGCATAGCGGGACACGGCCT  >  AE015451/3845477‑3845906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑