Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	210,799	T→C	12.3%	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	12.3%	104.1 / 2.2	40	D5435D (GAT→GAC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (21/13);  new base C (2/3);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAAC  >  AE015451/210572‑211012                                                                                                                                                                                                                                    |                                                                                                                                                                                                                      gaGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCa                                                                                                                                                                                                               <  1:4771/236‑1 (MQ=35)        gTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGc                                                                                                                                                                                                        <  1:702999/236‑1 (MQ=37)             ccAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAg                                                                                                                                                                                                              <  1:506482/225‑1 (MQ=21)                          gCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcg                                                                                                                                                                                         <  2:342427/233‑1 (MQ=34)                                          cACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGAt                                                                                                                                                                                                            >  2:782423/1‑198 (MQ=33)                                              ctCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAACCCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACt                                                                                                                                                                                           <  1:214011/210‑1 (MQ=21)                                                     ccgACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCgag                                                                                                                                                           >  1:634799/1‑236 (MQ=37)                                                                          agaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCaa                                                                                                                                       >  1:536036/1‑235 (MQ=39)                                                                                                 gTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGt                                                                                                                        >  1:255808/1‑227 (MQ=39)                                                                                                  ttCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGacac                                                                                                               <  1:742601/235‑1 (MQ=40)                                                                                                     aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGAcc                                                                                                                    >  2:456134/1‑227 (MQ=40)                                                                                                     aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCGTCGGTTCAAGTGACCGAcaccac                                                                                                            >  1:104209/1‑235 (MQ=38)                                                                                                      cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGc                                                                                                                                                                     <  2:79506/177‑1 (MQ=35)                                                                                                          aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAAc                                                                                                                                                        >  2:549008/1‑186 (MQ=38)                                                                                                          aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGAAGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGAAGCCACTGGTGCGACCTTCgagacc                                                                                                                                                        >  1:48541/1‑184 (MQ=16)                                                                                                             cAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGaa                                                                                                     >  2:417149/1‑234 (MQ=40)                                                                                                                       cTCACGACCTGGTCTTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACTCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACAACACCGACGAAATGGTGGCCa                                                                                           <  2:422424/234‑1 (MQ=255)                                                                                                                         cACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCtc                                                                                                                                                     >  1:295294/1‑173 (MQ=38)                                                                                                                         cACGACCAGGCCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGCCAGATCAGCCCGGGc                                                                                                                                                                                                  >  2:93462/1‑129 (MQ=11)                                                                                                                                   tCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGAcg                                                                             >  1:352549/1‑235 (MQ=40)                                                                                                                                   tCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAGGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcggcg                                                                                                                                                                                      >  2:140198/1‑131 (MQ=21)                                                                                                                                     gTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAgt                                                                                                   >  2:317464/1‑212 (MQ=255)                                                                                                                                                caacaaTGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCCggg                                                                                                                                                                                                   >  2:97488/1‑105 (MQ=17)                                                                                                                                                          ggggTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACcg                                                                                                          <  1:163819/182‑1 (MQ=39)                                                                                                                                                                aCCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCAc                                                <  1:554981/236‑1 (MQ=40)                                                                                                                                                                   aTCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                                                                    <  1:701081/101‑1 (MQ=38)                                                                                                                                                                   aTCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                                                                    >  2:701079/1‑101 (MQ=38)                                                                                                                                                                      gaggcaggcGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTg                                                                                  <  1:648772/195‑1 (MQ=255)                                                                                                                                                                          ctggCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACCGGGGCGAccct                                                                                                                                                               >  2:163817/3‑113 (MQ=11)                                                                                                                                                                          ccGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTcg                                                                                                                                                                                         >  1:434053/1‑89 (MQ=39)                                                                                                                                                                               gagaCCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGa                                                                                                                                                                                 >  2:262889/1‑92 (MQ=37)                                                                                                                                                                                agaCCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGCGCGACCCTCGAGAACCTGCa                                                                                                                                                   >  2:186537/1‑121 (MQ=255)                                                                                                                                                                                 gaCCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACGTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGCCACCGAAGGCGCCGAGATCACCTACACCAt                                       >  2:432684/1‑228 (MQ=255)                                                                                                                                                                                                         tGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGAcaccac                                                                                                            >  2:337814/1‑135 (MQ=40)                                                                                                                                                                                                             cAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCt               <  2:634797/224‑1 (MQ=40)                                                                                                                                                                                                              aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATca       >  2:296757/1‑231 (MQ=255)                                                                                                                                                                                                              aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCaacaac  >  1:341015/1‑236 (MQ=255)                                                                                                                                                                                                                         gtctACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATcaacaa   <  1:232401/221‑1 (MQ=255)                                                                                                                                                                                                                          tttACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACGGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGc             <  2:255806/213‑1 (MQ=39)                                                                                                                                                                                                                                gcGATGGTGGCCGGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGCCCtt                                                                                                                                                               <  2:257521/60‑1 (MQ=11)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                      GAGCCAGGTACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAAC  >  AE015451/210572‑211012

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