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%	92.6 / NA	28	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/13);  total (15/13)

GGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATG  >  NC_000913/1652089‑1652574                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                                                                                                                                                                                               <  1:261858/250‑1 (MQ=255)                ataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacagaa                                                                                                                                                                                                                               <  1:78453/251‑1 (MQ=255)                          tCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACTATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaaa                                                                                                                                                                                                                      <  1:190554/250‑1 (MQ=255)                                 ggTGAGCGGTGTGCACGTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAATACACAGGCGGGAACCAATCGGGTCATTCATCTTATGAAGCCATCAATCGACGCTCTCCGGATTCAGATGACATTAACGAGGCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAcc                                                                                                                                                                                                                <  2:288318/249‑1 (MQ=255)                                    gAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGt                                                                                                                                                                                                           <  2:355213/251‑1 (MQ=255)                                                cTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCCTTCATCTTATTAAGCCAGCAATCGACGCCCTCCGGAGTCAGATGACATTAAAGAGACTGAGCCAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAAGCACCTTTTTTTTT‑Caaacttc                                                                                                                                                                                              >  2:76105/1‑246 (MQ=255)                                                        ggaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTTAGCCAGCAATCGACGCTCCCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTTTTTTT‑CAACCTGAAGTGTCa                                                                                                                                                                                       >  1:246226/1‑250 (MQ=255)                                                           ggaCATAGATCTGAAAAATGGAACAATGATGCTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATAGGGTCATTTCTCTTATTAAAACCGCAATCGACGCTCTACGGAGTCAGATGACATTAACGAGGCTGAGCAAAGAGCATATCATTGATGTTCAACTCAGAGGGCATGGCCGGACAGAGAAAAAAAAATTCCCCTTTGTTTTTTCAAC‑TGAAGAGTCAGca                                                                                                                                                                                    >  1:26252/1‑249 (MQ=255)                                                           ggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCga                                                                                                                                                                                   <  1:120254/251‑1 (MQ=255)                                                              cATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGagag                                                                                                                                                                                <  1:309456/251‑1 (MQ=255)                                                                           aaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTAAACCAACAATCGACGCTCTCCGGGGGCAGATGGAATTAACGAGACTGAGCAAAGAGCCTATCATTGGTGTTTACCTCAGAGGGGATGGCAGAAAAGAGAAACAAAAATACACCTTTGTTTTT‑CCACCTTAAGTGAAAGCGGGAGTAAAAAAATATgg                                                                                                                                                                   >  2:250639/1‑251 (MQ=255)                                                                                  cAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAt                                                                                                                                                            <  1:263606/251‑1 (MQ=255)                                                                                  cAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCCTCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAAAGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCCGAGAGTATGGCAGAACAGAAAAACAAAAAATCACCTTTTGTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAt                                                                                                                                                            >  1:284508/1‑251 (MQ=255)                                                                                   aaTGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACAAATCGGGTCATTCATCTTATTAAGCCAGCAAACGACGCCCTCCGGAGTCAGATGACATTAAACAGACTGAGCAAAGAGCATATCTTTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCACCGAGAGAAAAAAAATATGGAGACCaat                                                                                                                                                           >  1:242805/1‑249 (MQ=255)                                                                                                        aGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCATCGAAAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAgg                                                                                                                                       <  2:49978/250‑1 (MQ=255)                                                                                                            aaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTTTTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGGATATCATTGGTGTTCACCTCAGAGAGt                                                                                                                                                                                                                                             >  1:354536/1‑145 (MQ=255)                                                                                                            aaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGTCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                                                                                                                                                                                             <  2:354536/145‑1 (MQ=255)                                                                                                                                           cGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTc                                                                                                   >  1:200046/1‑251 (MQ=255)                                                                                                                                                  cAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                                            <  2:200046/251‑1 (MQ=255)                                                                                                                                                          tcattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                                     >  2:218791/1‑250 (MQ=255)                                                                                                                                                                   ttattaAGCCAGCAATCGCCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCg                                                                           >  2:63550/1‑251 (MQ=255)                                                                                                                                                                                aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                                                              >  2:262169/1‑251 (MQ=255)                                                                                                                                                                                    gACGCTCTCCGGAGTCAGATGAAATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGATGACCATTTTACCGTTGACTCAATAAGGCAGATGTAGGACGCAGCGATCAAACGTGCCGGACTCCACCATCGAAAATCATATCAGTCGAGACATACTTATGCctg                                                          >  2:282169/1‑251 (MQ=255)                                                                                                                                                                                                          cATTAACGAGACTGAGCAAAGAGCATATCATTGATGCTCACCTCAGAGAGTATGGCAGAACAGAAAAACAACAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCt                                    >  1:127085/1‑251 (MQ=255)                                                                                                                                                                                                                        aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                      >  2:10654/1‑251 (MQ=255)                                                                                                                                                                                                                         gCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtata                     <  1:218791/251‑1 (MQ=255)                                                                                                                                                                                                                           aaaGAGAATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGc                   >  2:8001/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATg  <  2:193400/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     GGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTT‑CAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATG  >  NC_000913/1652089‑1652574

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