Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	727,616	A→T	14.3%	N860Y (AAC→TAC)	rhsC →	Rhs family putative polymorphic toxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	727,616	0	A	T	14.3%	124.7 / 14.9	49	N860Y (AAC→TAC)	rhsC	Rhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base A (26/16);  new base T (5/2);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACC  >  CP009273/727594‑727763                       |                                                                                                                                                    gTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGtata                       <  1:695613/149‑1 (MQ=11) gTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGtata                       >  2:1859805/1‑149 (MQ=11)  tCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATAc                      <  1:553894/149‑1 (MQ=31)  tCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATAc                      <  2:571626/149‑1 (MQ=31)     gACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcac                   >  1:753260/1‑149 (MQ=32)      aCCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                  >  1:1564407/1‑149 (MQ=32)      aCCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                  <  2:4779296/149‑1 (MQ=32)      aCCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                  >  1:2202489/1‑149 (MQ=32)      aCCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                  >  1:2203171/1‑149 (MQ=32)        cgcgATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGa                >  1:3907553/1‑149 (MQ=32)        cgcgATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGa                >  1:1517648/1‑149 (MQ=32)         gcgATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAc               >  2:3078484/1‑149 (MQ=32)         gcgATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAc               >  2:3078332/1‑149 (MQ=32)         gcgATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAc               >  2:14620/1‑149 (MQ=32)           gATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAccc             <  2:2714317/149‑1 (MQ=32)           gATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAccc             <  2:2714288/149‑1 (MQ=32)           gATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAccc             <  1:1248901/149‑1 (MQ=32)            aTTACACCTGGTACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            <  1:4247967/149‑1 (MQ=17)            aTTACACCTGGTACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  2:4662303/1‑149 (MQ=17)            aTTACACCTGGTACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            <  1:3280829/149‑1 (MQ=17)            aTTACACCTGGTACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  2:4662165/1‑149 (MQ=17)            aTTACACCTGGTACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  2:1879061/1‑149 (MQ=17)            aTTACACCTGGTACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  2:4662702/1‑149 (MQ=17)            aTTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  1:517391/1‑149 (MQ=32)            aTTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  1:1637648/1‑149 (MQ=32)            aTTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  2:1766444/1‑149 (MQ=32)            aTTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  2:4155876/1‑149 (MQ=32)            aTTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  2:4155977/1‑149 (MQ=32)            aTTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg            >  1:3607437/1‑149 (MQ=32)             ttACACCTGGTACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCgg           >  2:2059312/1‑149 (MQ=17)             ttACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCgg           >  1:2627560/1‑149 (MQ=32)             ttACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCgg           >  1:234758/1‑149 (MQ=32)             ttACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCgg           <  1:2153392/149‑1 (MQ=32)              tACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGc          <  2:944250/149‑1 (MQ=32)               acacCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTCTACCACAGACCCGGCa         >  2:1028519/1‑149 (MQ=17)               acacCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCa         >  2:2117149/1‑149 (MQ=32)               acacCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCa         <  2:2803183/149‑1 (MQ=32)               acacCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCa         >  2:2115543/1‑149 (MQ=32)                  ccTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGt      <  2:488090/149‑1 (MQ=32)                    tGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTaa    >  2:4659623/1‑149 (MQ=32)                    tGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTaa    >  2:4659835/1‑149 (MQ=32)                     ggAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAc   <  1:3803263/149‑1 (MQ=32)                     ggAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAc   <  1:2910412/149‑1 (MQ=32)                     ggAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAc   <  1:1186506/149‑1 (MQ=32)                      gAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCCTCCCGTATACCACAGACCCGGCAGGTCAcc  >  1:3251508/1‑149 (MQ=17)                      gAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAcc  >  2:1688445/1‑149 (MQ=32)                      gAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAcc  <  2:1260633/149‑1 (MQ=32)                      gAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAcc  <  2:1260627/149‑1 (MQ=32)                      gAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAcc  >  1:3250735/1‑149 (MQ=32)                       |                                                                                                                                                    GTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACC  >  CP009273/727594‑727763

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