Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	5,681,378	A→C	5.1%	intergenic (‑96/+108)	PP_4986 ← / ← PP_4987	Channel protein, hemolysin III family/putative Chemotaxis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	5,681,378	0	A	C	5.1%	359.7 / 5.8	98	intergenic (‑96/+108)	PP_4986/PP_4987	Channel protein, hemolysin III family/putative Chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base A (53/40);  new base C (2/3);  total (55/43)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CCGCCGTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGC  >  AE015451/5681143‑5681603                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  ccgccgTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGa                                                                                                                                                                                                                                   >  1:961781/1‑236 (MQ=255)               gaAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGc                                                                                                                                                                                                                      <  1:1300697/235‑1 (MQ=255)                 aGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCt                                                                                                                                                                                                                     >  1:964616/1‑234 (MQ=255)                        gATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCAtt                                                                                                                                                                                                               <  1:111377/233‑1 (MQ=255)                                       gCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCc                                                                                                                                                                                             >  1:250624‑M1/1‑235 (MQ=255)                                                      ggCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGGTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCa                                                                                                                                                                                                      <  2:46329/212‑1 (MQ=255)                                                        ccACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATTCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGTGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCc                                                                                                                                                                                             >  2:771305‑M1/1‑218 (MQ=255)                                                           cGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagcccta                                                                                                                                                                         <  1:938628‑M1/236‑22 (MQ=255)                                                            gATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagcccta                                                                                                                                                                         <  1:1102174‑M1/235‑22 (MQ=255)                                                             aTCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctg                                                                                                                                                                                   <  1:446409‑M1/224‑12 (MQ=255)                                                                      aTGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGGGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCggcggttcatgcagtggcaatgaaagcctgggg                                                                                                                                                                                         <  2:1031461‑M1/209‑34 (MQ=255)                                                                         gCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGACGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacacggc                                                                                                                                                           >  1:106731‑M1/1‑201 (MQ=255)                                                                         gCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTg                                                                                                                                                                                                        >  2:1210356/1‑191 (MQ=255)                                                                         gCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCa                                                                                                                                                                                                      >  1:477946/1‑193 (MQ=255)                                                                          cagcaatAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGCCCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgac                                                                                                                                                                  <  1:589621‑M1/225‑29 (MQ=255)                                                                             cAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGc                                                                                                                                                                                                       <  1:149772/188‑1 (MQ=255)                                                                                tAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAAc                                                                                                                                                                                                 <  2:1139633/191‑1 (MQ=255)                                                                                   cagcGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagcccta                                                                                                                                                                         <  2:713568‑M1/212‑22 (MQ=255)                                                                                          gACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagacc                                                                                                                                          <  1:909650‑M1/236‑53 (MQ=255)                                                                                           aCGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGACGCGGCCCCAGGCTGTCATTGCCACTGCATGAAc                                                                                                                                                                                                 >  2:410725/1‑180 (MQ=255)                                                                                              gCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACCCAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactcc                                                                                                                                               <  2:760737‑M1/227‑48 (MQ=255)                                                                                                     ccAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCGGTTGTGATGCTGATGGCGTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                                                                                              <  2:49269/173‑1 (MQ=255)                                                                                                     ccAGGTGGGTCCAGGCGTGGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgct                                                                                                                                                                                    <  2:400509‑M1/183‑11 (MQ=255)                                                                                                      cAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaa                                                                                                                                                              <  1:161021‑M1/204‑33 (MQ=255)                                                                                                         gTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgc                                                                                                                            <  1:808237‑M1/235‑67 (MQ=255)                                                                                                            ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggcc                                                                                                                                  >  2:957892‑M1/1‑166 (MQ=255)                                                                                                            ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggcc                                                                                                                                  <  1:957894‑M1/226‑61 (MQ=255)                                                                                                            ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactc                                                                                                                                                >  2:687389‑M1/1‑166 (MQ=255)                                                                                                              tCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcaca                                                                                                                      <  1:388425‑M1/236‑73 (MQ=255)                                                                                                                 aGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCAtt                                                                                                                                                                                                               <  2:1017329/144‑1 (MQ=255)                                                                                                                 aGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCAtt                                                                                                                                                                                                               >  1:1017332/1‑144 (MQ=255)                                                                                                                 aGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggct                                                                                                                                                                                       >  2:1511406‑M1/1‑161 (MQ=255)                                                                                                                          aGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaa                                                                                                                                                              >  2:286642‑M1/1‑152 (MQ=255)                                                                                                                                cACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccacgccagaccaccaggccctgcgcccca                                                                                                                        >  2:1335481‑M1/1‑146 (MQ=255)                                                                                                                                   cGTAATACATGCACACACAGACCTCCTGAGGTGGGCTTGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccg                                                                                                                                                         <  1:441266‑M1/180‑38 (MQ=255)                                                                                                                                      aaTACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcca                                                                                                                                                                   >  2:141235‑M1/1‑140 (MQ=255)                                                                                                                                       aTACATGCAAACACAGCCCTCCTGAGGTGGGCTGGGTTCCTGTGGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCggcggttcatgcagtggcaatgaaagcctgggg                                                                                                                                                                                         <  2:620166‑M1/144‑34 (MQ=255)                                                                                                                                        tACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTTTGTCGCGAAAGGGGCGCGAAGCGGCCCCAgg                                                                                                                                                                                                                       >  2:15252/1‑113 (MQ=255)                                                                                                                                            tGCAAACACAGCCCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagacc                                                                                                                                          <  1:598175‑M1/186‑53 (MQ=255)                                                                                                                                              cAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACcgc                                                                                                                                                                                              >  1:710370/1‑132 (MQ=255)                                                                                                                                                aaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtccggttgcagatccgcggca                                                                                     >  1:173513‑M1/1‑130 (MQ=255)                                                                                                                                                aaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtca                                                                                     >  1:173115‑M1/1‑130 (MQ=255)                                                                                                                                                aaCACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCCTTCATTGCCACTGCATGaa                                                                                                                                                                                                  >  2:424951/1‑126 (MQ=255)                                                                                                                                                aaCACAGACCTCCTGAGGAGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatc                                                                                                                                                                       >  1:167576‑M1/1‑130 (MQ=255)                                                                                                                                                 acacAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcac                                                                                    <  1:1210359‑M1/235‑107 (MQ=255)                                                                                                                                                  cacaGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcccctcacat                                                                                                                     >  2:505066‑M1/1‑128 (MQ=255)                                                                                                                                                         ctcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccg                                                                                                                                                         >  2:598173‑M1/1‑121 (MQ=255)                                                                                                                                                          tcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACt                                                                                                                                                                                                         >  2:874727/1‑109 (MQ=255)                                                                                                                                                          tcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagccc                                                                          >  1:1486592‑M1/1‑120 (MQ=255)                                                                                                                                                          tcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgc                                                                                  <  1:1152654‑M1/228‑109 (MQ=255)                                                                                                                                                             tGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcacc                                                                               <  2:964614‑M1/228‑112 (MQ=255)                                                                                                                                                              gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaag                                                                                                           >  2:663994‑M1/1‑116 (MQ=255)                                                                                                                                                              gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccaga                                                                                                                                            >  2:463526‑M1/1‑116 (MQ=255)                                                                                                                                                                ggTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttcca                                                                     >  1:441078‑M1/1‑114 (MQ=255)                                                                                                                                                                    ggCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacaccaggccgctcccaccc                                                                                                                                                >  2:452003‑M1/1‑110 (MQ=255)                                                                                                                                                                      cTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGCACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccaga                                                                                                                                            >  1:76978‑M1/1‑108 (MQ=255)                                                                                                                                                                      cTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttcc                                                                      >  2:1046213‑M1/1‑108 (MQ=255)                                                                                                                                                                       tGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcc                                                                                                                                                     >  1:828207‑M1/1‑107 (MQ=255)                                                                                                                                                                         ggTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCCTGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaa                                                                                                            >  2:346384‑M1/1‑105 (MQ=255)                                                                                                                                                                          acagctgttGTGATGCTGATGGCTTCGTTCTTTGTGGGGGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAGCCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcga                                                               <  1:596002‑M1/227‑128 (MQ=255)                                                                                                                                                                            tCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgct                                                                                                                                                                                    <  2:167575‑M1/112‑11 (MQ=255)                                                                                                                                                                              cTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCggcggttcatgcagtggcaatgaaagcctggggccgcttcgcgcccctttcgcgacacaa                                                                                                                                                              >  1:431250‑M1/1‑72 (MQ=255)                                                                                                                                                                              cTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCggcggttcatgcagtggcaatgaaagcctggggccgcttcgcgcccctttcgcgacacaa                                                                                                                                                              <  2:431248‑M1/132‑61 (MQ=255)                                                                                                                                                                               tgttgtGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcac                                                                                    <  2:1422530‑M1/205‑107 (MQ=255)                                                                                                                                                                                gttgtGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccg                                                      <  2:902457‑M1/234‑137 (MQ=255)                                                                                                                                                                                  tgtgATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcc                                                                                                                                                     >  2:689000‑M1/1‑96 (MQ=255)                                                                                                                                                                                    tgATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgc                                                                                         <  2:1486589‑M1/195‑102 (MQ=255)                                                                                                                                                                                      aTGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGATAGGGCTGCGCCGCAGCCCCggcggttcatgcagtggcaatgaaagcctggggccgcttcgcgcccctttcgcgacacaaggccg                                                                                                                                                         >  2:1511057‑M1/1‑64 (MQ=255)                                                                                                                                                                                            aTGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAACGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGCACCGCcggggctgctgcgcagccctatcgcgacacaa                                                                                                                                                              >  2:14305‑M1/1‑86 (MQ=255)                                                                                                                                                                                            aTGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCCTGAACCGCcg                                                                                                                                                                                            >  2:765714‑M1/1‑86 (MQ=255)                                                                                                                                                                                            aTGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaac                                                                                                               >  2:596000‑M1/1‑86 (MQ=255)                                                                                                                                                                                              ggCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccacccggccctgcgcctcacatcctacccaaggtcccgttgccgatccgcgtca                                                                                     >  2:1184535‑M1/1‑84 (MQ=255)                                                                                                                                                                                                    gTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccac                                                                                                                                        >  2:1065312‑M1/1‑78 (MQ=255)                                                                                                                                                                                                         ttGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATCGCCACAGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaagcccgcacc                                                                                                                                                     >  1:947605‑M1/1‑73 (MQ=255)                                                                                                                                                                                                              ggAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGACTGAAcc                                                                                                                                                                                                >  2:501920/1‑66 (MQ=255)                                                                                                                                                                                                                aGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgc                    <  1:8787‑M1/236‑171 (MQ=255)                                                                                                                                                                                                                aGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgc                    <  1:943136‑M1/236‑171 (MQ=255)                                                                                                                                                                                                                  cGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgcccacccacctgcaccg                     >  2:829338‑M1/1‑64 (MQ=255)                                                                                                                                                                                                                  cGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgcggcgcagccctatcgcgacacaaggccgctcccacaaagaacgaagccatcagcat                                                                                                                             >  2:1511057‑M2/1‑64 (MQ=255)                                                                                                                                                                                                                     ccTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgc                 >  2:1024725‑M1/1‑61 (MQ=255)                                                                                                                                                                                                                       ttGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccacaaagaacgaagccatcagcatcacaacag                                                                                                                     <  2:431248‑M2/132‑74 (MQ=255)                                                                                                                                                                                                                       ttGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccacaaagaacgaagccatcagcatcacaacag                                                                                                                     >  1:431250‑M2/1‑59 (MQ=255)                                                                                                                                                                                                                       ttGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccac                                                                                                                                                  >  2:386682‑M1/1‑59 (MQ=255)                                                                                                                                                                                                                          tgtCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgccagctcc          <  1:829340‑M1/236‑181 (MQ=255)                                                                                                                                                                                                                            tCGCGAAAGGGGCGCGAAGCTGCCCGAGGCTTGCAGTGCCACTGCATGAACCGCcggggctgctgcgcaggcctagcgcgacacaa                                                                                                                                                              <  2:33721‑M1/86‑33 (MQ=255)                                                                                                                                                                                                                            tCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactcccgacccc                                                                                                                                        >  2:483906‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                             cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacatggccgctcccactccagaccaccaggccctgcgcc                                                                                                                           >  2:797620‑M1/1‑53 (MQ=255)                                                                                                                                                                                                                             cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgcca               <  1:1511409‑M1/228‑176 (MQ=255)                                                                                                                                                                                                                             cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctg                  >  2:507788‑M1/1‑53 (MQ=255)                                                                                                                                                                                                                             cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggt                                                    <  2:1491864‑M1/191‑139 (MQ=255)                                                                                                                                                                                                                             cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggt                                                    >  1:1491867‑M1/1‑53 (MQ=255)                                                                                                                                                                                                                             cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcc                                                                                                                           <  1:797622‑M1/120‑68 (MQ=255)                                                                                                                                                                                                                               cgAAAGGGGCGCGAAGCGGCCCCGGGTTTTCATTGCCACTGCATGACCCGCcggggctgctgcgcagccctatcgcgacacaa                                                                                                                                                              <  2:1394981‑M1/83‑33 (MQ=255)                                                                                                                                                                                                                                 aaaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgccagctccaacctg    <  1:988748‑M1/235‑187 (MQ=255)                                                                                                                                                                                                                                  aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtccggttgcagatcc                                                                                           >  2:227855‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                  aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgccagctccaacctgcc  >  1:1225298‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                  aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggtacccgcgccacctgctcacccacctgcaccgctgccagctccaacctg    >  2:1271642‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                  aaGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagccttt                                                                        >  2:946395‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                                         gcgAAGCGGCCCCAGGCTTTCATTGCCCCTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgc                                                                 <  1:293977‑M1/166‑126 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  CCGCCGTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGC  >  AE015451/5681143‑5681603

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