Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	2,550,125	G→A	43.7%	noncoding (1436/1522 nt)	PP_16SF →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	2,550,125	0	G	A	43.7%	48.9 / 85.9	64	noncoding (1436/1522 nt)	PP_16SF	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base A (16/12);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATC  >  AE015451/2550052‑2550358                                                                          |                                                                                                                                                                                                                                          ccGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                         >  1:655125/1‑236 (MQ=11)   gggCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAcc                                                                       >  1:687274/1‑236 (MQ=31)       cTTGTACACACCGCCCGTCACATCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                         <  2:1175080/230‑1 (MQ=11)       cTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGt                                                                   <  1:1222342/236‑1 (MQ=31)         tGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGt                                                                   <  2:776899/234‑1 (MQ=31)             cacacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGa                                                                         <  2:687271/224‑1 (MQ=11)             cacacCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAgg                                                             <  1:1078983/236‑1 (MQ=31)                  cGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTg                                                                    <  2:1195431/224‑1 (MQ=31)                   gCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                       <  1:126136/236‑1 (MQ=14)                     ccGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                       >  2:1227854/1‑234 (MQ=31)                     ccGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATGGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTg                                                         >  1:320716/1‑232 (MQ=14)                      cGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAAcgcgc                                                                         >  1:11083/1‑214 (MQ=11)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTtata                                                               >  2:186101/1‑224 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                       <  1:735852/232‑1 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                       >  2:735849/1‑232 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAg                                                      >  2:516980/1‑233 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAg                                                      >  2:993063/1‑233 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAg                                                      >  1:414598/1‑233 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCt                                                    >  1:337519/1‑235 (MQ=31)                       gTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTc                                                   >  1:875138/1‑236 (MQ=14)                        tCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGc                                                     >  2:908865/1‑233 (MQ=31)                        tCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCa                                                  >  1:243616/1‑236 (MQ=31)                             ccATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGt                                                                   <  1:256548/214‑1 (MQ=31)                               aTGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCt                                                    >  2:559735/1‑227 (MQ=14)                                 gggAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTag                                         >  1:1192917/1‑236 (MQ=31)                                 gggAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTc                                                           >  1:228067/1‑218 (MQ=14)                                  ggAGTGGGTTGCACCAGAAGTAGCTAGTCTATCCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                        <  1:130389/236‑1 (MQ=14)                                  ggAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTCGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                        <  1:769123/236‑1 (MQ=14)                                   gAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCCGTTGGTTagag                                       >  2:113237/1‑236 (MQ=14)                                   gAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTagag                                       >  1:387388/1‑236 (MQ=31)                                    aGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                        <  2:399735/234‑1 (MQ=31)                                      tGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                    >  1:452214/1‑236 (MQ=31)                                      tGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                    <  1:781583/236‑1 (MQ=31)                                         gTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGATCTCACACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCAcc                                 >  1:769289/1‑236 (MQ=14)                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCTACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAggac                                                           >  2:1179711/1‑206 (MQ=14)                                           tGCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTa                                                       >  2:1241586/1‑212 (MQ=14)                                            gCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCt                              <  1:870601/236‑1 (MQ=14)                                            gCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCt                              <  1:520460/236‑1 (MQ=14)                                            gCACCAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCt                              <  1:309881/236‑1 (MQ=14)                                              aCCAGGAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCt                              <  2:835044/234‑1 (MQ=14)                                               ccAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTaga                                        <  2:1114643/223‑1 (MQ=14)                                                cAGGAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCt                                                    <  1:657606/210‑1 (MQ=14)                                                cAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATa                          >  1:423455/1‑236 (MQ=31)                                                cAGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTAAGTTGGTTAGAGCGCACCCCTGATa                          >  1:1279285/1‑236 (MQ=14)                                                 aGAAGTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATaa                         >  1:336267/1‑236 (MQ=14)                                                     gTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGGTCAGTTGGTTag                                         >  2:1119053/1‑216 (MQ=14)                                                     gTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAggg                      >  1:904358/1‑235 (MQ=14)                                                     gTAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGt                     >  1:428965/1‑236 (MQ=14)                                                      tAGCTAGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGATg                    >  1:590335/1‑236 (MQ=14)                                                           aGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAggg                      >  2:1245575/1‑229 (MQ=31)                                                           aGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTg                                                         <  2:948311/194‑1 (MQ=17)                                                           aGTCTAACCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTg                                                         >  1:948314/1‑194 (MQ=14)                                                                aaCCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGAACCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAg          >  1:498112/1‑236 (MQ=14)                                                                aaCCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCa           >  1:882876/1‑235 (MQ=14)                                                                 aCCTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATAGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTg                                                                    <  2:1112746/177‑1 (MQ=14)                                                                  ccTTCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGTCGATGCTGCAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCa           <  2:655122/233‑1 (MQ=14)                                                                   gttCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGt                                                                   <  2:176399/175‑1 (MQ=14)                                                                   cTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGt                                                <  2:387385/195‑1 (MQ=32)                                                                    ttCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCa      >  1:911222/1‑236 (MQ=31)                                                                    ttCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGtt        >  2:872536/1‑234 (MQ=14)                                                                    ttCGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCa      <  1:1258151/236‑1 (MQ=14)                                                                    tccGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGACCCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTtata                                                               <  2:152005/177‑1 (MQ=14)                                                                     tCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgc                                      <  1:149184/203‑1 (MQ=31)                                                                      cGGGAGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTta                                                                 <  1:421051/175‑1 (MQ=14)                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAAcgcg                                                                          >  2:1166228/1‑163 (MQ=11)                                                                         ggggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATc  >  1:871220/1‑235 (MQ=31)                                                                          gggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAgcgc                                    >  2:1123399/1‑200 (MQ=31)                                                                          gggACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATc  >  1:700295/1‑234 (MQ=31)                                                                          |                                                                                                                                                                                                                                          CCGGGCCTTGTACACACCGCCCGTCACACCATGGGAGTGGGTTGCACCAGAAGTAGCTAGTCTAACCTTCGGGGGGACGGTTACCACGGTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTAAAAGACGATGCTGTAACGCGACCCTGTTATAGGTCTGTAGCTCAGTTGGTTAGAGCGCACCCCTGATAAGGGTGAGGTCGGCAGTTCAAATC  >  AE015451/2550052‑2550358

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