Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	3,544,639	A→C	7.8%	F266V (TTC→GTC)	PP_3130 ←	conserved exported protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	3,544,639	0	A	C	7.8%	213.8 / 9.8	64	F266V (TTC→GTC)	PP_3130	conserved exported protein of unknown function
Reads supporting (aligned to +/- strand):  ref base A (34/25);  new base C (3/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

GCTTGCGCAGCGATTCGGTGCTGAGCTGCTGGGTATCGCGCACCACGGTGCCGTCGCCAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCTTGAGCTGCTGGCGGCCGGCGCT  >  AE015451/3544411‑3544858                                                                                                                                                                                                                                     |                                                                                                                                                                                                                             gCTTGCGCAGCGATTCGGTGCTGAGCTGCTGGGTATCGCGCACCACGGTGCCGTCGCCAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAAc                                                                                                                                                                                                                             >  1:198505/1‑230 (MQ=255)                      tGAGCTGCTGGGTATCGCGCACCACGGTGCCGTCGCCAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGc                                                                                                                                                                                                               <  2:818655/222‑1 (MQ=255)                                ggTATCGCGCACCACGGTGCCGTCGCCAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACa                                                                                                                                                                                                                            >  1:919455/1‑200 (MQ=255)                                              cGGTGCCGTCGCCAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACaacaccccgg                                                                                                                                                                                                 <  1:1040770/212‑10 (MQ=255)                                                    cGTCGCCAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGtg                                                                                                                                                                           >  1:282280/1‑227 (MQ=255)                                                         ccAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCcggcgg                                                                                                                                                                              >  2:504408/1‑220 (MQ=255)                                                             aGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTg                                                                                                                                                             <  2:609690/233‑1 (MQ=255)                                                                     gggCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAg                                                                                                                                                  <  1:416109/236‑1 (MQ=255)                                                                            cccAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGGTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAAc                                                                                                                                                                                                                             >  2:51987/1‑155 (MQ=255)                                                                                   cGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAg                                                                                                                                     >  1:529382/1‑235 (MQ=255)                                                                                      gTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGttttt                                                                                                                                                                        <  1:1242462/197‑1 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGCAGCGTACCGGGGTGTTGTCGGCGaa                                                                                                                                                                                                                              >  1:28192/1‑142 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCggggt                                                                                                                                                                           >  1:316441/1‑190 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCt                                                                                                                                                                   <  2:516563/200‑1 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCt                                                                                                                                                                   >  1:516566/1‑200 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGtt                                                                                                                                                              >  1:599637/1‑205 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGtt                                                                                                                                                              <  2:599634/205‑1 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCTAGGTgcg                                                                                                                                >  1:423348/1‑235 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGGTTTTCCGCTggg                                                                                                                                                                >  1:326351/1‑201 (MQ=255)                                                                                        cGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGCTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGAGAGGGTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑Tcgt                                                                                                                                                                                                                       >  2:136051/1‑146 (MQ=255)                                                                                          tggtgATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGACCagttgat                                                                                                                                                                                                                      <  1:273536/148‑7 (MQ=255)                                                                                              gATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGGTTTTc                                                                                                                                                                       >  1:267417/1‑190 (MQ=255)                                                                                                tGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAg                                                                                                                       <  1:470543/236‑1 (MQ=255)                                                                                                 gCAGGCGCCCCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCcg                                                                                                                                                                                  <  1:221004/176‑1 (MQ=255)                                                                                                    ggCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTg                                                                                                                                  <  2:1173131/221‑1 (MQ=255)                                                                                                     gcgcgcCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCTAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGcc                                                                                                                                                                                                              <  1:695351/144‑1 (MQ=255)                                                                                                       gcACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGACCagttggg                                                                                                                                                                                                                      >  2:273535/1‑129 (MQ=255)                                                                                                          ccAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCAc                                                                                                                                            <  2:529379/205‑1 (MQ=255)                                                                                                           cAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCAc                                                                                                                                            <  1:312419/204‑1 (MQ=255)                                                                                                                      ccAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCg                                                                                                                                                                     >  2:1130454/1‑168 (MQ=255)                                                                                                                           gCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCATCTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGgttttg                                                                                                                                                                       >  1:298917/1‑160 (MQ=255)                                                                                                                           gCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCggg                                                                                               >  2:981086/1‑233 (MQ=255)                                                                                                                           gCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTGTTTCCGCTGCGTTGAGGAACGACAGGTcacccg                                                                                                                                          >  2:556242/1‑185 (MQ=255)                                                                                                                            ccGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGtttttt                                                                                                                                                                       >  1:232629/1‑159 (MQ=255)                                                                                                                                    cGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGt                                                                                             <  1:1007341/226‑1 (MQ=255)                                                                                                                                      tttGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCg                                                                                  <  1:51987/235‑1 (MQ=255)                                                                                                                                       ttGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGcc                                                                                                                                                                                                              >  2:728652/1‑110 (MQ=255)                                                                                                                                            ttgttgGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGGTTGTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGa                                                                                                                                      >  2:1275119/1‑177 (MQ=255)                                                                                                                                            gtgttgGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAg                                                                                                                       <  1:215281/191‑1 (MQ=255)                                                                                                                                                 ggCGGAGTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTg                                                                      <  1:145978/236‑1 (MQ=255)                                                                                                                                                          cgtcgtcgGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGcgctc                                                                          <  1:1217701/223‑1 (MQ=255)                                                                                                                                                            tcgtcgGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCg                                                                                                                                                                >  1:410946/1‑135 (MQ=255)                                                                                                                                                               tcgGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCCACTGGTCGCCGACAGCCTTGACGACGCCCCTCTTGCcggcgg                                                                                                                                                                              >  2:300045/1‑118 (MQ=255)                                                                                                                                                                  gggCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCcggcgg                                                                                                                                                                              >  2:51720/1‑115 (MQ=255)                                                                                                                                                                          gAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGc                                             <  1:407576/236‑1 (MQ=255)                                                                                                                                                                               gCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGtt                                                      <  2:765368/222‑1 (MQ=255)                                                                                                                                                                                     aGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACg                                     >  2:1080981/1‑233 (MQ=255)                                                                                                                                                                                        ccGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGa                                 <  1:504411/234‑1 (MQ=255)                                                                                                                                                                                               gTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGa                                                         >  2:1222392/1‑203 (MQ=255)                                                                                                                                                                                               gTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACgg                           >  2:504277/1‑233 (MQ=255)                                                                                                                                                                                                ttGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCtt                       <  1:856888/236‑1 (MQ=255)                                                                                                                                                                                                     gCCCCAGAAGCGTACCGGGGTGTTGTCGGCGACCAGTTGA‑TGTTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGttttt                                                                                                                                                                        <  2:625397/86‑1 (MQ=21)                                                                                                                                                                                                     gCCCCAGAAGCGTACCGGGGTGTTGTCGGCGACCA‑GTTGATGTTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGttttt                                                                                                                                                                        >  1:625400/1‑86 (MQ=25)                                                                                                                                                                                                         cAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCTTGAg                    >  2:772786/1‑230 (MQ=255)                                                                                                                                                                                                            aaGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTGGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTTCGCCAGTCCAGCCGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCg                                                          >  1:22227/1‑189 (MQ=255)                                                                                                                                                                                                             aGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCTTGa                     >  2:859409/1‑225 (MQ=255)                                                                                                                                                                                                                 tACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGcc                                            <  2:534397/198‑1 (MQ=255)                                                                                                                                                                                                                    cGGGGTGTTGTCGGCGACCA‑GTTGATGTTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGc                               >  2:699893/1‑208 (MQ=255)                                                                                                                                                                                                                    cGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCTTGAGCTGCTGGCGGCCGgcg    >  1:1047660/1‑235 (MQ=255)                                                                                                                                                                                                                     ggggTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCg                                           >  1:394056/1‑195 (MQ=255)                                                                                                                                                                                                                      gggTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCTTGAGCTGCTGGCGGCCGGCGCt  <  1:807170/235‑1 (MQ=255)                                                                                                                                                                                                                         tgttgtCGGCGAACA‑TCAACGGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAAc                                                                                                                                                       <  2:832555/83‑1 (MQ=255)                                                                                                                                                                                                                            tgtCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCGGCGt                                                                                                                                                               >  2:147490/1‑72 (MQ=255)                                                                                                                                                                                                                             gtCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGGCATCAGGCCAGGTGGCGGGGGCTGGCAGGCCGAAGCGCTCGCTGGgt                                                                   >  1:644665/1‑163 (MQ=255)                                                                                                                                                                                                                              tCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGtta                                                     >  2:1024309/1‑175 (MQ=255)                                                                                                                                                                                                                                    aaCA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCg                                                                                                                                       >  1:322179/1‑88 (MQ=255)                                                                                                                                                                                                                                    aaCA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCg                                                                                                                                       <  2:322178/88‑1 (MQ=255)                                                                                                                                                                                                                                    aaCA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTg                                                                                                                                  >  2:423318/1‑93 (MQ=255)                                                                                                                                                                                                                                    aaCA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTg                                                                                                                                  <  1:423321/93‑1 (MQ=255)                                                                                                                                                                                                                                    aaCA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGTGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCTTGAGCTGCTGGcggcc        >  1:659325/1‑215 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                             GCTTGCGCAGCGATTCGGTGCTGAGCTGCTGGGTATCGCGCACCACGGTGCCGTCGCCAAAGACGTTGGGGCTGACCCAGGGCGAGTCGTGGTGATGCAGGCGCACCAGGTTGAAACCCAGTGCCGAAAGGCGTTTGGCTTGTTGGCGGATTTCGTCGTCGGGGCTTTTGAACAGCGAGTAGGCCGACAGGTTGGTGCCCCAGAAGCGTACCGGGGTGTTGTCGGCGAACA‑TCAACTGGTCGCCGACAGCCTTGACGAAGCCACGCTTGCCGGCGGGTTTTTCCGCTGCGTTGAGGAACGACAGGTCCACCGGCGAGGTGCGCCAGTCCAGCTGGTCATCAGGCCAGGTGGCGGGGTCTGCCAGGCCGAAGCGCTCGCTGGGTGTCGGGCCGAGTTCGATGTCGCCGGTCACGTTGAGCACGGCCTTGAGCTGCTGGCGGCCGGCGCT  >  AE015451/3544411‑3544858

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