Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	197,834	T→C	10.8%	L1114L (TTG→CTG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	197,834	0	T	C	10.8%	132.7 / 6.1	46	L1114L (TTG→CTG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base C (3/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCAC  >  AE015451/197610‑198036                                                                                                                                                                                                                                 |                                                                                                                                                                                                           aCGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGGGGCGTGATCACCTACGCCGTGACCATGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACgg                                                                                                                                                                                                  <  1:848407/235‑1 (MQ=255)                     aCTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTc                                                                                                                                                                                                        <  2:588226/209‑1 (MQ=255)                                                                                  cGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                           <  2:366500/177‑1 (MQ=255)                                                                                  cGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGc                                                                                                                >  1:327953/1‑236 (MQ=255)                                                                                                                                tGCGCCCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                           <  2:142367/131‑1 (MQ=255)                                                                                                                                  cGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAAc                                                                      >  1:257057/1‑230 (MQ=255)                                                                                                                                            cGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTCACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTCCCAGACCCCGGCCAATGACg                                                                                                                              >  2:312921/1‑164 (MQ=255)                                                                                                                                             gCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTCACCCTGAGCAATCCTGCCCACACGCCGGTAACCGTGCCCTTCTCCACCGACCAAAACATTACCGTTGAAGCGGGTGACACCCAGAGCAg                                                                                                                                                           >  2:33951/1‑134 (MQ=17)                                                                                                                                                   ccgaGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAAc                                                                                                                                                                                                    <  2:444962/84‑1 (MQ=255)                                                                                                                                                   ccAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGt                                                                                                                                                                                                         <  1:779160/82‑1 (MQ=255)                                                                                                                                                       ccGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                                                                                                <  2:84025/102‑1 (MQ=18)                                                                                                                                                       ccGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                           <  1:753617/107‑1 (MQ=255)                                                                                                                                                         gtcaCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACgg                                                                                                                                                                                                  <  2:314880/80‑1 (MQ=18)                                                                                                                                                                   gcggcgTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGACCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGAc                                            >  2:866206/1‑223 (MQ=255)                                                                                                                                                                     ggcgTGATCACCTACACCGTGACCCTGAGCAATCCGGCCCAGACGCCGGTAACCGTTACCTTGTCCAAAGACCAAAc                                                                                                                                                                                            >  2:297099/1‑77 (MQ=14)                                                                                                                                                                         taaTCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGtt                                                                                                                                                                                  <  2:4065/81‑1 (MQ=38)                                                                                                                                                                                   cacCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTc                                   >  2:666235/1‑216 (MQ=255)                                                                                                                                                                                   cacCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTc                                   <  1:666238/216‑1 (MQ=255)                                                                                                                                                                                       gTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAAc                                                                                                                                                                                                    <  2:600942/51‑1 (MQ=35)                                                                                                                                                                                         gACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATa         <  1:475975/236‑1 (MQ=255)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACgg                                                                                                                                                                                                  >  1:860563/1‑49 (MQ=38)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACgg                                                                                                                                                                                                  <  2:860560/49‑1 (MQ=38)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGc                                                                                                                                                                                                 >  1:852269/1‑50 (MQ=40)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGc                                                                                                                                                                                                 <  2:852266/50‑1 (MQ=38)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                  <  1:404318/113‑1 (MQ=40)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                  >  2:404315/1‑113 (MQ=40)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCAc                                                                                >  2:293519/1‑163 (MQ=38)                                                                                                                                                                                           cccTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGCCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACt                               >  2:528412/1‑212 (MQ=11)                                                                                                                                                                                             cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCgg                                                                                                                                                                           <  1:411344/70‑1 (MQ=40)                                                                                                                                                                                             cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTa                                                                                                                          >  1:171203/1‑119 (MQ=39)                                                                                                                                                                                             cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGggcg                                                                           >  2:193246/1‑166 (MQ=40)                                                                                                                                                                                                 gCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGt     >  1:275935/1‑232 (MQ=255)                                                                                                                                                                                                 cagaTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCAc  <  1:891021/232‑1 (MQ=255)                                                                                                                                                                                                  cAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGCCGATTTCCAGGCCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGggcgg                                                                          <  2:225462/162‑1 (MQ=255)                                                                                                                                                                                                   aaTCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCac                                                                                            >  2:94908/1‑142 (MQ=39)                                                                                                                                                                                                      ccTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACt                                                       <  2:327950/177‑1 (MQ=255)                                                                                                                                                                                                        tGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTACAACCCAGAGCAGCGTCGATTTCCAGACCCCGGCCAAt                                                                                                                                  >  2:702902/1‑100 (MQ=255)                                                                                                                                                                                                             aGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGAc                                            <  2:718070/181‑1 (MQ=255)                                                                                                                                                                                                             aGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGAc                                            >  1:718073/1‑181 (MQ=255)                                                                                                                                                                                                               ccGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGa                                 <  2:612880/189‑1 (MQ=255)                                                                                                                                                                                                                  ccGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGGGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGACCATCTGACTACGAATc                                                >  2:487038/1‑172 (MQ=255)                                                                                                                                                                                                                  ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCt                                                           >  2:891018/1‑161 (MQ=37)                                                                                                                                                                                                                   cGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCt                                    <  2:257054/183‑1 (MQ=255)                                                                                                                                                                                                                       aaCCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAg                                                                                                     >  2:188274/1‑114 (MQ=39)                                                                                                                                                                                                                       aaCCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGCCGGCAACTTCGAACAACTGACt                                                       >  1:103483/1‑160 (MQ=11)                                                                                                                                                                                                                                cTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGa                                                                 >  2:154799/1‑141 (MQ=35)                                                                                                                                                                                                                                 |                                                                                                                                                                                                           ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCAC  >  AE015451/197610‑198036

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