Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	588,045	A→G	72.4%	D711D (GAT→GAC)	nfrA ←	bacteriophage N4 receptor, outer membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	588,045	0	A	G	72.4%	45.2 / 17.4	29	D711D (GAT→GAC)	nfrA	bacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (2/6);  new base G (11/10);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACG  >  W3110S.gb/587977‑588104                                                                     |                                                            tGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                           <  1:134511/71‑1 (MQ=255) tGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                           <  1:910489/71‑1 (MQ=255) tGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                           <  1:1575184/71‑1 (MQ=255) tGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                           <  1:3507610/71‑1 (MQ=255) tGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATc                                                           <  1:415923/71‑1 (MQ=255)  gCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                           >  1:132150/1‑70 (MQ=255)  gCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                           >  1:3257/1‑70 (MQ=255)   cAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCaa                                                         >  1:2784644/1‑71 (MQ=255)     aaCGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTc                                                           <  1:67841/67‑1 (MQ=255)           ggAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATc                                                           <  1:832428/61‑1 (MQ=255)            gAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCatcaa                                                <  1:2042445/71‑1 (MQ=255)            gAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCatcaa                                                <  1:2452441/71‑1 (MQ=255)            gAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCatcaa                                                <  1:2229645/71‑1 (MQ=255)                     gCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGcc                                       >  1:2222283/1‑71 (MQ=255)                     gCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGcc                                       >  1:343975/1‑71 (MQ=255)                          gttgATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCg                                  >  1:3022877/1‑71 (MQ=255)                           ttgATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGt                                 <  1:2757004/71‑1 (MQ=255)                              aTTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGt                              <  1:3001551/71‑1 (MQ=255)                               ttttGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGt                              <  1:951633/70‑1 (MQ=255)                                                  gggTTATCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCa          >  1:2852568/1‑71 (MQ=255)                                                     ttATCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCg                     <  1:1991973/57‑1 (MQ=255)                                                       aTCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATc        >  1:666428/1‑68 (MQ=255)                                                       aTCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATc        >  1:1391116/1‑68 (MQ=255)                                                       aTCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATc        <  1:143224/68‑1 (MQ=255)                                                        tCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAg     <  1:134774/70‑1 (MQ=255)                                                        tCAGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGa    >  1:2643960/1‑71 (MQ=255)                                                          aGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACg  >  1:2167967/1‑71 (MQ=255)                                                          aGCGCCTGATTGTCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACg  >  1:1894422/1‑71 (MQ=255)                                                             gcCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCAt         >  1:1718004/1‑61 (MQ=255)                                                                     |                                                            TGCAAACGGCGGAAATTGAAGCGTTGTTGATTTTGTTCTGGGGTCAGTGGGGTTATCAGCGCCTGATTATCAATGTCATCAATCACCAGCCGGGCGTAGTGCTGCGTCGCAGGCATGTCATCCAGACG  >  W3110S.gb/587977‑588104

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