Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,110,903	(T)7→8	coding (285/786 nt)	ECB_01999 →	putative phage protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,110,896	1	.	T	97.1%	114.2 / ‑2.4	35	coding (278/786 nt)	ECB_01999	putative phage protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (18/16);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GAATCAATGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑TTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACAATAATTGTAAT  >  REL606/2110821‑2110965                                                                             |                                                                      gAATCAATGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTa                                                                        <  1:8232/76‑1 (MQ=255)  aaTCAATGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑t                                                                      >  1:322988/1‑76 (MQ=255)    tCAATGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑ttt                                                                    <  2:343637/76‑1 (MQ=255)      aaTGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑ttttt                                                                  <  1:1129999/76‑1 (MQ=255)       aTGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑tttttt                                                                 >  1:981275/1‑76 (MQ=255)       aTGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑tttttt                                                                 <  1:370227/76‑1 (MQ=255)       aTGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑tttttt                                                                 >  1:484812/1‑76 (MQ=255)       aTGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGACTGATTGATGGAAAACTAAAAACATCGTA‑tttttt                                                                 >  1:981237/1‑76 (MQ=255)               ggCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACt                                                          <  2:843865/76‑1 (MQ=255)                  aTCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGc                                                       <  2:191822/76‑1 (MQ=255)                    cGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCtt                                                     >  1:838608/1‑76 (MQ=255)                     gAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTc                                                    >  1:575574/1‑76 (MQ=255)                         aTGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTggg                                                >  1:73512/1‑76 (MQ=255)                               gaATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACa                                          <  2:1138060/76‑1 (MQ=255)                                 aTGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATa                                        >  1:255105/1‑76 (MQ=255)                                  tGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATaa                                       >  2:209323/1‑76 (MQ=255)                                   gATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAAc                                      >  2:706721/1‑76 (MQ=255)                                        agcagGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAac                                 <  1:493838/76‑1 (MQ=255)                                         gcagGGGCTGATTGATGGAAAACTAAATACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAacc                                <  2:89975/76‑1 (MQ=255)                                          cagGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAacca                               <  2:359948/76‑1 (MQ=255)                                          cagGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAacca                               >  2:197330/1‑76 (MQ=255)                                          cagGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAacca                               >  2:1139199/1‑76 (MQ=255)                                          cagGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAacca                               >  2:1139188/1‑76 (MQ=255)                                           agGGGCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAaccac                              >  2:1064819/1‑76 (MQ=255)                                            ggggCTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAaccacc                             <  1:154115/76‑1 (MQ=255)                                                cTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACcaccacgc                         >  2:709804/1‑74 (MQ=255)                                                cTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACcacac                         <  1:213756/76‑1 (MQ=255)                                                cTGATTGATGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACcacac                         <  2:1037175/76‑1 (MQ=255)                                                  gattgatGGAAAACTAAAAACATCGTA‑TTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATa                      >  2:1058226/1‑76 (MQ=255)                                                   attgatGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATa                      <  2:830830/76‑1 (MQ=255)                                                    ttgatGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATaa                     <  1:242094/76‑1 (MQ=255)                                                     tgatGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAAt                    >  1:468816/1‑76 (MQ=255)                                                      gatGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAAtt                   <  2:278930/76‑1 (MQ=255)                                                      gatGGAAAACTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAAtt                   <  2:278936/76‑1 (MQ=255)                                                             aaCTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACan            <  1:346808/76‑2 (MQ=255)                                                             aaCTAAAAACATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACaa            >  2:1039523/1‑76 (MQ=255)                                                                 aaaaaCATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACaataat        >  2:447525/1‑76 (MQ=255)                                                                   aaaCATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACAATAATTg      >  2:321767/1‑76 (MQ=255)                                                                     aCATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACAATAATTGTa    >  1:87214/1‑76 (MQ=255)                                                                     aCATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACAATAATTGTa    >  1:564088/1‑76 (MQ=255)                                                                     aCATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACAATAATTGTa    <  2:1094406/76‑1 (MQ=255)                                                                      cATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACAATAATTGTaa   <  1:873887/76‑1 (MQ=255)                                                                      cATCGTATTTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCA‑CACACATAATTAAAACAAATATTTTAAt  >  1:662260/1‑76 (MQ=255)                                                                             |                                                                      GAATCAATGTATCAGGCATCGAAGATGTGAGAATGATACAGCAGGGGCTGATTGATGGAAAACTAAAAACATCGTA‑TTTTTTTCTTACTCGCTTCTGGGGTGACATAACAAAACCACCACACATAATTAAAACAATAATTGTAAT  >  REL606/2110821‑2110965

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