Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,189	A→G	44.7%	intergenic (+396/+248)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,189	0	A	G	44.7%	‑3.9 / 31.4	27	intergenic (+396/+248)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (8/7);  new base G (10/2);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

ATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCA  >  NC_000913/4295902‑4296467                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                        tttttGGAGGGCGGATATGACGTTTAGTCCTCATGCGAGGTCAACGCCGTATAGAACCCTTAACTCTTCGTATCAGGGCAACGCCAAACGACCCAATAGCCAGATTTAGCGTGATTTTGTAGGGCGGGTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                                                                                                                                                                                                        <  1:317574/288‑1 (MQ=11)                tAAGGATTTTATTCCGCTTCAGACCTCACAGCATGATGGGACTCTAAACGCGTCTTATCAGGCCACCGCCCGACAGGGCAATTGCCTGATTTAGCGTGATTTTGTAGGACGCATAAGGCGTTTACGCCGCATCCCACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                         <  1:1502966/289‑1 (MQ=11)                             ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGAGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCGTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACa                                                                                                                                                                                                                                                             <  2:967731/288‑1 (MQ=25)                                    cgatgtccaCGCCTGATGGGACCATGTAGGTGTCTTATCAGGCTTACGCCCTACAGCGTATTAGTCTGATTTCGGGTGATTTTTTAGGTCGGATGAGGCGTTTACACCGCATCCGAAATAAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                                                                                                                                       <  1:1478684/279‑1 (MQ=11)                                             cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTCGCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCATCGCCTGACGCGACGCTTACCGCGTCTTACAAGGCCTCCCCCAGACAGCGCTTTACCCTGAtt                                                                                                                                                                                                                                             >  1:1341507/1‑288 (MQ=11)                                                 gtaTGCGACGCTTCCCACGGCTTATCAGGCCTACGCCCGATGTGACACTAGCCTGCTTTAGCGTGATTTTTTGTGTCGGCTAAGGCGGTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                        <  1:376577/287‑1 (MQ=11)                                                   aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGCTCCGATACGTCGTTTATGCCGCATCCGACATCCACGCCTGCTGCTACGCTTACCGCGTCTTATCAGGCCTTCGCCAgcccacgccctcgcatgattctgtgtg                                                                                                                                                                                                                                      >  1:284788/1‑263 (MQ=11)                                                                     cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATTCGACGCTTAACGCGTCTTATCAAGCCTACCCCCGACTGCGTAACAGCCTGCTTTAGCGGTATTTTCTAGTTccggttc                                                                                                                                                                                                                    >  1:1400711/1‑283 (MQ=11)                                                                                       gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTGCTGCGACGCTTAACGCGTCTTATCAGGCCTCCCCCAGACATCGCAACAGCCTGATTTCGCGTGATTTTCTCGATCGGAtccggagattccgccgcagc                                                                                                                                                                                                   >  1:388709/1‑269 (MQ=11)                                                                                       gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTATGCCAAACAGCGCAATAGCCCGATTTCGCGTGATTCTGTAGTTCGTATAATGCGCTTCAGCCCCATcc                                                                                                                                                                                                  >  1:383354/1‑289 (MQ=11)                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                                                                                                                                                                                                        >  1:1252199/1‑198 (MQ=34)                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                                                                                                                                                                                                        <  2:1252199/198‑1 (MQ=34)                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCTGATTCGACGCTTATCGCGTCTTATCAGGCCTACGCACGACAGTGCACTATCCTGCTTTAGCTTGACTCTGTAGGTCGTCTAAGGcttgtacgctccattccact                                                                                                                                                                                              >  1:1278102/1‑269 (MQ=11)                                                                                                        aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                         <  1:1422193/201‑1 (MQ=34)                                                                                                                  aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                                                                                               >  2:1422193/1‑201 (MQ=34)                                                                                                                   ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                                                                                                                                                                                                           <  2:1383283/236‑1 (MQ=33)                                                                                                                             cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                                                                                                                                                                                 >  1:1383283/1‑236 (MQ=33)                                                                                                                                                              cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                                                                                                                                       <  2:933688/165‑1 (MQ=11)                                                                                                                                                                                                                        gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATACGTCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCGTGCGACCTCAACCCCTGATTCGACACTTACCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATCGCCTGAGTCAGAGTTGTATTGTAGGTCGGATCAGTCACTTACGACGCCGCCTACA‑CAATGTCCTGAt                                                                >  1:440122/1‑289 (MQ=11)                                                                                                                                                                                                                                      ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATT‑AGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaatgcctgatgcgacgcttgccgcg                                                  >  2:542616‑M2/1‑264 (MQ=255)                                                                                                                                                                                                                                              cGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaatgcctgatgcgacgcttgccgcgtcttat                                            >  2:70181‑M2/1‑257 (MQ=255)                                                                                                                                                                                                                                                  gagtgggtttacgCCGCATCCGACATCAACGCCGGATGCGACGCTTGCCGCGTCTGATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGACATCaatgcctgatgcgacgcttgccgcgtcttatcaggc                                       >  1:1364933‑M2/7‑253 (MQ=255)                                                                                                                                                                                                                                                          ttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTa                              >  2:1139354/1‑289 (MQ=17)                                                                                                                                                                                                                                                           tacgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CATGATGCGACGCTTGCCGCCCCTTATCAGGCCTATCTTaa                             >  1:1062364/4‑289 (MQ=255)                                                                                                                                                                                                                                                           tacgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTa                              >  2:695948/4‑288 (MQ=17)                                                                                                                                                                                                                                                                cGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGAATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGACTATCTTAACCGtt                        >  2:678224/1‑289 (MQ=255)                                                                                                                                                                                                                                                                 gCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                       >  2:666933/1‑289 (MQ=37)                                                                                                                                                                                                                                                                               cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa         <  2:86127/289‑1 (MQ=255)                                                                                                                                                                                                                                                                               cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCTAGCCAGCGCGATAGCCTGATTTAGCGTGATTTTGTAGCTCGGATAAGGCGTTTACGCGGCATCCGACATCAATG‑CCTGATGCGTCGCTTGCCGCGTCTTACAAGGCCTATCCTAACCATCGGGTAATTTTTCaa          >  1:91959/1‑288 (MQ=11)                                                                                                                                                                                                                                                                                       gCGACGCTTGCCGCGTCTTATCAGGCCTACGCCGGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca  <  2:284788/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                                aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                      <  1:1130940/148‑1 (MQ=14)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                        ATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCA  >  NC_000913/4295902‑4296467

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