Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,428	G→A	intergenic (+634/+8)	gltP → / ← yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,428	0	G	A	100.0%	81.8 / NA	31	intergenic (+634/+8)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/18);  total (13/18)

TATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTTGGCG  >  NC_000913/4296199‑4296670                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTTATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTgg                                                                                                                                                                                                                                  <  2:249818/250‑1 (MQ=255)  ggcatgccttcgccaatccgagcAAGAGCCTGTGTTATCGTGTTTTTGTGGGTCGGATAAGGCGTTTAGGCTGCATCCGACACCAACGTCTGATGCGACGCTTAGCGCGTCTTATCAGTCCAACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTGGCCGCGTCTTAGCAAGCCTATCTTAACCGTTGGt                                                                                                                                                                                                                                 <  2:190284/229‑1 (MQ=11)      ggCCTACGCCAGGCAGCGCAATAGCCTGATTTAGCGTGATTTTGTTGGGCGGATAAGGCGTTTTCTCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCGGATGCGTCGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAAtt                                                                                                                                                                                                                            <  2:19056/251‑1 (MQ=11)           aCGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTACCCGTTGGTTAATTTTTCa                                                                                                                                                                                                                       <  1:49084/251‑1 (MQ=255)                  aCAGCACAATAGCCTGATTGAGCGTGATTTGGTCGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACTC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTGAAATTTTCAAACTCtt                                                                                                                                                                                                                <  2:133274/251‑1 (MQ=14)                              ccTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACACCAGCCAGCGCAATAGACTGATTTATCTTTATTTCGTACATCTGATCCGGCGTTAACGC‑‑CGCATCCGACCTCATTGCGTGCTTCTCCGTTGGCCGCGTCTTTTCAAGCCTATCTTCACCGTTGGTTAATCTTTCAAACTCATCACACCCGGTa                                                                                                                                                                                                     >  2:228333/1‑250 (MQ=2)                                                tGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCCAGc                                                                                                                                                                                  >  2:49562/1‑251 (MQ=255)                                                            aaggagtttacgcctcagcctccATCAACGCCTGATGCGACGCTTAACGCGTTTTTTCAGTCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCTCGTCTTATCAAGCCTATCTTAACCGGTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCGATACATCTCAGGTTCAAc                                                                                                                                                                       <  2:171194/228‑1 (MQ=17)                                                            aaGGCGTTTATGACGCATCCGCCATCAACACCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAACCGGTATCAAACCCTTCCATACAGCTCAGGTTAAAcc                                                                                                                                                                      >  1:196282/1‑251 (MQ=18)                                                             aGGCGTTTATGCCGAATCCGACATCAACGCCATATGCGACGCTTAACGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTTGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa                                                                                                                                                                     <  1:144926/251‑1 (MQ=18)                                                                                   aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa                                                                                                                                                                        >  1:208977/1‑226 (MQ=255)                                                                                   aTCAACGCCTGATGCGACGCTTAACGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGGTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa                                                                                                                                                                        <  2:208977/226‑1 (MQ=18)                                                                                        cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                                                                                                          >  2:199511/1‑251 (MQ=255)                                                                                             gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATa                                                                                                                                      <  1:187038/250‑1 (MQ=255)                                                                                               tGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaaa                                                                                                                                    >  1:211553/1‑250 (MQ=255)                                                                                               tGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCCCACCCGGTATCAACCCCTTCCATACAGCTCAGGTTCCCCCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaaa                                                                                                                                    >  2:192923/1‑250 (MQ=21)                                                                                                 cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAcc                                                                                                                                 <  2:185661/251‑1 (MQ=255)                                                                                                   aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACctt                                                                                                                               <  1:148577/251‑1 (MQ=255)                                                                                                                                 aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                                                                                 <  2:143884/251‑1 (MQ=255)                                                                                                                                  gACAGAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGACTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                                                                                 >  1:262120/1‑250 (MQ=255)                                                                                                                                     aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGt                                                                                              >  2:115098/1‑250 (MQ=255)                                                                                                                                         cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CGACGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCg                                                                                         >  2:100435/1‑251 (MQ=255)                                                                                                                                           aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGa                                                                                       <  1:199083/251‑1 (MQ=255)                                                                                                                                              gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAAt                                                                                     >  1:58529/1‑250 (MQ=255)                                                                                                                                                             attgtGTAGGTCGGTTAAGGCGTTTACGC‑‑CGCAGCCGGCATTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGTGAAAt                                                                     <  2:227742/247‑1 (MQ=255)                                                                                                                                                                      gTCGGATAAGGCGTTTACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCt                                                            <  1:56996/251‑1 (MQ=255)                                                                                                                                                                                                      aTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCAtc                              <  2:2584/251‑1 (MQ=255)                                                                                                                                                                                                          aTGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCAtcgtcg                          >  1:18312/1‑251 (MQ=255)                                                                                                                                                                                                                         cTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACg           >  1:209028/1‑251 (MQ=255)                                                                                                                                                                                                                           tGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACg           <  2:126934/249‑1 (MQ=255)                                                                                                                                                                                                                                     tATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTTGGCg  <  1:13004/248‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   TATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTTGGCG  >  NC_000913/4296199‑4296670

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