BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	5,681,373	G→C	5.0%	intergenic (‑91/+113)	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,373	0	G	C	5.0%	346.3 / 7.9	99	intergenic (‑91/+113)	PP_4986/PP_4987	Channel protein, hemolysin III family/putative Chemotaxis protein
Reads supporting (aligned to +/- strand): ref base G (54/40); new base C (2/3); total (56/43)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 103 total aligned reads displayed.

TGCCGCCGTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGC  >  AE015451/5681141‑5681603

tGCCGCCGTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGgcgc                                                                                                                                                                                                                                     >  1:902459/1‑236 (MQ=255)

  ccgccgTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGcgcg                                                                                                                                                                                                                                    <  1:243747/235‑1 (MQ=255)

  ccgccgTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGa                                                                                                                                                                                                                                   >  1:961781/1‑236 (MQ=255)

   cgccgTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGgcgc                                                                                                                                                                                                                                     >  2:1300694/1‑233 (MQ=255)

                gaAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGc                                                                                                                                                                                                                      <  1:1300697/235‑1 (MQ=255)

                  aGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCt                                                                                                                                                                                                                     >  1:964616/1‑234 (MQ=255)

                     tGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGgcgc                                                                                                                                                                                                                                     >  1:1139636/1‑215 (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)

                                                                                                        aGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAATGGgcg                                                                                                                                                                                                                                      >  2:746026/1‑131 (MQ=255)

                                                                                                          gTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgc                                                                                                                            <  1:808237‑M1/235‑67 (MQ=255)

                                                                                                             ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggcc                                                                                                                                  <  1:957894‑M1/226‑61 (MQ=255)

                                                                                                             ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggcc                                                                                                                                  >  2:957892‑M1/1‑166 (MQ=255)

                                                                                                             ggTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactc                                                                                                                                                >  2:687389‑M1/1‑166 (MQ=255)

                                                                                                               tCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcaca                                                                                                                      <  1:388425‑M1/236‑73 (MQ=255)

                                                                                                                  aGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCAtt                                                                                                                                                                                                               >  1:1017332/1‑144 (MQ=255)

                                                                                                                  aGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCAtt                                                                                                                                                                                                               <  2:1017329/144‑1 (MQ=255)

                                                                                                                  aGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggct                                                                                                                                                                                       >  2:1511406‑M1/1‑161 (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)

                                                                                                                                                               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                                                                                                                     >  1:431250‑M2/1‑59 (MQ=255)

                                                                                                                                                                                                                        ttGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccacaaagaacgaagccatcagcatcacaacag                                                                                                                     <  2:431248‑M2/132‑74 (MQ=255)

                                                                                                                                                                                                                        ttGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccac                                                                                                                                                  >  2:386682‑M1/1‑59 (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                                                    >  1:1491867‑M1/1‑53 (MQ=255)

                                                                                                                                                                                                                              cgcgAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCcggggctgctgcgcagccctatcgcgacacaaggccgctcccactccagaccaccaggccctgcgcctcacatcctaaccaaggtcaggttgcagatccgcgtcacgcaccagcctttccaatgcgaccaggtccggt                                                    <  2:1491864‑M1/191‑139 (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)

TGCCGCCGTAGATGGAGAAGCTGACGATCTTCCACGGGTCGCCTTGCAGGCCCGCGGCCACGATCAGCCAGATGGCACCAATAGCAGCGAGGACGGCACCGACCAGGTGGGTCCAGGCGTTGAAGCGTTCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCACCTGCTCACCCACCTGCACCGCTGCCAGCTCCAACCTGCCCAGC  >  AE015451/5681141‑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: `‑`