Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	331,879	G→A	intergenic (+3/‑492)	yahV → / → pdeL	protein YahV/DNA‑binding transcriptional regulator/c‑di‑GMP phosphodiesterase PdeL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	331,879	0	G	A	100.0%	87.7 / NA	27	intergenic (+3/‑492)	yahV/pdeL	protein YahV/DNA‑binding transcriptional regulator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/12);  total (15/12)

CAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCAAAAAT  >  NC_000913/331643‑332123                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     cAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAAt                                                                                                                                                                                                                                        >  2:6735/1‑251 (MQ=255)    aaCCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATg                                                                                                                                                                                                                                       >  1:160770/1‑249 (MQ=255)                gggtgtttgggaaggtCAATAAAGTTAATCTTATTTTAAGTTAATAATGCTGAACAATTTGAA‑TGACGGATTTATCAAGGGAAAAGTGTAACGTAAAACAACTTAAAATTAAGGCA‑CCATTTTGTATTTTAAGGTT‑AATTTACTGTGTTTGAAAATTTATTTGAAT‑TTTTAAATATTGGTTTTTAT‑GGGATTGCTTTAA‑ATATGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAg                                                                                                                                                                                                                         <  1:385882/236‑1 (MQ=255)                          aaTGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCAt                                                                                                                                                                                                               >  1:489660/1‑251 (MQ=255)                                cAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGt                                                                                                                                                                                                          <  2:528377/250‑1 (MQ=255)                                cAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTa                                                                                                                                                                                                         <  2:160770/251‑1 (MQ=255)                                                    aaTGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GCAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAAcc                                                                                                                                                                                      >  1:68702/1‑250 (MQ=255)                                                             atTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGGTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCttt                                                                                                                                                                             <  2:68702/250‑1 (MQ=255)                                                                  tGATCC‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAACGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCAt                                                                                                                                                                       >  2:276025/1‑251 (MQ=255)                                                                       a‑TTTGAAATGACGCATTATTCATGAGCAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGt                                                                                                                                                                  >  2:330292/1‑251 (MQ=255)                                                                                tGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GACATTAACGCCACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGGGATATTTTATT‑GAATGTTTTAAGTATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTc                                                                                                                                                                         >  2:89918/1‑236 (MQ=255)                                                                                tGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTAACTGTTTGTGATATTTTATT‑TAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCGTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAAGGATACGTATGCATAACCAACATCTTTACTc                                                                                                                                                                         <  1:89918/236‑1 (MQ=255)                                                                                   cGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACGCATTATGTCATTGAATGtt                                                                                                                                                       <  2:536759/251‑1 (MQ=255)                                                                                      attattCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGAc                                                                                                                                                    >  2:327578/1‑251 (MQ=255)                                                                                                      gtgtATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGa                                                                                                                                     <  1:421453/250‑1 (MQ=255)                                                                                                        gtATCGTACATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAACGTATAATTATCTGTGTGTGATATTTAACT‑GAATGTTTTCAATATTG‑TTTTTATTGGAATTGCTCTAATAT‑TGGTTATCCTTTGCTGAATAGCTTCAGTCTTAATGAGTGGGTTTTTAAGCGACAGGAATAGAGTTATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGTggg                                                                                                                                  >  1:57658/1‑251 (MQ=255)                                                                                                            cGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGc                                                                                                                                                                                            >  2:99705/1‑189 (MQ=255)                                                                                                            cGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGc                                                                                                                                                                                            <  1:99705/189‑1 (MQ=255)                                                                                                                   cAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTAtttt                                                                                                                       >  1:515711/1‑251 (MQ=255)                                                                                                                            aTTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGAt                                                                                                               <  2:515711/251‑1 (MQ=255)                                                                                                                                                                               ttATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCa                                                              <  1:395771/250‑1 (MQ=255)                                                                                                                                                                                            ttAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTTGTggg                                                  <  2:426035/250‑1 (MQ=255)                                                                                                                                                                                                        ttttATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTCAGGGAAAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAat                                       >  1:544883/1‑250 (MQ=255)                                                                                                                                                                                                        ttttATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAat                                       <  1:327578/250‑1 (MQ=255)                                                                                                                                                                                                                                        aTTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTAcc        >  1:318962/1‑250 (MQ=255)                                                                                                                                                                                                                                        aTTTGCTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGAATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTAcc        >  1:303006/1‑250 (MQ=255)                                                                                                                                                                                                                                             cTGAATAGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCAAAAAt  >  1:398450/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     CAAAACCTACATTTGGGCTGTTGCGAATGTTCAATAAGTTTAGTCTTATTTAATGTAAATATTGCTGATCA‑TTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACT‑GAAATTAACGCAACCATTT‑GTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATT‑GAATGTTTTAAATATTG‑TTTTTATTGGCATTGCTATAATAT‑TGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCTTTGCTCATGAATGGATGCCAGTATGTAGTGGGAAATTATAAATATTGAAATAGTCCAACTACTTCTTTATTACCAAAAAT  >  NC_000913/331643‑332123

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