Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,177	G→A	intergenic (+2192/+551)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,177	0	G	A	100.0%	87.5 / NA	23	intergenic (+2192/+551)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/12);  total (11/12)

AGTTTGGTGTTGATAACGGTTACCTGGCAGATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACA‑GAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAA‑GTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTT  >  NC_000913/1651939‑1652423                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             aGTTTGGTGTTGATAACGGTTACCTGGCAGATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAAt                                                                                                                                                                                                                                                  >  2:91263/1‑251 (MQ=255)         gtTGATAACGGGTACCTGGCAGATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCaa                                                                                                                                                                                                                                           <  2:182810/250‑1 (MQ=255)               aaCGGTTACCTGGCAGATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGAcc                                                                                                                                                                                                                                     <  1:56815/250‑1 (MQ=255)                  ggTTACCTGGCAGATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGttt                                                                                                                                                                                                                                 >  1:160072/1‑251 (MQ=255)                                          aaGGGGATCTCGGCATTACAGGAAACAAGAACCATTCCGGCTCCTCTTTCGAGGGAAGAAGTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAACATTTATGGTGTGTTTCTGTTTCTACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTTGATGGGAGGACATAGATCTGAAAAATGGAACAATGATACTCAGGAGAAATTTATCAAAAGACCGTTTCACGGTACCAAAAACACAGGCggg                                                                                                                                                                                                         >  2:22171/1‑251 (MQ=255)                                                                              ccGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCaa                                                                                                                                                                     >  1:112360/1‑250 (MQ=255)                                                                                       ctTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCTTCTTATTAATCCAGCAATCGACCctcc                                                                                                                                                           >  2:365222/1‑250 (MQ=255)                                                                                         ttCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGctc                                                                                                                                                            <  1:91263/248‑1 (MQ=255)                                                                                                           ccGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTa                                                                                                                                       <  1:184711/250‑1 (MQ=255)                                                                                                                          gTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGc                                                                                                                           <  2:333356/247‑1 (MQ=255)                                                                                                                          gTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AAACAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGc                                                                                                                           >  1:333356/1‑247 (MQ=255)                                                                                                                             cAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGCAAAGAG‑c                                                                                                                   <  1:173358/251‑1 (MQ=255)                                                                                                                                         agcaaAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATg                                                                                                       <  1:177840/251‑1 (MQ=255)                                                                                                                                                                  gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCCCTTCAGAGAGGATGGcagaa                                                                               >  1:340458/1‑250 (MQ=255)                                                                                                                                                                  gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGcag                                                                                 <  1:77585/248‑1 (MQ=255)                                                                                                                                                                          tGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGAACAAGAGGGGACAT‑CATTTTTTTAAGCCCGCAATCGTCGCCCTCCCGC‑GGCAGAAGACATTAAAGAGACGTAACCAAGAG‑CATATCACTGGTTTTCTCTTCAGAGAGTATGGCAGAACAAGaaaa                                                                       >  2:100244/1‑250 (MQ=255)                                                                                                                                                                           ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AAACAAACGGGTCATTCATCTTATTAAGCCCTCAATCCACGCCCTCCGGG‑GTGCGAAGGAATTAACGGGACTGAGAAAAGAG‑GATATTATTAATGTTCATTTCAGAGCGAAAGGCAGAAAA‑GAAAAAc                                                                     >  1:394613/1‑250 (MQ=255)                                                                                                                                                                            gAGTTCGCCCTGGTTAGCTGTGTGCACTTGGATGGGAGGACATAGATCGGAAAAATGGAACAATGAGAAACAGGAGAAATTTAGCAAAAGACCGTGTCACGGTACCAAAAACACAGGCGGT‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGAAAAGAG‑CATATCATTGATGTTCACTTCGGAGAGTATGGCAGAACA‑GAAAAACaa                                                                   <  2:169762/251‑1 (MQ=255)                                                                                                                                                                                               gtgggCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCTAAAGACCGTTGCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGGCGCTCTCCGGA‑GTCAGATGACATGAACGAGACTGAGCAAATAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCATAACA‑GAAAAGCAAAAT‑TGCACCTtttttttt                                                <  2:43092/246‑1 (MQ=255)                                                                                                                                                                                                     aCTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCACTCGGGGCATTCATCTTATTAAACCAGCAATCGACGCTCTCCGGA‑GGCAGATGACATTAACGCGACTGAGCAAAGAG‑CCAATCATTGATGTTCACTTCAGAGAGTATGGCAGAACA‑GAAAAACAAAAA‑TGCACCTTTGTTTTTTAACCTg                                         >  1:80314/1‑251 (MQ=255)                                                                                                                                                                                                                  gaCATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAACAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑ACCCAATCGGGGCCATCATCCTTTTAAGCCCGCCATCGACGGCCTCCGGA‑GTCCGAGGACATTAACCAGACT‑AGAAAAAAGGCATATCATTGATGGTCACATCAGAGAGTATGGCAGAACA‑AAAAAAAAAAAA‑TACACCTTTGATTTTCAACCTTAA‑gttgcagcg                             >  1:275247/1‑242 (MQ=255)                                                                                                                                                                                                                  gaCATAGATCTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACGTTAACGAGACTGAGCAAAGAG‑Cat                                                                                                                 <  2:180330/168‑1 (MQ=255)                                                                                                                                                                                                                           cTGAAAAATGGAACAATGATAATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAAGCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGcagaaca‑g                                                                           <  1:304012/196‑1 (MQ=255)                                                                                                                                                                                                                                           tgataaTCAGGAGACATTTAGCAAAAGCCCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCAGGTCCTTCATCTTATTTAGACAGCAAACCACGATCTCCGGAAGT‑AGATGACATTAACAAGACTGAAGAAAGAG‑CATATCCTTGATTTTCACATCTGAGAGTATGGCCGAACA‑GAAAAAAAAAAAATG‑ACCTTTTTTTTTCAACCTGCA‑GGGTCAGCGAGAGTAAAAAATTATTGTGCACAtttt  >  2:66163/1‑251 (MQ=255)                                                                                                                                                                                                                                           tgataaTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACCCAGGCGGG‑AAACAAACGGGGCAGTCATCTTATTAATCCAGCAATCGGCGCTCTCCGGA‑GTCAAAACACCTTTACGAGACTGGACAAAAAG‑CATATCATTGATGTTCAATTTAGAGGGAATGGCAGAACA‑GAAAAACAAAAA‑TGCACCTTTTTTTTTCAAACTTAACGGGTC‑GCGAGAGTAACAAAATAGGGTGGCCAtttt  >  1:421161/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                             AGTTTGGTGTTGATAACGGTTACCTGGCAGATAATCCGTTTAAGGGGATCTCGCCATTAAAGGAATCAAGAACCATTCCGGATCCTCTTTCGCGGGAAGAATTTATCCGTCTTATCGATGCGTGCAGAAATCAGCAAGCAAAAAATTTATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGG‑AACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGA‑GTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACA‑GAAAAACAAAAA‑TGCACCTTTGTTTTTCAACCTGAA‑GTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTT  >  NC_000913/1651939‑1652423

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