Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,429,013	C→T	intergenic (+29/‑36)	lomR → / → stfR	Rac prophage; protein LomR_1/Rac prophage; putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,013	0	C	T	100.0%	68.8 / NA	22	intergenic (+29/‑36)	lomR/stfR	Rac prophage; protein LomR_1/Rac prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/11);  total (11/11)

ACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACC  >  NC_000913/1428884‑1429141                                                                                                                                  |                                                                                                                                 aCCCGACCCAATCCGCGGACATTGATATTGTTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTccg                                                                                                                               <  1:186973/133‑1 (MQ=255)       ccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGtt                                                                                                                         >  2:467567/1‑133 (MQ=255)       ccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGtt                                                                                                                         >  1:128746/1‑133 (MQ=255)          aaTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAg                                                                                                                      <  1:467567/133‑1 (MQ=255)                      tttatatTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACTTAGTGTTATGACTGTCCGCCGGTTCAGGCGGGCtttttt                                                                                                          <  2:81828/130‑1 (MQ=255)                           atTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                     <  2:302076/133‑1 (MQ=255)                               cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                     <  2:84496/129‑1 (MQ=255)                               cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                     <  2:476805/129‑1 (MQ=255)                               cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                     >  1:476805/1‑129 (MQ=255)                               cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                     >  1:84496/1‑129 (MQ=255)                                   tGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAatat                                                                                             >  1:466176/1‑133 (MQ=255)                                    ggaGGCCCCCTCAGTAGCGACTGTCGTACTGAAGGTTTCATGGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATg                                                                                            <  2:248854/131‑1 (MQ=255)                                    gAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATg                                                                                            >  1:22845/1‑133 (MQ=255)                                     aaGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATgg                                                                                           <  2:466176/133‑1 (MQ=255)                                         ccccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGGGAATATGGCAGt                                                                                       >  2:363998/1‑133 (MQ=255)                                          cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTa                                                                                      <  1:421882/133‑1 (MQ=255)                                          cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTa                                                                                      <  1:347977/133‑1 (MQ=255)                                                                          cATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTCCTGAAAGACGGCACAgg                                                      <  2:128746/133‑1 (MQ=255)                                                                                 ggTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCg                                               >  1:84203/1‑133 (MQ=255)                                                                                    gtCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCAAAGGAAAACCGGTa                                            <  2:404463/133‑1 (MQ=255)                                                                                                  cTGATTAGCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAt                              >  2:88169/1‑133 (MQ=255)                                                                                                         gCCAGGTAACACAGTGTTATGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTg                       <  2:476186/133‑1 (MQ=255)                                                                                                                            tGACAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGca    >  2:102324/1‑133 (MQ=255)                                                                                                                              aCAGTCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCTGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGcacc  >  2:27978/1‑133 (MQ=255)                                                                                                                              aCAGTCCGCCGGGTCAGGCGGGCTTTTTTGTGGGGTGGATATGGCAGTAAAGATGTCAGGTTAACTGAAACACGGCACAGGAAAACCGGTACAGAACTGCACAATTCAGCTGAAAGCAAAACGGAACAgaatc  >  2:373031/1‑129 (MQ=255)                                                                                                                                  |                                                                                                                                 ACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACC  >  NC_000913/1428884‑1429141

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