Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	589,934	0	T	G	44.1%	37.8 / 19.0	34	T341P (ACC→CCC)	nfrA	bacteriophage N4 receptor, outer membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (4/15);  new base G (15/0);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 2.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  NC_000913/589786‑590072                                                                                                                                                     |                                                                                                                                             aTCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCCAGCATTTCATTGGCGgggggggggg                                                                                                                                            >  1:340198/1‑151 (MQ=255)    aGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGa                                                                                                                                         <  1:554164/151‑1 (MQ=255)     gTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGTGGCGAg                                                                                                                                        <  2:393747/151‑1 (MQ=255)          cAGGTTAGTTGATCCAGGCGGGGAAGAATTGCCGGGTCTTTCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTTTTTCGGGTCGCCCCGCTGGCAGCATAACCCTCCTTAAGCATTTTATTGGCGgggggggggggggggg                                                                                                                                      >  1:531566/1‑142 (MQ=255)                  ttGATCCAGGCGGGTAAGATTTGCCGGGTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAAAGCTCCTCAAGCCTTTTATTGGCGggggggggggggggggggtttt                                                                                                                                >  1:819771/1‑134 (MQ=255)                      tCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATTGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGGCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGgggggggggggggggggggtttt                                                                                                                               >  1:296406/1‑130 (MQ=255)                            cGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGgtttt                                                                                                                                >  1:484014/1‑136 (MQ=255)                              ggTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGa                                                                                                               <  2:788630/151‑1 (MQ=255)                                 aaGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGAGCAGGTTTTTCGCTGCGTCGTACTGAcc                                                                                                             >  2:300212/1‑150 (MQ=255)                                 aaGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc                                                                                                             <  2:121110/150‑1 (MQ=255)                                 aaGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGGGCAGGTTTTTTGCTggcgcgtc                                                                                                                   >  1:785661/1‑137 (MQ=255)                                     tttGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGTTTTTTGGCTGCGGCGGACTGACcttttt                                                                                                        >  2:875763/1‑151 (MQ=255)                                     tttGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGGGCAGTGTTTTCGCTGCGTCGTACTGACcctctt                                                                                                        >  1:42718/1‑151 (MQ=255)                                        gCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTa                                                                                                     <  2:872368/151‑1 (MQ=255)                                         ccGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGCGTGTTTTTGCTGCGTCCGACTGAACCTCttttt                                                                                                     >  1:670788/1‑149 (MQ=255)                                         ccGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTTGCTGCGTCGTTCTTACCCTCtttt                                                                                                      >  1:567463/1‑149 (MQ=255)                                               cTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCGCCGCGTCGGACTGACCTtttttttaaaaca                                                                                               >  2:62708/1‑143 (MQ=255)                                                 tGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGGTTTTTTGCTGCGTGCTACTTACCTTCTTTTTAAAACAcc                                                                                             >  1:761657/1‑150 (MQ=255)                                                  gCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGGGGAGGTTTTTCGCTGCGGTGTACTGACCCTCtttttataaca                                                                                               >  1:510407/1‑140 (MQ=255)                                                     gATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAg                                                                                        <  2:609980/151‑1 (MQ=255)                                                      aTACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGc                                                                                       <  1:577603/151‑1 (MQ=255)                                                       tACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAACACCGGGAGCg                                                                                      >  2:820929/1‑151 (MQ=255)                                                                cGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTGTGCTGCGGCGTACTGACCTTCTTTTTAAAGCACCCGTAGGGGCGCGcc                                                                               >  1:276031/1‑149 (MQ=255)                                                                cGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCcaa                                                                             <  2:646192/151‑1 (MQ=255)                                                                       gACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATa                                                                      <  1:500755/151‑1 (MQ=255)                                                                          gCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTg                                                                   <  1:698414/151‑1 (MQ=255)                                                                           cAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTgg                                                                  <  2:484024/151‑1 (MQ=255)                                                                            agagCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGGTTTTTTGCTGCGTCGTACTTACCCTCTTTTTACAGCACCGGGAGCGTCGCGCCCACAACAAACTgg                                                                  >  1:506248/1‑150 (MQ=255)                                                                                    aGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGGAGCGTCGCGCCAACAACATAATTGC‑gggt                                                            >  2:548578/1‑145 (MQ=255)                                                                                      ccTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGGTTTTTTGCGGGGTCGTACTGACCTTCTTTTTACAGCACCGGGAGGGGGGCGCCCCCAAAATACAGgg‑gggggt                                                          >  1:483253/1‑143 (MQ=255)                                                                                                            gACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGCGGGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGTCGCGCCCACAACAAACTGGCGGG‑TTGTTGGCAAAACTTACCCTATAAtttcc                                >  2:22843/1‑148 (MQ=255)                                                                                                                gCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCAAACTGTACCGTATAATTCGCCAAct                           <  1:364114/151‑2 (MQ=255)                                                                                                                 cATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCAAACTGTACCGTATAATTCGCCAACGc                          <  1:506990/151‑1 (MQ=255)                                                                                                                                        ttGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GGTTTTTGGGCAAAC‑GTACCGCATAAATTGCCCACGCCCGAACGGGGGTatggccggatt   >  2:787902/1‑141 (MQ=255)                                                                                                                                          ggCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTAttt  <  1:884174/150‑1 (MQ=255)                                                                                                                                          ggCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTAttt  <  1:62707/150‑1 (MQ=255)                                                                                                                                           gCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCCAACTTGACCCTATAATACGCCCACGCCTGCACCGGGGGTGGGCTGTAtt   >  1:65251/1‑148 (MQ=255)                                                                                                                                           gCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCCAACCGTACCGGATAATTTGCCAACCCCCGCACGGGGGTAGCGCTGTAttt  >  1:914388/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                             ATCAGTTGCCAGGTTAGTTGATCCAGGCGGGTAAGATTTGCCGGTTCTTGCTGATACAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGC‑GG‑TTGTCGGCAAACTGTACCGTATAATTCGCCAACGCCTGAACGGGGTTAGCGCTGTATTT  >  NC_000913/589786‑590072

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