Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	197,872	A→G	20.6%	K1126K (AAA→AAG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	197,872	0	A	G	20.6%	70.3 / 8.0	33	K1126K (AAA→AAG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (20/6);  new base G (5/2);  total (25/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAA  >  AE015451/197656‑198101                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      cacaGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACAGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                                                                                    >  1:408478/1‑236 (MQ=255)                                      aCGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                                                                                        <  2:301052/211‑1 (MQ=255)                                                               aCTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTGGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAgc                                                                                                                                                                                                              >  1:173483/1‑179 (MQ=255)                                                                          accacTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCa                                                                                                                                                  >  1:489647/1‑229 (MQ=255)                                                                             acaaccgCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATc                                                                                                                                                            >  2:298228/7‑216 (MQ=255)                                                                                 ctGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTc                                                                                                                                    >  1:391242/1‑236 (MQ=255)                                                                                  tGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                                                                                    >  1:375379/1‑155 (MQ=255)                                                                                            aCCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCg                                                                                                                                                                              >  1:255129/1‑183 (MQ=255)                                                                                            aCCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGGAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                                                                                        >  2:129402/1‑157 (MQ=255)                                                                                              cGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCTGCCaa                                                                                                                                                                                                    >  1:267810/1‑159 (MQ=255)                                                                                               gCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGc                                                                                                                                                                     >  1:31945/1‑189 (MQ=255)                                                                                                      aaaGCGTGACCGAAGTCGGCGTGATCACCTACACCGTGACCCTGAGCATTCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTATAACCCAGGGCAGCGTCTATTTCCAGACCCCCGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCg                                                                                                                                                           <  2:267809/191‑1 (MQ=255)                                                                                                              aCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTacaa                                                                                                                                                                                        <  2:375378/155‑1 (MQ=35)                                                                                                                                          tGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCACTGACGTCTACAACAACGGCTCGACCGTCAGCGTCCCCATCGAGAGCGCCACGGGCGGCAACTTCGAACAt                                                                                                                              >  2:44266/1‑185 (MQ=255)                                                                                                                                            aCCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATAc                                                                         <  2:391241/236‑1 (MQ=255)                                                                                                                                               cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                                                                                       >  1:270682/1‑139 (MQ=255)                                                                                                                                               cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCa                                                                                                                                                                       <  2:270681/139‑1 (MQ=255)                                                                                                                                               cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCCGACTCCGAATCCGACCCCAGCCTCGACTg                                                                                               >  1:47356/1‑211 (MQ=255)                                                                                                                                               cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACTGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGAACATCTGACTCCGAACCCGACCCCAg                                                                                                        >  2:483864/1‑202 (MQ=11)                                                                                                                                                          tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGACCGCCACGGGCGGCaagt                                                                                                                                      >  2:268237/1‑159 (MQ=11)                                                                                                                                                                    ccGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCCGCGTCACCATCGAGAgcgc                                                                                                                                                    >  2:416577/1‑137 (MQ=38)                                                                                                                                                                    ccGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGtcatca                                                                                         >  1:222632/1‑196 (MQ=255)                                                                                                                                                                         gacagTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGa                                                                                     >  2:599031/5‑195 (MQ=255)                                                                                                                                                                         aaCCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGa                                                                                                                          >  1:220344/1‑158 (MQ=40)                                                                                                                                                                           cagTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAgc                                                   <  1:599032/225‑1 (MQ=255)                                                                                                                                                                                       ccAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGt                                                                      <  2:222632/196‑1 (MQ=255)                                                                                                                                                                                        cAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTg                                                        >  1:22051/1‑209 (MQ=255)                                                                                                                                                                                                   ccATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCCATGGCAATGTGACCGAAGACCAGCAGccttc                   >  1:152029/1‑235 (MQ=255)                                                                                                                                                                                                   ccATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGa                                                                                                  >  2:536284/1‑156 (MQ=255)                                                                                                                                                                                                               aaGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGt                                                                      >  2:38892/1‑172 (MQ=255)                                                                                                                                                                                                               aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGAc                                                                                                             <  2:643916/133‑1 (MQ=255)                                                                                                                                                                                                               aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACAGGCGGCAACTTCGAACATCTGACTCCGAATCCGAc                                                                                                             >  1:643917/1‑133 (MQ=255)                                                                                                                                                                                                                 gCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCCAGAGCGCCACGggc                                                                                                                                             >  1:88074/1‑99 (MQ=18)                                                                                                                                                                                                                     gcaaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTaa  <  1:609654/232‑1 (MQ=255)                                                                                                                                                                                                                       aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTCAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTAc             >  1:124832/1‑221 (MQ=255)                                                                                                                                                                                                                       aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGa                                                                                                  >  2:664713/1‑136 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      CACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAACTTCGAACATCTGACTCCGAATCCGACCCCAGCCTCGACTGTCATCAATGACAGCATCGATACCGTCACCGTAAGTATTGTCAGCAATGGCAATGTGACCGAAGACCAGCAGCCTTCCTTTACTGTGAAAGTAA  >  AE015451/197656‑198101

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