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

ACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGT  >  AE015451/210591‑211018                                                                                                                                                                                                                 |                                                                                                                                                                                                                            aCGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGc                                                                                                                                                                                                        >  1:279243/1‑230 (MQ=1)                                                        gaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAAc                                                                                                                                            <  1:235861/235‑1 (MQ=39)                                                                                            gCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCATGGGCACCAACTCGGCGGTAGACGCCACTg                                                                                                                                                                               >  1:19718/1‑164 (MQ=18)                                                                                                         gACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCGGCGCAGGGCGATGACGTTTACAACGACCCTGGCCAGATCAGCCTGGGCATCAACTcgg                                                                                                                                                                                              >  2:437185/1‑136 (MQ=255)                                                                                                                  gTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGCTGCGACCTCCGAGAACGTGCAACTGggcgg                                                                                                                                                >  1:419976/1‑173 (MQ=255)                                                                                                                        cTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGt                                                                           >  1:277661/1‑236 (MQ=255)                                                                                                                             caacaaTGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTc                                                                                                                                       >  2:404012/1‑171 (MQ=255)                                                                                                                              aacaaTGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGtg                                                                                                                                                                                         >  2:368380/1‑119 (MQ=21)                                                                                                                                      ccAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAggcggc                                                              <  1:404016/235‑1 (MQ=38)                                                                                                                                        aGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGCTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCAcc                                                                       >  1:454987/1‑224 (MQ=255)                                                                                                                                             aCCAGCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCGGCGCAGGGCGATGACGTTGACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACg                                                                                                                                                                                     <  2:65972/109‑1 (MQ=255)                                                                                                                                                aTCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGACCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCt                                                    >  1:292440/1‑235 (MQ=255)                                                                                                                                                                  aGCGCGCCGTACACGCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCt                                                                                                                                                                      <  2:75512/103‑1 (MQ=25)                                                                                                                                                                  aGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGGGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGGGACCGACACCACCGACGAAGGGGGGGCCAAGCTGACt                                                                                         >  1:99381/1‑180 (MQ=255)                                                                                                                                                                                ccACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCACCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGacac                                                                                                                     <  1:263075/138‑1 (MQ=255)                                                                                                                                                                                    gCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACcaaca                               >  1:145774/1‑220 (MQ=39)                                                                                                                                                                                      tGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGtg                                                                                                                                                                                         <  2:191837/64‑2 (MQ=18)                                                                                                                                                                                      tGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGac                                                                                                                       >  2:366234/1‑130 (MQ=40)                                                                                                                                                                                           aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCACCCcgacgg                                                                                                                                                                                           >  1:303646/1‑57 (MQ=11)                                                                                                                                                                                           aGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCaacaac        <  1:403636/236‑1 (MQ=255)                                                                                                                                                                                             ggCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAAcc       <  1:498225/235‑1 (MQ=40)                                                                                                                                                                                                gATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTg                                                                                                                                                                               <  1:340487/64‑1 (MQ=37)                                                                                                                                                                                                    aCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCAGTGGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATc              <  1:482643/221‑1 (MQ=255)                                                                                                                                                                                                        ttACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGt  <  1:162366/229‑1 (MQ=255)                                                                                                                                                                                                                 tGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGGTCAAGTGACCg                                                                                                                         >  2:98671/1‑101 (MQ=255)                                                                                                                                                                                                                 |                                                                                                                                                                                                                            ACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGATGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTAGACGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACTGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCAATCAACAACCACAGT  >  AE015451/210591‑211018

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