Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	206,782	T→C	9.5%	A4096A (GCT→GCC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	206,782	0	T	C	9.5%	177.4 / 2.5	68	A4096A (GCT→GCC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (5/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACA  >  AE015451/206550‑207011                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      aGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTgg                                                                                                                                                                                                                                     <  1:1119908/235‑1 (MQ=33)         aCGCTGGACAAGACCCCGGTCAGCACTACCGTCGCCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAg                                                                                                                                                                                                                            <  1:115143/236‑1 (MQ=17)           gCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGCTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAgcg                                                                                                                                                                                                                          <  1:610062/236‑1 (MQ=18)                   aGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTg                                                                                                                                                                                                                                      >  2:651749/1‑216 (MQ=33)                       cccGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACAcc                                                                                                                                                                                                               <  1:881007/235‑1 (MQ=33)                             cAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCt                                                                                                                                                                                                         <  1:947903/235‑1 (MQ=33)                               gCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAACCCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTgc                                                                                                                                                                                                       <  1:876968/235‑1 (MQ=17)                                    aCCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCggtgaaa                                                                                                                                                                                                                                <  2:422162/205‑6 (MQ=17)                                        tCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTgc                                                                                                                                                                                                       <  2:493118/226‑1 (MQ=32)                                            cgacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCt                                                                                                                                                                                                         <  2:586374/220‑1 (MQ=33)                                             gacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGt                                                                                                                                                                                        >  1:93264/1‑236 (MQ=32)                                             gacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGt                                                                                                                                                                                        >  1:948963/1‑236 (MQ=33)                                             gacgaGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGAc                                                                                                                                                                                          >  2:915312/1‑234 (MQ=18)                                                gggCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCa                                                                                                                                                                                                                             <  2:1036664/194‑1 (MQ=18)                                                               ggCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCa                                                                                                                                                                      <  1:696792/236‑1 (MQ=32)                                                                 cAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCAGAGCAACAACGCCCAGGTCACCATCAAGGCTGGc                                                                                                                                                                                                                                    >  1:219447/1‑172 (MQ=18)                                                                      cGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCCCCACCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACGTGCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAgcg                                                                                                                                                                                                                          <  2:123349/177‑1 (MQ=9)                                                                          cggcAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGc                                                                                                                                                                            <  2:880954/219‑1 (MQ=31)                                                                           cgcagcaACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGAc                                                                                                                                                                                          <  1:595283/203‑1 (MQ=14)                                                                                caACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAgcgcg                                                                                                                                                                                                                        <  2:60907/169‑1 (MQ=31)                                                                                    gAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACt                                                                                                                                                  >  1:946893/1‑235 (MQ=31)                                                                                          gACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAggg                                                                                                                                                                                                 >  2:654085/1‑182 (MQ=32)                                                                                          gACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                           <  1:280842/236‑1 (MQ=31)                                                                                           aCCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGt                                                                                                                                           <  1:651750/235‑1 (MQ=31)                                                                                           aCCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTa                                                                                                                                                                                                                   <  1:404188/163‑1 (MQ=14)                                                                                           aCCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTTCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTa                                                                                                                                                                                                                   >  2:404187/1‑163 (MQ=9)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCCTCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAac                                                                                                                                                                                 >  1:53353/1‑179 (MQ=14)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGc                                                                                                                                                                                                                                    >  2:75406/1‑128 (MQ=32)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTgc                                                                                                                                                                                                       >  2:30715/1‑157 (MQ=31)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTgc                                                                                                                                                                                                       <  1:30715/157‑1 (MQ=31)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGAcc                                                                                                                        >  1:1067767/1‑236 (MQ=31)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGAc                                                                                                                         >  1:976345/1‑235 (MQ=31)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTCCGCAGGGTGATGACGTCTa                                                                                                                                                                                     >  1:1087713/1‑175 (MQ=14)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCa                                                                                                                                                                                                                             >  2:997129/1‑135 (MQ=12)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCa                                                                                                                                                                      >  1:538586/1‑190 (MQ=9)                                                                                                             tCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACAcc                                                                                                                                                                                                               >  2:729075/1‑149 (MQ=17)                                                                                                                  gCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGTTGCGACCTTCg                                                                                                                    >  1:492389/1‑235 (MQ=14)                                                                                                                        cAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTTTACAACGACGCTGGCCAGAt                                                                                                                                                                   >  1:966239/1‑182 (MQ=14)                                                                                                                             ttCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTg                                                                                                            <  2:673070/232‑1 (MQ=9)                                                                                                                             gtCGGGGGCCGAGAACGTCAACCCGACGTTCCCCGTACACGTCACCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTa                                                                                                                                                                                     <  2:180244/158‑1 (MQ=9)                                                                                                                                           aCGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCcgcg                                                                                          <  1:915313/236‑1 (MQ=35)                                                                                                                                             gTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGaa                                                                                        <  1:1128596/236‑1 (MQ=35)                                                                                                                                                 aaCCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAgg                                                                                      <  1:714864/234‑1 (MQ=37)                                                                                                                                                                             gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCg                                                        >  2:594420/1‑236 (MQ=37)                                                                                                                                                                                gCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCg                                                        >  2:1115057/1‑233 (MQ=37)                                                                                                                                                                                    gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCCACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGACCCTGGCGCTGGGCGGTGCCGCGCAGGTggg                                                                                  >  2:86189/1‑202 (MQ=255)                                                                                                                                                                                    gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCaa                                                 <  1:841042/236‑1 (MQ=35)                                                                                                                                                                                    gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCa                                                  >  1:1082907/1‑235 (MQ=35)                                                                                                                                                                                    gCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGAGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCaa                                                 <  1:898103/236‑1 (MQ=25)                                                                                                                                                                                       cACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCaa                                                 <  2:948962/233‑1 (MQ=35)                                                                                                                                                                                         cGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCt                                                                                                       <  1:33742/177‑1 (MQ=33)                                                                                                                                                                                         cGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGtt                                              <  2:976344/234‑1 (MQ=35)                                                                                                                                                                                                  ccagACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTgg                                                                                   <  1:313124/185‑1 (MQ=37)                                                                                                                                                                                                   gTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCCTCCAGAACCTGGAGCTGGGCGGTGCCGCGAAGGTGGATGTCAcc                                                                           >  2:494161/1‑195 (MQ=18)                                                                                                                                                                                                         cTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGATGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTg                                                                                    >  2:627242/1‑180 (MQ=39)                                                                                                                                                                                                                    cGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAg                                                           >  2:425770/1‑194 (MQ=255)                                                                                                                                                                                                                    cGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAg                                                           <  1:425771/194‑1 (MQ=255)                                                                                                                                                                                                                       ccAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTg                                                                                                      >  2:855718/1‑148 (MQ=34)                                                                                                                                                                                                                       ccAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc               >  1:883410/1‑235 (MQ=35)                                                                                                                                                                                                                            tCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATc         <  1:1115058/236‑1 (MQ=34)                                                                                                                                                                                                                              aCCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAAGTCGGCGGTGGACGCCACCGGTGCGACCTTCGAGAACCTGCCACTGGGCGGTGCcg                                                                                            >  1:549203/1‑151 (MQ=255)                                                                                                                                                                                                                              aCCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCaa       >  2:782575/1‑235 (MQ=21)                                                                                                                                                                                                                              aCCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc       >  1:83780/1‑236 (MQ=33)                                                                                                                                                                                                                               ccATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCAc       <  1:1002981/235‑1 (MQ=34)                                                                                                                                                                                                                                   cAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGaa                                                                                        >  2:34826/1‑150 (MQ=37)                                                                                                                                                                                                                                   cAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGa                      >  2:959938/1‑216 (MQ=34)                                                                                                                                                                                                                                   cAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTaca  >  1:495676/1‑236 (MQ=35)                                                                                                                                                                                                                                      ggCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTaca  <  2:1067766/233‑1 (MQ=37)                                                                                                                                                                                                                                         tGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGCCGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAggcggc               <  1:447065/217‑1 (MQ=21)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      AGCAACTGACGCTGGACAAGACCCCGGTCAGCACTACCGTCACCGACGAGCCTGGCACCCCAGGCAATCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTACACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTTTACAACGACGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTGGCGCTGGGCGGTGCCGCGAAGGTGGATGTCACCGACACCACCGACGAAGTCGTGGCCAAGTTGACCGCCACCCCGTCGGTCACCGAAGGCGGCGAGATCACCTACA  >  AE015451/206550‑207011

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