Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	210,799	T→C	9.1%	D5435D (GAT→GAC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	210,799	0	T	C	9.1%	120.5 / 2.3	43	D5435D (GAT→GAC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (24/15);  new base C (2/2);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGA  >  AE015451/210579‑211025                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   gTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGc                                                                                                                                                                                                                     <  1:337893/236‑1 (MQ=37)         cgCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCAt                                                                                                                                                                                                             <  1:345644/235‑1 (MQ=18)                                  tCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACAGGt                                                                                                                                                                                    >  1:102983/1‑236 (MQ=21)                                        tCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGc                                                                                                                                                                                                                     >  2:282958/1‑197 (MQ=33)                                                   cAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGGTTACAACGATGCTGGCCAGATc                                                                                                                                                                                                                        >  1:358915/1‑183 (MQ=14)                                                             tGGCCGAGAACGTCAAACCGACCTTCACTGGTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCAt                                                                                                                                                                                                             <  2:138106/184‑1 (MQ=18)                                                                                           ttCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGacac                                                                                                                            <  1:218995/235‑1 (MQ=38)                                                                                               cGTCATCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                                                                                 <  2:102983/162‑1 (MQ=18)                                                                                                   aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                                                                                 >  2:183263/1‑158 (MQ=33)                                                                                                       aGCCGCTGGCTCACGACCTGGTCGTGACCCGGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCaa                                                                                                                                                    <  2:237874/199‑1 (MQ=255)                                                                                                        gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGGAGACGCCACTGGt                                                                                                                                                                                    >  2:246670/1‑166 (MQ=21)                                                                                                         ccGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAgtggt                                                                                                             <  1:183263/236‑1 (MQ=38)                                                                                                          cGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAgtggt                                                                                                             <  1:405819/235‑1 (MQ=40)                                                                                                               gCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCaa                                                                                                                                                                                                          >  2:473281/1‑137 (MQ=12)                                                                                                                       ccTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcggc                                                                                                                                                                                                    >  2:524497/1‑135 (MQ=34)                                                                                                                             cGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACt                                                                                                                                                                                       >  2:620379/1‑142 (MQ=35)                                                                                                                                    cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCCTCgggt                                                                                                                                          >  1:236680/1‑178 (MQ=255)                                                                                                                                        gCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcggcg                                                                                                                                                                                                   >  2:631235/1‑119 (MQ=17)                                                                                                                                                    aGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACGTTCGAGAACCtcg                                                                                                                                                                 >  2:705089/1‑139 (MQ=2)                                                                                                                                                         aCCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCgaga                                                                                                                                                                       <  1:377113/130‑1 (MQ=38)                                                                                                                                                         aCCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCgaga                                                                                                                                                                       >  2:377112/1‑130 (MQ=38)                                                                                                                                                            aTCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTa                                                          >  1:733132/1‑236 (MQ=255)                                                                                                                                                              cAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGaa                                                                           >  1:278470/1‑217 (MQ=40)                                                                                                                                                                   ccGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAGCTGggcg                                                                                                                                                        >  1:456543/1‑135 (MQ=255)                                                                                                                                                                   ccGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATc                                                   >  1:670812/1‑236 (MQ=255)                                                                                                                                                                          gaCCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGGGCGACCTTCGAGAACCTGCAACTGGGCGGCa                                                                                                                                                     >  1:229405/1‑131 (MQ=255)                                                                                                                                                                                 gagcCGTCCACCCACGCCGCGCAGGGTGATGACGTGTACACCGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACGGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCACCGCTTCGGTTCAAGTGACCGACACCACCGACGAAgtggtg                                                                                                            <  2:137226/163‑1 (MQ=11)                                                                                                                                                                                                  tGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCga                                                                   >  2:690575/1‑189 (MQ=39)                                                                                                                                                                                                    cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACt                                                                                                >  2:593015/1‑158 (MQ=255)                                                                                                                                                                                                    cgcAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACt                                                                                                <  1:593016/158‑1 (MQ=255)                                                                                                                                                                                                    cgcAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCAGCGAGAt                                                                >  1:613565/1‑190 (MQ=39)                                                                                                                                                                                                    cgcAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGGCGGTCACCGAAGGCGGCGAGATCACCCACACCATCACGCTGACCAACCAAGACGGTCTGCCAATcaac                  >  1:125802/1‑236 (MQ=255)                                                                                                                                                                                                       tgggCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCa                                                                                                                                                                                                           >  2:524310/2‑48 (MQ=38)                                                                                                                                                                                                       aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCCACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAgtgggg                                                                                                            >  1:74014/1‑141 (MQ=255)                                                                                                                                                                                                       aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCaaa                                                                                                                                                   >  2:222007/1‑103 (MQ=255)                                                                                                                                                                                                       aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                             <  1:698951/190‑1 (MQ=255)                                                                                                                                                                                                       aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                             >  2:698950/1‑190 (MQ=255)                                                                                                                                                                                                          gCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCaacaac               <  1:103150/233‑1 (MQ=255)                                                                                                                                                                                                          gCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCa                                                  >  2:672547/1‑198 (MQ=40)                                                                                                                                                                                                               gACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGc       <  1:284937/236‑1 (MQ=255)                                                                                                                                                                                                               gACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGc       <  1:111503/236‑1 (MQ=40)                                                                                                                                                                                                                  gTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATTAACAACCACAGTGCAt     <  1:302322/235‑1 (MQ=37)                                                                                                                                                                                                                    ttACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGa  <  1:672548/236‑1 (MQ=255)                                                                                                                                                                                                                             tGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGAACTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGGCACCgcc                                                                           >  2:111503/1‑152 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   GTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGTGCATTGA  >  AE015451/210579‑211025

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