Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,414,019	C→T	G259R (GGG→AGG)	recT ←	Rac prophage, recombination and repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,414,019	0	C	T	100.0%	60.4 / NA	19	G259R (GGG→AGG)	recT	Rac prophage, recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/10);  total (9/10)

ACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATATTTGAACAGGCGACGAATAGCCGTTTTCTTTGCCATTTCT  >  NC_000913/1413886‑1414152                                                                                                                                      |                                                                                                                                      aCAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGt                                                                                                                                    <  2:446267/137‑1 (MQ=255) aCAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGt                                                                                                                                    >  1:446267/1‑137 (MQ=255)                           cATTTTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCCTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATGACACTGTATTCCCTGTTTAATACAGAGGAATCTGCAgaaccgctt                                                                                                        <  1:239698/139‑9 (MQ=255)                                tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAATACAGAGGa                                                                                                                        <  1:351179/118‑1 (MQ=255)                                tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAATACAGAGGa                                                                                                                        >  2:351179/1‑118 (MQ=255)                                tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAg                                                                                                   <  2:158386/139‑1 (MQ=255)                                             tttttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAACGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCa                                                                                      <  1:100336/139‑1 (MQ=255)                                                 tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGt                                                                                                                                    <  1:173366/89‑1 (MQ=255)                                                 tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGt                                                                                                                                    >  2:173366/1‑89 (MQ=255)                                                                                         ccAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAata                                          >  1:147368/1‑139 (MQ=255)                                                                                           aGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTg                                                                  >  1:318838/1‑113 (MQ=255)                                                                                           aGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTg                                                                  <  2:318838/113‑1 (MQ=255)                                                                                                              aaTTATCGATTACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCa                                                                       >  1:300933/1‑89 (MQ=255)                                                                                                              aaTTATCGATTACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCa                                                                       <  2:300933/89‑1 (MQ=255)                                                                                                                        tACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGc                                                                    <  1:120580/82‑1 (MQ=255)                                                                                                                        tACACTGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGc                                                                    >  2:120580/1‑82 (MQ=255)                                                                                                                             tGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCaa                                                          >  2:244647/1‑87 (MQ=255)                                                                                                                             tGTATTCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCaa                                                          <  1:244647/87‑1 (MQ=255)                                                                                                                                 ttCCCTGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATATTTGAACAGGCGACGAATAGCCGTTTTCTTTGCCATttct  >  2:550708/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                      ACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATATTTGAACAGGCGACGAATAGCCGTTTTCTTTGCCATTTCT  >  NC_000913/1413886‑1414152

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