Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	205,801	C→G	10.2%	V3769V (GTC→GTG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	205,801	0	C	G	10.2%	108.3 / 2.8	38	V3769V (GTC→GTG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (17/17);  new base G (2/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAA  >  AE015451/205588‑206019                                                                                                                                                                                                                      |                                                                                                                                                                                                                           ttCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACg                                                                                                                                                                                                      <  2:523214/236‑1 (MQ=34)                  acacCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGCACCAACGAcc                                                                                                                                                                                       >  2:397125/1‑234 (MQ=17)                                             tGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGACAACGGAACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatcgatc                                                                                                                                                                      <  1:263244/224‑1 (MQ=37)                                              gACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGcc                                                                                                                                                                                                                      >  2:330822/1‑175 (MQ=35)                                                          gTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGCCCTTCACGCTGAGCGACGGCAACACCGTCATCCCCGTGCCGGCCAACGGCGCCGTGGGCACTGCTACCGTCAATGCCCCGGACAACGTCt                                                                                                                                                                                                         <  2:382740/176‑1 (MQ=255)                                                                  ccGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGATGTCGGCAAGTTc                                                                                                                                     <  1:357766/236‑1 (MQ=37)                                                                               gaTCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGCTCAACAACCACAGTGCATTGACCTTCACGCTGAGCGCCGGCAAAACCGTCATCACCGTGCCGCCCAACGGCCCCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCgg                                                                                                                                                                                                  <  1:494966/162‑1 (MQ=255)                                                                                            aTCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatc                                                                                                                                                                          <  1:16834/173‑1 (MQ=9)                                                                                                cGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTAc                                                                                                                                                                                               <  2:658633/148‑1 (MQ=35)                                                                                                cGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTAc                                                                                                                                                                                               >  1:658634/1‑148 (MQ=35)                                                                                                cGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGAcccc                                                                                                                                                                                     >  2:146819/1‑157 (MQ=18)                                                                                                cGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGATGAAGGt                                                                                                                                                       >  2:159765/1‑188 (MQ=37)                                                                                                 gCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTTCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGcc                                                                                                                                                                                                                      <  1:340204/124‑1 (MQ=17)                                                                                                 ccTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGcc                                                                                                                                                                                                                      <  1:500890/123‑1 (MQ=18)                                                                                                      ccAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTg                                                                                                                      >  2:670091/1‑215 (MQ=38)                                                                                                        aacaaAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTc                                                                                               >  1:386456/1‑236 (MQ=39)                                                                                                                     tGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGG‑CAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGta                                                                                                                                                             >  2:480800/1‑159 (MQ=14)                                                                                                                                  ccACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCa                                                                                                      >  2:485786/1‑203 (MQ=35)                                                                                                                                    acaGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCACATCGATCGCAACCGTTGAAGGCGCGGACGTCGGCAAGTTCGAGCAACTGACGc                                                                                                                        >  1:152235/1‑183 (MQ=25)                                                                                                                                    acaGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAg                                                                                                                <  1:548583/191‑1 (MQ=35)                                                                                                                                            ttGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGa                                                              >  2:381077/1‑233 (MQ=255)                                                                                                                                                 cTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGt                                                                            <  1:291125/214‑1 (MQ=255)                                                                                                                                                  tccACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAgg                                                           <  1:229369/228‑1 (MQ=255)                                                                                                                                                     aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTg                                                                                                                            >  2:287856/1‑162 (MQ=38)                                                                                                                                                     aCGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTg                                                                                                                            <  1:287857/162‑1 (MQ=38)                                                                                                                                                       gCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGt                                                  >  2:455738/1‑234 (MQ=255)                                                                                                                                                              gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAg                                                                                                                                                         <  2:361547/124‑1 (MQ=32)                                                                                                                                                              gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACa                                                                                                                  >  2:430414/1‑163 (MQ=35)                                                                                                                                                              gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGAAGCTGGACAAGACGCCGGTCAGCActt                                                                                                >  1:27221/1‑179 (MQ=21)                                                                                                                                                                                       cGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGt                                                                                                                                                       >  1:196942/1‑101 (MQ=35)                                                                                                                                                                                       cGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGt                                                                                                                                                       <  2:196941/101‑1 (MQ=35)                                                                                                                                                                                                 cACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGACCGTCACCGACGAGCCAGGCACCCCGGGCa                                                                  >  2:538542/1‑176 (MQ=255)                                                                                                                                                                                                     gTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTc    <  2:667297/234‑1 (MQ=255)                                                                                                                                                                                                       gggCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCaa  <  2:386455/234‑1 (MQ=255)                                                                                                                                                                                                            cTGCAACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGt                                                                                                                                       <  2:515140/96‑1 (MQ=18)                                                                                                                                                                                                              gcaaCCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACg                                                                         <  2:269202/153‑1 (MQ=255)                                                                                                                                                                                                              gcaaCCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACg                                                                         >  1:269203/4‑156 (MQ=255)                                                                                                                                                                                                                  ccGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGACGGTCACCGACGAGCCAgg                                                                             >  1:150936/1‑148 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                           TTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAA  >  AE015451/205588‑206019

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: ‑