Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	353,067	A→G	intergenic (+34/+84)	hisF → / ← cbcV	Imidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	353,067	0	A	G	100.0%	86.3 / NA	25	intergenic (+34/+84)	hisF/cbcV	Imidazole glycerol phosphate synthase subunit HisF/choline / betaine / carnitine ABC transporter ‑ ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/14);  total (11/14)

CTGACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTT  >  AE015451/352788‑353332                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                          ctgACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTC‑GGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                                                                                                                                                                  >  2:563934/1‑288 (MQ=255) ctgACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTC‑GGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                                                                                                                                                                  <  1:563955/288‑1 (MQ=255)    aCCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                                                                                                                                                    <  1:389646/300‑1 (MQ=255)       aGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                                                                                                                                                     >  2:911069/1‑280 (MQ=255)            ggACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATGGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGt                                                                                                                                                                                                                                                      <  1:911095/290‑1 (MQ=255)              aCCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATAAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtc                                                                                                                                                                                                                                           <  1:523228/299‑1 (MQ=255)                  ggACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCGATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                                                                                                                                                           <  2:371314/279‑1 (MQ=255)                                aCGGCTTCGACCTGGGCGTTACCCGGGCTACAAGTGAGGCGCTGGGGAGTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGt                                                                                                                                                                                                                          <  2:536081/298‑1 (MQ=255)                                                                           ggggATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTTCAGCCCTTCGCGGGc                                                                                                                                                                                                                                                                        >  1:768440/1‑209 (MQ=255)                                                                                           gCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCg                                                                                                                                                                                                                                                                 <  2:600592/200‑1 (MQ=255)                                                                                                                           ggCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCcca                                                                                                                             <  2:1406/300‑1 (MQ=255)                                                                                                                             cAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCccacc                                                                                                                           >  2:366034/1‑300 (MQ=255)                                                                                                                                               ggCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTACTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCCACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGGATATCCCCACCACCCTGTCGTTATCCt                                                                                                           >  1:494821/1‑298 (MQ=255)                                                                                                                                                                    cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTTCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                                                                                                                                                               >  1:735091/1‑225 (MQ=255)                                                                                                                                                                                 ccACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagcagc                                                                                                                                                               <  2:735065/212‑1 (MQ=255)                                                                                                                                                                                                      gCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGc                                                     <  1:366046/297‑1 (MQ=255)                                                                                                                                                                                                        cGGAGGCCAAGGCCTACATGCCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGATTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCCAACGGCCGCCAAGCGGCCCAAATTGCTTCAACCGTGGTTCTTGcc                                                                                                                                                                      >  1:859049/1‑182 (MQ=38)                                                                                                                                                                                                                   gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAg                                                                   <  1:903945/270‑1 (MQ=255)                                                                                                                                                                                                                                     gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTcc                  >  1:784997/1‑301 (MQ=255)                                                                                                                                                                                                                                        aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTcc                  >  1:519002/1‑298 (MQ=255)                                                                                                                                                                                                                                          cgtcgtTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCAcc                                                                                                   <  1:358023/215‑1 (MQ=255)                                                                                                                                                                                                                                                      aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCtt  >  2:536749/1‑300 (MQ=255)                                                                                                                                                                                                                                                          aTTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCgg                                            <  2:19381/254‑1 (MQ=255)                                                                                                                                                                                                                                                                    gctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCATTTCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCgg                            >  2:496732/1‑260 (MQ=255)                                                                                                                                                                                                                                                                              aGCCCTTCGCGGGAACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCcagca              <  1:296514/264‑1 (MQ=255)                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                          CTGACCAGCATGGACCAGGACGGCATGAAGAACGGCTTCGACCTGGGCGTTACCCGGGCTATCAGTGATGCGCTGGGGATTCCGGTGATTGCCTCGGGGGGGGTGGGTAACCTGCAGCACCTGGCAGACGGCATTCTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCCACCACCCTGTCGTTATCCTGCAGCACCAGCTTGTTGCCGGTTTGATAACGAATCTGCAGTGCCTCACGCATGCCGATGTCGGCGTGCACCACGGTCGGCCTGCGCTCCAGCAGCTCCACATCTT  >  AE015451/352788‑353332

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