Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	205,789	C→T	9.3%	G3765G (GGC→GGT)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	205,789	0	C	T	9.3%	105.8 / 4.9	43	G3765G (GGC→GGT)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base T (2/2);  total (21/22)
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.43e-02

CTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCC  >  AE015451/205568‑206008                                                                                                                                                                                                                              |                                                                                                                                                                                                                            cTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGccactg                                                                                                                                                                                                                   <  2:479842/232‑6 (MQ=255)          ggcggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                                                                                                   <  2:147756/223‑1 (MQ=34)                     ttCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTACTGCCACCg                                                                                                                                                                                                                   <  2:171110/212‑1 (MQ=255)                     ttCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACg                                                                                                                                                                                           <  1:335294/236‑1 (MQ=35)                      tCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCa                                                                                                                                                                                                                            <  1:400605/202‑1 (MQ=34)                                    cGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGa                                                                                                                                                                              <  2:167065/234‑1 (MQ=37)                                                              aGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGt                                                                                                                                                   <  1:467855/235‑1 (MQ=38)                                                                  gACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGa                                                                                                                                                                                                      >  2:500188/1‑180 (MQ=35)                                                                        gACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCtgtcg                                                                                                                                                                      >  1:269337/1‑204 (MQ=11)                                                                           ggtcgcGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGAc                                                                                                                                      >  1:61568/6‑235 (MQ=25)                                                                            ccGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGt                                                                                                                                    <  1:411745/236‑1 (MQ=35)                                                                             cGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCg                                                                                                                                                                                                        <  2:470705/167‑1 (MQ=33)                                                                             cGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCg                                                                                                                                                                                                        >  1:470709/1‑167 (MQ=33)                                                                                         aaGGCGGCGAGATCACCTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTGACTGCCCCGGACAACGTCTACGTTGGCACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGATGTCGGCAAg                                                                                                                             >  2:304501/1‑230 (MQ=11)                                                                                                                aTCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg                                                                                                                                                                                                                   <  1:158844/121‑1 (MQ=32)                                                                                                                    cGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACACCgc                                                                                                                                                                                                >  2:75867/1‑135 (MQ=17)                                                                                                                    cGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTg                                                                                                                 <  1:185510/215‑1 (MQ=39)                                                                                                                    cGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTg                                                                                                                 >  2:185506/1‑215 (MQ=38)                                                                                                                     gCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCaa                                                                                                                    >  1:337253/1‑211 (MQ=37)                                                                                                                       tGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGt                                                                                                                                    <  2:418130/193‑1 (MQ=34)                                                                                                                            aacaaAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTCATCATATCGATCGCAACCGTTGGAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACCCCGGTCAGCa                                                                                        >  1:415652/1‑232 (MQ=16)                                                                                                                            aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatcga                                                                                                                                                             <  1:280345/163‑1 (MQ=21)                                                                                                                            aacaaAGACGGCCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatcga                                                                                                                                                             >  2:280341/1‑163 (MQ=21)                                                                                                                                      gTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGAcc                                                                                                                                                                            >  2:281410/1‑138 (MQ=32)                                                                                                                                          gCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTc                                                                                                                                                                                                                 >  2:415129/1‑97 (MQ=33)                                                                                                                                           ccGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAAcc                                                                                                                                                     >  2:504847/1‑156 (MQ=31)                                                                                                                                            cGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGCAACCGGGGGCACTGCTACCGTCACCGCCCCGCACAACg                                                                                                                                                                                                 >  2:18000/1‑111 (MQ=11)                                                                                                                                                              cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGt                                                                                                                                    <  2:347288/154‑1 (MQ=33)                                                                                                                                                              cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTc                                                                                                                          <  1:457858/164‑1 (MQ=35)                                                                                                                                                               aTTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCAc                                                                                       >  1:110396/1‑198 (MQ=34)                                                                                                                                                                    ccTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGc                                               <  2:478331/233‑1 (MQ=255)                                                                                                                                                                      ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGTACCGTGGGTACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGGTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGAc                                                                                                        >  1:41328/1‑174 (MQ=255)                                                                                                                                                                      ttCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCg                                              <  1:3348/232‑1 (MQ=255)                                                                                                                                                                                  gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAg                                                                                                     <  1:381553/165‑1 (MQ=35)                                                                                                                                                                                  gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGGCAc                              >  1:73404/1‑236 (MQ=255)                                                                                                                                                                                  gACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAg                                                                                                     >  2:381549/1‑165 (MQ=18)                                                                                                                                                                                         aaaCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGt                                                                 <  2:148468/194‑1 (MQ=255)                                                                                                                                                                                            ccGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGt                                                                 >  2:302922/1‑191 (MQ=255)                                                                                                                                                                                                  tcaCCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGAc                 >  2:491226/1‑233 (MQ=255)                                                                                                                                                                                                           cGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTAcc                                                                             <  2:418560/164‑1 (MQ=21)                                                                                                                                                                                                           cGGCCAACGGCACCGTGGGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTAcc                                                                             >  1:418564/1‑164 (MQ=21)                                                                                                                                                                                                               taatggCACCGTGGGCACTGCCACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCg                                                                                                                                         <  2:244860/96‑1 (MQ=12)                                                                                                                                                                                                               cAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGcc  <  1:302926/235‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                            CTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGACGTCGGCAAGTTCGAGCAACTGACGCTGGACAAGACGCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCC  >  AE015451/205568‑206008

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