Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	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,331	0	T	C	100.0%	94.6 / NA	27	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/12);  total (15/12)

GTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652094‑1652573                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   ggttCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                                                                                                                                                                                                        <  2:14811/248‑1 (MQ=255)                          cTGGTGCGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTCATAAGCCAGCAATCGACGCTCTCCGGAGACAGATGAAATTAACGAGACTGAGCAAAGAGCATATAATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAcc                                                                                                                                                                                                              >  1:204604/1‑251 (MQ=255)                                 gCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAACGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCACTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTTCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCttttgtt                                                                                                                                                                                                       <  2:41817/251‑5 (MQ=255)                                              ggATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgt                                                                                                                                                                                          <  1:361837/251‑1 (MQ=255)                                                  gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                                                                                                                                                                                       <  1:256265/250‑1 (MQ=255)                                                   ggagcacatatgtCTGAAAAATGGAACAATGATGATCAGGAGAAAGTTAGCAAAAGACCGTTTCACGGTACCAAAAACTCAGGCGGGAACCAATCTGGTCATTCATCTTATTAAGCCAGCAATCGTCGCTATCCGTAGTCATACGGCATTAACGAGACTGAGCAAATAGCATATCATTGATGTTCACCTCTGAGAGTGTGGCAGAACAGAAAAACAAAAATGTGCCTTTGTTTTGCAACCTGAAGTGTCa                                                                                                                                                                                      <  2:494515/238‑1 (MQ=255)                                                               tCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGGCATTCAACTTATTAAGGCAGCAATCGACGCCCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGAATATCATTGAAGTTCACCTCAGAGCGGAATGCAAAAAAGAAAAAAAAAAATCCACCTTTTTTTTTCAACCTGAAACGTTAGCGAGAGTaaaa                                                                                                                                                                          >  1:261764/1‑250 (MQ=255)                                                                 tGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAt                                                                                                                                                                       >  1:268548/1‑251 (MQ=255)                                                                         ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGCCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGACCAGAAAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGTGTCAGCGAGAGGAAAAAATTATGGTg                                                                                                                                                                >  1:449173/1‑250 (MQ=255)                                                                            aCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAACCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAccc                                                                                                                                                            >  1:238145/1‑250 (MQ=255)                                                                                       cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACCAGTCATTTTTCTTATTAAACCAGCAATCGACGCTCGCCCGAGGCAGATGACATTAAAGAGACTGAGCCAAAAGCATAACATTGGTGTTTACCACAGAGAGACTGGCAGAACAGAAAAACAAAAATGCACCTTTTTTTTTCAACCTGAAGAGGCAACGAGAATAAAAAATTATGGTGTCCCTTTTTCCgatg                                                                                                                                                 >  2:33078/1‑248 (MQ=255)                                                                                                       aaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCTGGTCAGTCATCTTATTAAGCCAGCAATCGACGCTCGCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCATAGAGTTTGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGTTTACCATTTTACCTTTGTCTGTATAAGGCAGAt                                                                                                                                 <  2:256351/251‑1 (MQ=255)                                                                                                                  cACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAg                                                                                                                      <  2:401848/251‑1 (MQ=255)                                                                                                                   aCGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTAAGACAGCAATCGACGCTCTCCGGAGTCAGAGGAAATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGAACAGAAAAAAAAAAATTCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGAAAAAAATTATGGGGACCATTTTACCGTTGGCTCAATAAGGCAGATGTGGGACGCAgg                                                                                                                     >  2:175719/1‑250 (MQ=255)                                                                                                                                cacaGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGAGAAAACGTGcc                                                                                                        <  2:45877/251‑1 (MQ=255)                                                                                                                                                 cGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaa                                                                                        >  2:452905/1‑250 (MQ=255)                                                                                                                                                    gTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGCAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTCCCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGCAAATc                                                                                    >  1:380680/1‑251 (MQ=255)                                                                                                                                                          caacTTATTAAGCCAGCAAACGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATc                                                                               <  2:88733/247‑1 (MQ=255)                                                                                                                                                                    aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACACAAATGCACCTTTGTTTTTCAACCTGCAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGGGGGACGCAGCGATAAAACGTGCCGGACTACGCCATCGCAAATAATATCCGTCGAGACAt                                                                    >  2:175138/1‑251 (MQ=255)                                                                                                                                                                      ccAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAc                                                                  <  1:175719/251‑1 (MQ=255)                                                                                                                                                                        aGCAAACGACGATCTCCGGAGACAGATGACATTAACGCGACTGAGCAAAGAGCATCTCATTGATGTTCACCTCAGAGCGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTAAACCTGAAGTGTCAGCAAGAGTAAAAAATTATGGTGACCATTTTACCGTTGAATCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCTTATCAGTCGAGACATACtt                                                                >  2:226229/1‑251 (MQ=255)                                                                                                                                                                           aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                                                             <  1:182229/251‑1 (MQ=255)                                                                                                                                                                                aCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgc                                                        >  2:481407/1‑251 (MQ=255)                                                                                                                                                                                               aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCt                                         >  1:370293/1‑251 (MQ=255)                                                                                                                                                                                                                        agagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGCACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCa                 >  1:70846/1‑250 (MQ=255)                                                                                                                                                                                                                                    tGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCAGCGAAAATCATATCAGGCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATATCAAACCAGATGGGc     <  2:204604/250‑1 (MQ=255)                                                                                                                                                                                                                                      aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAt  >  1:222798/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   GTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652094‑1652573

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