Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,763,869	(T)7→9	coding (1285/1617 nt)	abpA ←	CP4‑57 prophage; anti‑bacteriophage protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,763,862	1	.	T	89.5%	57.7 / 0.5	19	coding (1292/1617 nt)	abpA	CP4‑57 prophage; anti‑bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (6/11);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01
*	NC_000913	2,763,862	2	.	T	89.5%	57.6 / 0.5	19	coding (1292/1617 nt)	abpA	CP4‑57 prophage; anti‑bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (6/11);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCT  >  NC_000913/2763785‑2763932                                                                               ||                                                                       cATCGAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑ttt                                                                     >  1:497789/1‑81 (MQ=255)     gAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑ttttttt                                                                 >  1:497811/1‑81 (MQ=255)     gAAAGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑tttttt                                                                  <  1:497769/80‑1 (MQ=255)        aGGTTGGCTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑tttttttttg                                                              >  1:497804/1‑78 (MQ=255)               cTGCCATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTTtgcgcgcgg                                                       >  1:497807/1‑78 (MQ=255)                  ccATCAAGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTAt                                                      <  1:497770/81‑1 (MQ=255)                       aaGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTGCGCGCGGTATCGATGc                                                >  1:497805/1‑80 (MQ=255)                       aaGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGAt                                                  <  1:497772/80‑1 (MQ=255)                       aaGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGAt                                                  <  1:497773/80‑1 (MQ=255)                       aaGAATCTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGAt                                                  <  2:497771/80‑1 (MQ=255)                                   aaTGTCGTGGCGAAGCAGGTAGTTGTCGTCTTTAATGTCGAGAATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAAt                                      <  1:497774/80‑1 (MQ=255)                                         gTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCg                               >  1:497806/1‑81 (MQ=255)                                                 gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACt                        >  2:497816/1‑80 (MQ=255)                                                 gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTg                       <  1:497776/81‑1 (MQ=255)                                                 gCAGGTAGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTg                       <  2:497775/81‑1 (MQ=255)                                                       aGTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCg                 >  2:497808/1‑81 (MQ=255)                                                        gTTATCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCg                 <  1:497777/80‑1 (MQ=255)                                                            tCGTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTAtct            >  1:497825/1‑81 (MQ=255)                                                              gTCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTcc          <  2:497778/81‑1 (MQ=255)                                                               tCTTTAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCg         <  2:497779/81‑1 (MQ=255)                                                                   tAATGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGc       <  1:497780/79‑1 (MQ=255)                                                                     aTGTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATc   >  2:497814/1‑81 (MQ=255)                                                                       gTCGAGGATTTTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCt  >  2:497818/1‑80 (MQ=255)                                                                               ||                                                                       CATCGAAAGGTTGGCTGCCATCAAGAATTTCGTCAATGTCGTGGCGAAGCAGGTAGTTATCGTCTTTAATGTCGAGGA‑‑TTTTTTTGCGCGCGGTATCGATGCAATCCGAAATGTCCTCGTACAACTGATCCCGTATCTCCGGCAATCT  >  NC_000913/2763785‑2763932

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑