Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	781,790	C→T	50.0%	L273L (TTG→TTA)	hlyE ←	hemolysin E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	781,790	0	C	T	50.0%	‑3.4 / 72.8	80	L273L (TTG→TTA)	hlyE	hemolysin E
Reads supporting (aligned to +/- strand):  ref base C (24/16);  new base T (9/31);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 1.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTTTCTGATACTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGTTTTTAT  >  minE/781736‑781852                                                       |                                                               tCTTTTCTGATACTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTaa                                                  >  1:838197/1‑69 (MQ=255) tCTTTTCTGATACTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTa                                                   >  1:3383130/1‑68 (MQ=255) tCTTTTCTGATACTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTa                                                   >  1:1876518/1‑68 (MQ=255)            aCTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATcatcat                                          <  1:2476686/66‑1 (MQ=255)            aCTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATcatcat                                          <  1:1160714/66‑1 (MQ=255)            aCTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATcatcat                                          <  1:1813925/66‑1 (MQ=255)            aCTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATcatcat                                          <  1:556700/66‑1 (MQ=255)            aCTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATcatcat                                          <  1:3144326/66‑1 (MQ=255)            aCTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATcatcat                                          <  1:1616208/66‑1 (MQ=255)                        tgtTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATc                           >  1:3017716/1‑69 (MQ=255)                        tgtTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATc                           >  1:47478/1‑69 (MQ=255)                        tgtTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATc                           >  1:1824370/1‑69 (MQ=255)                        tgtTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATc                           >  1:793426/1‑69 (MQ=255)                        tgtTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATc                           >  1:299376/1‑69 (MQ=255)                        tgtTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATc                           >  1:1053354/1‑69 (MQ=255)                            aaTCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATcatcat                                          <  1:3237945/50‑1 (MQ=255)                               cATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAg                               <  1:855219/58‑1 (MQ=255)                                aTTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGgttgt                    >  1:2534923/1‑68 (MQ=255)                                aTTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGgttg                     >  1:3264992/1‑67 (MQ=255)                                   tttttGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCa                <  1:2817200/69‑1 (MQ=255)                                   ttttgGGCCGCTTC‑‑TTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:807714/69‑1 (MQ=255)                                    ttttGGCAGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCa                <  1:2986791/68‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2209746/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2803612/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:3110586/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2510380/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2514163/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2520538/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2585054/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:27821/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2158755/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:907919/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:903050/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:1478534/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:1674504/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:764921/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:661650/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:557489/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:192508/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:3370397/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:3365376/66‑1 (MQ=255)                                      tgggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:203615/66‑1 (MQ=255)                                       gggCCGCTTCTTTTAGTAAAGAAAGCCTTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2270223/66‑1 (MQ=255)                                       gggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:1337081/66‑1 (MQ=255)                                       gggCCGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2240281/66‑1 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:3158230/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:995164/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:1799705/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:2001681/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:1260589/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:1600014/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:2198830/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGaa                             >  1:1368550/1‑48 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGa                              >  1:1699422/1‑47 (MQ=255)                                           cGCTTCTTTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:1594557/63‑1 (MQ=255)                                              ttcttTTAGTAAAGAAAGCATTAAATCATCATAATCATCGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:1845476/69‑1 (MQ=255)                                              ttcttTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:2718069/69‑1 (MQ=255)                                              ttcttTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:43425/69‑1 (MQ=255)                                              ttcttTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:2593312/69‑1 (MQ=255)                                              ttcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:2375911/69‑1 (MQ=255)                                              ttcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:3434722/69‑1 (MQ=255)                                              ttcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:181343/69‑1 (MQ=255)                                              ttcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:571337/69‑1 (MQ=255)                                               tcttTTAGTAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGt              <  1:2006585/59‑1 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:1300265/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:918358/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    <  1:9858/69‑1 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:1151050/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:1901554/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:213490/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:2233434/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:306774/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGttttt    >  1:2607116/1‑69 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     >  1:1714877/1‑68 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     >  1:3353759/1‑68 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     >  1:2020506/1‑68 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     >  1:3145023/1‑68 (MQ=255)                                               tcttTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGtttt     >  1:1241504/1‑68 (MQ=255)                                                  tttAGTAAAGAAAGCATTAAATCATCATAATCAACGAAGAATCTGGTTGTTTCAGTTTCCGtttt     <  1:695084/65‑1 (MQ=255)                                                  tttAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGTTTTTAt  <  1:2705276/68‑1 (MQ=255)                                                       |                                                               TCTTTTCTGATACTCATTACAGGTGTTAATCATTTTTTTGGCCGCTTCTTTTAGCAAAGAAAGCATTAAATCATCATAATCAACGTAGAATCTGGTTGTTTCAGTTTCCGTTTTTAT  >  minE/781736‑781852

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