BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	196,234	T→C	12.1%	T580T (ACT→ACC)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	196,234	0	T	C	12.1%	313.9 / 21.0	126	T580T (ACT→ACC)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand): ref base T (49/60); new base C (8/8); total (57/69)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

TGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCA  >  AE015451/196009‑196457

tGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGAcc                                                                                                                                                                                                                        >  1:992084/1‑235 (MQ=255)

          ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGTAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAg                                                                                                                                                                                                             >  1:905379/1‑236 (MQ=255)

          ccGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGCCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAg                                                                                                                                                                                                             >  1:952808/1‑236 (MQ=255)

               aaCTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGAcc                                                                                                                                                                                                                        >  1:27333/1‑220 (MQ=255)

                      aCCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGCCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGt                                                                                                                                                                                                    <  2:952806/233‑1 (MQ=255)

                        ccTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGAc                                                                                                                                                                                                                         <  1:1215095/211‑1 (MQ=255)

                        ccTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGAc                                                                                                                                                                                                                         >  2:1215092/1‑211 (MQ=255)

                             ccAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGAcc                                                                                                                                                                                                                        >  2:1085231/1‑206 (MQ=255)

                               aaCGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAg                                                                                                                                                                                            >  1:617201/1‑232 (MQ=255)

                               aaCGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGc                                                                                                                                                                                                                >  1:149999/1‑212 (MQ=255)

                                    ccAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAGACCACCACcg                                                                                                                                                                                                                             >  2:510417/1‑194 (MQ=255)

                                        aCGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                                                                                                                                                               <  1:465548/236‑1 (MQ=255)

                                         cGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGGTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACt                                                                                                                                                                              >  1:644817/1‑236 (MQ=255)

                                             tACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTaca                                                                                                                                                                          <  1:750864/236‑1 (MQ=9)

                                                       gCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGACGGCGGCGTGATCACCTACACCGTGACCCt                                                                                                                                                                >  1:371456/1‑236 (MQ=255)

                                                         gggTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCt                                                                                                                                                                <  1:1017526/234‑1 (MQ=40)

                                                            caaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGcc                                                                                                                                                                                                           <  2:1268458/187‑1 (MQ=39)

                                                             aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                                                                                                               >  2:1473254/1‑215 (MQ=21)

                                                              aaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCaa                                                                                                                                                          <  1:1006513/235‑1 (MQ=37)

                                                                        gaGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCTGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTCAGCAATCCt                                                                                                                                                      <  2:783329/228‑1 (MQ=255)

                                                                        cAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGa                                                                                                                                               >  1:783331/1‑236 (MQ=40)

                                                                           cGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTaca                                                                                                                                                                          <  2:688935/206‑1 (MQ=39)

                                                                              cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTg                                                                                                                                                                                                                     >  1:150373/1‑160 (MQ=40)

                                                                              cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAg                                                                                                                                                             >  1:1137335/1‑216 (MQ=38)

                                                                              cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACAcctg                                                                                                                                                                      >  2:965835/1‑205 (MQ=255)

                                                                              cGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCAGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           >  2:1113906/1‑234 (MQ=255)

                                                                               ggtttCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTg                                                                                                                                                                                                                     <  2:13629/157‑1 (MQ=16)

                                                                                 tttCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGa                                                                                                                                      <  1:433817/236‑1 (MQ=39)

                                                                                  gtCCAGACCCCGCCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGt                                                                                                                                                                                                        <  2:472313/168‑1 (MQ=255)

                                                                                   tCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCTGCGTGATCACTTACACCGTGAccc                                                                                                                                                                 <  1:1382047/207‑1 (MQ=255)

                                                                                      aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTg                                                                                                                                 <  1:1126864/236‑1 (MQ=255)

                                                                                             ggCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGt                                                                                                                          <  1:1481787/236‑1 (MQ=255)

                                                                                             ccagaaTGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                                                                                                                                                     <  1:620189/173‑1 (MQ=9)

                                                                                               ccAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                                                                                                                                                                                               >  1:386141/1‑165 (MQ=40)

                                                                                                     aCGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa                                                                                                                                                              <  1:504799/192‑1 (MQ=38)

                                                                                                     aCGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa                                                                                                                                                              >  2:504797/1‑192 (MQ=40)

                                                                                                      cGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                                                                                                       <  1:638839/166‑1 (MQ=40)

                                                                                                          tACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACGCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                                                                                                       <  2:394283/162‑1 (MQ=16)

                                                                                                               caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAA‑TTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTa                                                                                                        >  1:1168383/1‑236 (MQ=16)

                                                                                                                aacGGTTCGACGGTCAGCGTCACCAACGAGAACGCCACAGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                                                                                                                                                         <  1:453292/154‑1 (MQ=255)

                                                                                                                 acGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACTGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                                                                                                                                                     <  2:1139369/157‑1 (MQ=38)

                                                                                                                 acGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGAccc                                                                                                                             >  1:492632/1‑213 (MQ=255)

                                                                                                                     ttCGACGGTCAGCGTCACCATCGAGAACGCCACTGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGa                                                                                                                                      <  2:980318/200‑1 (MQ=255)

                                                                                                                     ttCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGAccc                                                                                                                                                                 >  1:1450119/1‑173 (MQ=255)

                                                                                                                     ttCGACCGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGACCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                                                                                                                                                                                               <  2:329885/143‑1 (MQ=255)

                                                                                                                       cGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGaa                                                                                                >  1:1342164/1‑236 (MQ=255)

                                                                                                                          ccgTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACGTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGc                                                                                                                  <  2:535643/213‑1 (MQ=255)

                                                                                                                           ggTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGGGCAATCCTGCCCAGACGCCGGTAACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  1:520088/236‑1 (MQ=255)

                                                                                                                                gCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca                                                                                                                                                                          >  2:1474665/1‑153 (MQ=255)

                                                                                                                                  gTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAc                                                                                     >  1:1295320/1‑236 (MQ=37)

                                                                                                                                  gTCACCATCGAGAACGCCACGGGCGGTAA‑CTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGt                                                                                                                          <  2:836647/199‑1 (MQ=255)

                                                                                                                                  atCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTaca                                                                                                                                                                          <  2:566739/150‑1 (MQ=38)

                                                                                                                                   tCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  2:1342161/228‑1 (MQ=17)

                                                                                                                                   tCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGCCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                                                                                                                                                     <  2:61483/139‑1 (MQ=21)

                                                                                                                                     aCCATCGAGAACGCCACAGGCGGTAA‑CTTCGAGCAACTGACCCCGACCCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCt                                                                                                                                                                <  2:287504/158‑1 (MQ=9)

                                                                                                                                           gagaACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGa                                                                                                                                      <  2:417984/178‑1 (MQ=40)

                                                                                                                                                 gCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGt                                                                                                                                                                                                        <  2:652000/106‑1 (MQ=255)

                                                                                                                                                 gCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTacac                                                                                                                                                                         >  1:1240064/1‑137 (MQ=40)

                                                                                                                                                   cACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGc                                                                                                                                            <  2:1137332/164‑1 (MQ=40)

                                                                                                                                                    aCGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAGCGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAAc                                                                                     >  2:694988/1‑218 (MQ=25)

                                                                                                                                                    aCGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTg                                                                                                                                                                     >  2:314549/1‑138 (MQ=21)

                                                                                                                                                    aCGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTg                                                                                                                                                                     <  1:314551/138‑1 (MQ=21)

                                                                                                                                                      gggCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTaaa                                                                                       >  1:900483/1‑214 (MQ=255)

                                                                                                                                                        gCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAAc                                                                                                             <  2:617199/190‑1 (MQ=40)

                                                                                                                                                        gCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGtt                                                                                                   <  2:1320720/200‑1 (MQ=255)

                                                                                                                                                                tgtCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGGTCAACGACTCGGTCGACACCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                                                                                                                                                                                         <  1:1005044/105‑1 (MQ=255)

                                                                                                                                                                                 cGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACGCCACTACCGCGACCCTGAccg                                                                                                                                                                                                                 >  2:675315/1‑64 (MQ=255)

                                                                                                                                                                                    aCCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTg                                                                                                                                 <  1:382360/143‑1 (MQ=39)

                                                                                                                                                                                     cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTc                                                                                                                         <  2:366022/150‑1 (MQ=39)

                                                                                                                                                                                       cGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACAGTGAcc                                                                                                                              >  1:1193418/1‑143 (MQ=21)

                                                                                                                                                                                          cGCCGGCTCAGACCACGATCCACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:288460/127‑1 (MQ=25)

                                                                                                                                                                                           gCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGCCaacaacga                                     >  1:299497/1‑226 (MQ=21)

                                                                                                                                                                                            ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                                                                                                 >  1:359376/1‑87 (MQ=33)

                                                                                                                                                                                            ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                                                                                                 <  2:359374/87‑1 (MQ=33)

                                                                                                                                                                                            ccGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           <  2:42474/125‑1 (MQ=35)

                                                                                                                                                                                                   aGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                     <  2:583555/140‑1 (MQ=21)

                                                                                                                                                                                                   aGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                                     >  1:583557/1‑140 (MQ=21)

                                                                                                                                                                                                     aCCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgtg                    <  1:694990/235‑1 (MQ=37)

                                                                                                                                                                                                            tCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                                                                                                 <  1:219719/71‑1 (MQ=34)

                                                                                                                                                                                                            tCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                                                    <  2:817170/132‑1 (MQ=37)

                                                                                                                                                                                                            tCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                                                    >  1:817172/1‑132 (MQ=12)

                                                                                                                                                                                                                cGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAttt                                                                   >  1:1367620/1‑177 (MQ=25)

                                                                                                                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                                                                                                       >  1:980930/1‑59 (MQ=35)

                                                                                                                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                                                                                                       >  2:532577/1‑59 (MQ=38)

                                                                                                                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcggc                                                                                                                                                                                       <  1:532579/59‑1 (MQ=38)

                                                                                                                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                                           >  2:1203666/1‑103 (MQ=38)

                                                                                                                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTg                                                                                                  <  1:189529/144‑1 (MQ=39)

                                                                                                                                                                                                                  aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAgcgg                   >  1:224351/1‑220 (MQ=14)

                                                                                                                                                                                                                  aCTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTg                                                                                                                           <  2:362626/119‑1 (MQ=25)

                                                                                                                                                                                                                         cGACACCTCTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTg                                                                                                                                 <  2:817064/106‑1 (MQ=16)

                                                                                                                                                                                                                         cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt                                                                                                                                                                                     >  2:1456635/1‑54 (MQ=38)

                                                                                                                                                                                                                         cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                                                                                                               >  1:1140577/1‑60 (MQ=38)

                                                                                                                                                                                                                         cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 >  1:782855/1‑138 (MQ=40)

                                                                                                                                                                                                                         cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 <  2:1186501/138‑1 (MQ=21)

                                                                                                                                                                                                                         cGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                 <  2:1001871/138‑1 (MQ=18)

                                                                                                                                                                                                                           acacCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            >  1:1328679/1‑141 (MQ=38)

                                                                                                                                                                                                                           acacCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                                            <  2:1328676/141‑1 (MQ=17)

                                                                                                                                                                                                                              ccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                                                                                                 <  2:164668/53‑1 (MQ=35)

                                                                                                                                                                                                                              ccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATc                                                                                                                                                                                 >  1:164669/1‑53 (MQ=35)

                                                                                                                                                                                                                              ccacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCt                                                                                                                                                                >  1:1023766/1‑70 (MQ=35)

TGCCCAGACGCCGGTAACTGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAA‑TTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCA  >  AE015451/196009‑196457

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