BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,782,920	G→C	14.2%	G247G (GGC→GGG)	gpsA ←	glycerol‑3‑phosphate dehydrogenase (NAD+)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,782,920	0	G	C	14.2%	201.1 / 6.8	105	G247G (GGC→GGG)	gpsA	glycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand): ref base G (39/50); new base C (10/6); total (49/56)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 106 total aligned reads displayed.

TCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGC  >  NC_000913/3782685‑3783154

tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGAGACTGGTTGTCGGTACAGGTAAGCACCAGAACGccc                                                                                                                                                                                                                                                >  2:963021/1‑238 (MQ=255)

tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑GTACCCGGCCGGTTGTCGGTACCGGTAAGCACCAAAACGccc                                                                                                                                                                                                                                                >  2:979578/1‑238 (MQ=255)

    tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGG‑TTACGCGACTGGTTGTTGGTACAGGTAAGCACCAAATACCC‑aacc                                                                                                                                                                                                                                            >  1:1284619/1‑236 (MQ=255)

    tgctgcCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGGTTTCGGTACAGGGAAGGCCCAAAACGCC‑CAGc                                                                                                                                                                                                                                            >  2:424139/1‑238 (MQ=255)

         ccTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTAAGCCACTGGTTGTCGGTGCCGGGAAGCACCAGAACGCC‑AAGCCCCGc                                                                                                                                                                                                                                       >  1:118653/1‑238 (MQ=255)

          cTCGCGCTCGTTTTTTCCGCAATATAATACTTGATAATGTTCCTCTGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGCCTTCTTTCGTATTGCGGTCGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACGGGTAAGCACCAGATCGCA‑AAGCCCCGCCa                                                                                                                                                                                                                                     <  2:34697/239‑1 (MQ=255)

           ccgcgcgcgTTTTTTCCGAAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAGGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACACCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCAt                                                                                                                                                                                                                                    <  2:1253600/238‑1 (MQ=255)

                          ccGAAATATAATCCTTAATAAATTTCCCCGGTTATGGCCTTTTCAACGCCGAAGGGATGCGCCGTTTCGCGGATTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTGCCCAATCTCCTCCTGCGCCCTTTGTACATCCATCCCCTGACCGAGCATCATGCAAAAACGGCGG‑TTCCGCGCCTGGTTGTCGGTACAGGTAAGCCCCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGc                                                                                                                                                                                                                     <  2:1378922/239‑1 (MQ=255)

                              aaTATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATa                                                                                                                                                                                                                             <  1:1016856/227‑1 (MQ=255)

                                tataATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCCCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGACCAAACTTCTCCGGCGCTCTTTGTACATCCATGCCCTGACCGAGCATCAGGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCAt                                                                                                                                                                                                                                    <  1:1389937/218‑1 (MQ=255)

                                     taCTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACCCGCCTGGGTGTCGGGGCAGGTAAGCACCAAAACGCCAAACCCCCGCCA‑GCCCCCCAAGGGGGCAGGggcggcc                                                                                                                                                                                                           >  2:1152983/1‑237 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCATCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTACACACCGGGTTTTGGGCAC‑GGTAAGCACCAAATCCCC‑AAACCCCCCCCCGCCCATAAAAAGGGCAgggggagcgcccca                                                                                                                                                                                                      >  2:248159/1‑227 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCTTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGCCTGGGTGTCGGGACAGGTAAGCACCCAAACCCC‑CAACCCCCCCCCGCCCACAAAAGTGGCAGGGGCGGCAccccc                                                                                                                                                                                                      >  2:1204478/1‑236 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTACAC‑CCTGGTTTTCTGTACAAGTAAGCACCAAAAAGGC‑CAACCCCCCCACACCCATaaaagctgtctctttaacacaactc                                                                                                                                                                                                     >  2:1389937/1‑217 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCTCCTGGTTGTCGGTACAGGGAAGAACCAGATCCCC‑AAAACCCCCCCCGCCCCTAAAAGGGGCAGGGTCGGGACCCac                                                                                                                                                                                                      >  1:39806/1‑237 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGCATGGTTGTCGGGACAGGGAAACACCAGAACCCC‑AAACCCCGCCATGCCCATAAAGGTGGCAGGGGCGGCAccccc                                                                                                                                                                                                      >  1:1135680/1‑236 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTGCAGGTAAGCACCAAAAAGCC‑AAGCCCCGCCCTGCCCCTAAAAGGGGGAGGGGCGGGAccccg                                                                                                                                                                                                      >  2:487356/1‑236 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGGACAGGGAAACAACAAAACGCC‑AAACCCCCCCCATCCCCTAAAGGTGGCCAGGGCGGCAACCac                                                                                                                                                                                                      >  1:1409045/1‑237 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGGTTTCGGGACAGGGAAGCACCAAAACCCC‑CAACCCCCCCCTTCCCCTAAAGGGGGCAGGGGCgggcccccc                                                                                                                                                                                                      >  2:769259/1‑231 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACCGGTTTTCGGTACAGGTAAGCACCAGATCGCC‑AAACCCCCCCATGCCCCTAAAAGTGGCAGGGTCGGCACCCAg                                                                                                                                                                                                      >  1:203874/1‑238 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCCGCCGGTTTTCGGTACAGGTAAGCACCAGAACGCC‑CACCCCCGCCCTGCCCCTAAAGGTGGCAGGGGCGGCACCCCgc                                                                                                                                                                                                     >  2:378119/1‑239 (MQ=255)

                                          gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTAAGCGCCCGGTTGTCGGTACAAGTAAAACCCACAACCCC‑CAACCCCCCCCTTCCCCTAAAGAGGGCCGGGGCGGCAccccg                                                                                                                                                                                                      >  2:1067617/1‑236 (MQ=255)

                                            tAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTTCGCGCCTGGTTGTCGGTACAGGGAAGCACCAAATCCCC‑AAGCCCCGCCCTTCCCATAAAGGGGGCAGGGTCGGCACCCAgcc                                                                                                                                                                                                    >  1:932088/1‑237 (MQ=255)

                                                      cGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGa                                                                                                                                                                                         <  1:938984/239‑1 (MQ=255)

                                                         ttATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCACCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑cg                                                                                                                                                                                      <  1:419490/238‑1 (MQ=255)

                                                          tATTGGCATTTCAACGCCGGGGCGTTGCGCCGTTTCGCGTTCTTTTTTCGTATTGCGGTAGCCTCCCACCCCCGGGCCAATCTTCTCCGGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCAGGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑cgcg                                                                                                                                                                                    <  2:966754/239‑1 (MQ=255)

                                                          tATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑cgcg                                                                                                                                                                                    <  1:310616/239‑1 (MQ=255)

                                                                   ttCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCGTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCAt                                                                                                                                                                                                                                    <  1:85238/183‑1 (MQ=255)

                                                                   ttCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGGAAGCACCAAATCGCC‑CAGCCCCCCCATGCCCCTAAAAGTGGCAGGGGCGGCACCCAGCCCCGCACCAAAA‑CGCGACATTTc‑a                                                                                                                                                                           >  2:581323/1‑238 (MQ=255)

                                                                   ttCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTCTCGGTACAGGGAAGCACCAGATCGCC‑CAACCCCGCCATGCCCATAAAGGTGGCAGGGTCCGCACCCAGCGCCCCACCaaacaccgacaattct                                                                                                                                                                             >  1:1330029/1‑224 (MQ=255)

                                                                   ttCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGGTGTCGGGACAGGTAAGCACCAAAACGGC‑AAGCCCCCCCAAGCCCCTAAAGGGGGCAGGGGCGGCACCCCGCGCCCCACCCAGA‑CGCGACAATTc‑a                                                                                                                                                                           >  1:1217021/1‑238 (MQ=255)

                                                                   ttCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCCACCGGTTGTCGGGACAAGGAAGCACCAAATCCCC‑CAGCCCCCCCCTGCCCCTAAAAGTGGCCGGGTCGGCACCCAGCGCCGCCCCAAGA‑CGCGACATTTcc                                                                                                                                                                            >  2:725121/1‑237 (MQ=255)

                                                                    tCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTCCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTc‑agc                                                                                                                                                                         <  1:424139/239‑1 (MQ=255)

                                                                    gcAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTc‑agc                                                                                                                                                                         <  1:487356/238‑1 (MQ=255)

                                                                    caaaCGCGGAAGCGATGCGCCAGTTCGCGGCCTTTTTTCGTATTGCGGTCGCCTTCCACCCCCTGACAAATCTTCCCCTGCGCGCTTTGTACTTCCATGCCCGGACCGAGCATCATTCCAAAACGGCGG‑TTACGCGCCTGGTTGTCGGTAAAGGTAAGCACCAGATCGCC‑AGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTc‑agc                                                                                                                                                                         <  2:1015847/237‑1 (MQ=255)

                                                                           cGAAGCGATGCGCAGGTGCGCGGTCTTCTTTCGTATTGCGGTAGCCTCCCCCCACCTTCCCAACCTTCCCCGCCGCGCTTTGTACATCCATGCCCTAGCCGGGACTCAGGCCAAACCGGCGT‑TTACGCGATTGTTTTTCGGTACAGGTAACCCCCAGACCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCa                                                                                                                                                                  <  2:118653/239‑1 (MQ=255)

                                                                            gAAGCGATGCGCCCGTTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTGACCAATCTTCCCCTGCGCGTTTTGTACATCCATGCCCTGCCCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCAc                                                                                                                                                                 <  2:932088/239‑1 (MQ=255)

                                                                             aaGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACCCGACTGGGTTTCGGGACAGGCAAGGACCCGACCGCC‑AAGCCCCGCCCTGCCCCTacaa                                                                                                                                                                                                                          >  2:85238/1‑180 (MQ=255)

                                                                              aGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACgg                                                                                                                                                               <  1:63073/239‑1 (MQ=255)

                                                                                      gcCAGTCCGCGGTCTTTTTTCGTATTGCGGTTGCTTCCCCCCACCGTACCAACCTTCTCCGCGGCGCTTTGTACAGCCATGCCCTGACCGAGCATCATGCAAAAACGGCGG‑TTACGCGCGTGGTTGTCGGTACAGGAAAGCACCAGACCGCC‑AGGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAg                                                                                                                                                       <  2:1217021/239‑1 (MQ=255)

                                                                                        cagtcgccgGTCTTTTTTCGTATGGCGGTACCCTTCCACCCCCTGACCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGTTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAgcg                                                                                                                                                     <  2:451949/232‑1 (MQ=255)

                                                                                             cgcgGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGGTGTCGGTACAGGTAAGCACCAGATCGCC‑AAACCCCGCCACGCCCCTAAAGGCGGCAGGGGCGGCACCCCCCGCCGCCCCCAGA‑CCCCACATTTC‑AGCCAGCCCACGCGTGGTCAGCGCCgc                                                                                                                                                 >  2:463937/1‑237 (MQ=255)

                                                                                                           gTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCCCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCAcc                                                                                                                                  <  1:248159/239‑1 (MQ=255)

                                                                                                            tATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACCGGGTTTCGGGTCAGGGAAGGCCCAGAACGCC‑CAAACCCCCCACGCCCCTAAAAGGGGGAGGGGCGGCACCCCGCCCCGACCCCCAA‑CCCGGAATTTT‑AGCCACGCCACAGGGGGTCAGCGcccccagcgccctcgcccca                                                                                                                                 >  2:910755/1‑221 (MQ=255)

                                                                                                              ttGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGGCTTGTTGTCGGGACAGGTAAGCACCAGAACGCC‑AAGCCCCGCCCTGCCCATAAAAGTGGCAGGGGCGGCCCCCAGCGCCGCCCCAAGACCGCG‑CATTTC‑AGCCCGCCCCCCGGGGCACAGCGCCGTACGCGCGTTCGCCCCaca                                                                                                                               >  2:464782/1‑237 (MQ=255)

                                                                                                              ttGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCaaa                                                                                                                               <  1:102266/239‑1 (MQ=255)

                                                                                                                gCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACCCGACCGGTTGTTGGTGCAGGGAAACACCAGAACGCC‑CAACCCCCCCCTGCCCCTAAAGGTGGCAGGGGCGGCACCCACCCCCGCCCCCAGA‑CGCCACATTTC‑CCCCAGCCCCCCGGGGAGCCGCGCCGCACGCGCACTCGCCCCaaaa                                                                                                                              >  2:646801/1‑238 (MQ=255)

                                                                                                                    tacccTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGa                                                                                                                          <  1:769259/235‑1 (MQ=255)

                                                                                                                    tAGCCTTCCCCCACCTGACCAATCTTCCCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAACCGGCGG‑TTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGa                                                                                                                          <  1:988602/238‑1 (MQ=255)

                                                                                                                             cccaccTGACAAATCTCCTCCTGCGCGCTTGGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑a                                                                                                               <  1:378119/238‑1 (MQ=255)

                                                                                                                             ccacccTGACAAATCTTCTCCTGCGCGCTTTGTCCACCCATGCCCGGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGCCGGTACAGGTAGGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑a                                                                                                               <  1:774673/235‑1 (MQ=255)

                                                                                                                             acccccGGACCAATCTTCTCCTGCGGGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGACCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑a                                                                                                               <  1:421624/239‑1 (MQ=255)

                                                                                                                             accaccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCAAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑a                                                                                                               <  1:1204478/239‑1 (MQ=255)

                                                                                                                              ccaccTGACAAATCTTCCCCTGCGCGCTTTGTACACCCATGCCCTGCCCGAGCAGCATGCCAAACCGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑Ac                                                                                                              <  2:203874/239‑1 (MQ=255)

                                                                                                                                   tgcccAATCTTCCCCTGCGCGCTTTGTACACCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATccc                                                                                                         <  1:939665/236‑1 (MQ=255)

                                                                                                                                   tgcccAATCTCCCCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGCCTGGTTGCCGGTACAGGTAAGCCCCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATccc                                                                                                         <  1:963021/236‑1 (MQ=255)

                                                                                                                                    caCCAATCTTCTCCGCCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGCCTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATcccc                                                                                                        <  2:224797/238‑1 (MQ=255)

                                                                                                                                              ccccTGCGCGCTTTGTACATCCATGCCCTGACGGAGCATCAGGCCAAACCGGCGG‑TTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATgg                                                                                             <  1:725121/237‑1 (MQ=255)

                                                                                                                                                    cgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGGCTGGTTGTGGGGTAAGGTAAAGCCCAGATAGCC‑AAACCCCCCCCTGCCCATAAAAGTGGGAGGGTAGGCCCCCCCCCCCGCCCCCAGA‑AGCGGCAATTT‑AGCCCCCCCACGCGTGAGCAGCCCCCGACCCCCTTTCCGACCCCAACCGCTACCACCGG‑CCACACCCGCACC‑AACGGCACTa                                                                                        >  2:910190/1‑238 (MQ=255)

                                                                                                                                                                 tCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTTACCCACCGGTTGTCGGTACAGGGAAGCACCAAAACGCC‑CAGCCCCGCCCCCCCCCTAAAGGCGGCAGGGGCGGCACCCCCCGCCGCCCCACAA‑AGCGCCACTTCAACCCACCCCACGGGGGA‑CAGCGCCCCACCCGCACTCGCCCCAAACACCATAACGTCCG‑CCACCCCCCCCAC‑CACCGCAAAA‑AACTTTTTCAccc                                                                          >  2:251512/1‑237 (MQ=255)

                                                                                                                                                                  ccAGGCCCTGACCGAGCATCAGGCCAAAACGGCGG‑TTACGCGACTGGTTGCCGGTCCAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATGCCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACcgcg                                                                        <  1:646801/239‑1 (MQ=255)

                                                                                                                                                                     tGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAAAACGCC‑CAGCCCCGCCATGCCCATACAGGCGGCAGGGGCGGCACCCAGCGCCGCCCCCAGA‑CGCGGCATTTC‑AGCCAGCCCACGGGGGGTCAGCGCCGCACGCGCATTCGCCCCAAAACCGATACCGCCGG‑ACAACCCCCCCCC‑CACACCAATA‑ACGATTTTCACCGcgcc                                                                      >  2:1386684/1‑238 (MQ=255)

                                                                                                                                                                     tGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACCCGACTGGGTGTCGGGAAAGGTAAACACCAAAACGCC‑AAACCCCGCCCTTCCCCAAAAAGGGGGAGGGGCGGGACCCAGCGCCGCACCCCGA‑CGCCACAATTC‑AGCCAGCCCACGGGGGGTCAGCGCCGTACCTGCCTTCTCACCAAAAACCATACCCTCGG‑ACAACCCCCCACC‑ACTGGCACTa‑aaactttttaacgcccc                                                                      >  2:666231/1‑222 (MQ=255)

                                                                                                                                                                      gCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAgcgc                                                                                                                                                                                                   <  2:1098848/117‑1 (MQ=255)

                                                                                                                                                                      ccccTGCCCGACCATCATCCAAAACCGGCGT‑TTCCGCGACTGTGTGGCGGTACGGCTAAGCACCAGATCCCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGcgccg                                                                     <  2:1284619/237‑1 (MQ=255)

                                                                                                                                                                            aCCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTTTCGGGACAGGGAAACAACAAATCGCC‑AAGCCCCCCCAAGCCCCTAAAAGGGGGAGGGGCGGCACCCCCCGCCGCACCCAGA‑CCCGCCAATTC‑AGCCAGCCCACCGGGGATCAGCGCCGCACGCCCACTCCCACCCAAACCGACACCCTCGG‑ACAACCCCCCACC‑ACAGGCAATA‑ACGATTTTCACCGCGCCGCCGAGc                                                               >  2:776819/1‑238 (MQ=255)

                                                                                                                                                                            aCCGAGCATCATGCCAAAACGGCGG‑GTACGCGACTGGTTGTCGGTACAGGTAAGCACCAAATCGCC‑AAACCCCCCCCAGCCCCAAAAAGAGGCAGGGGCGGACCCCCGCGCCCCCCCAAGA‑CGCGACAATTT‑ACCCAGCCCACCGGGGAGCCGCCCCCTACACCCCCTCCCACCCAAACCGATAACGGCGGCACA‑CCCCGCACCAAACGGCAACACACGTTTTCCCCCGccccccaac                                                                 >  2:696631/1‑233 (MQ=255)

                                                                                                                                                                             ccgaccagagtGCAAAAAGGGCGG‑TTACGCGACTGGTGGCGGGAACAGGAAACCACCAGACCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTg                                                             <  2:39806/229‑1 (MQ=255)

                                                                                                                                                                             ccGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGCCGGTACAGGTAAGCCCCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTg                                                             <  1:1152983/239‑1 (MQ=255)

                                                                                                                                                                             ccGAGCATCATACCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACGGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTg                                                             <  1:50866/239‑1 (MQ=255)

                                                                                                                                                                             ccGACCATCATGCAAAAACGGCGG‑TAACGCGACGGTTTGTCGGTAAAGGTAAGCACCAGACCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTg                                                             <  1:910190/239‑1 (MQ=255)

                                                                                                                                                                                aGCATCATGCCAAAACGGCGG‑TTACGCGACTGGGTGTCGGGACAGGGAAGCACCAGATCGCC‑AAACCCCGCCATGCCCATAAAGGGGGCAGGGGCGGCACCCAGCGCCGCAcccaaa                                                                                                                                                                                         >  1:1098848/1‑113 (MQ=255)

                                                                                                                                                                                   atcagGCCAAACCGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGc                                                        <  1:954834/238‑1 (MQ=255)

                                                                                                                                                                                      atGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGGACACGTAAGCAACAAAACCCC‑AAACCCCGCCATGCCCATAAAAGTGGGAGGGGTGGCACCCAGCGCCGCCCCAAAA‑ACCGACATTTT‑AGCCCCCCCCCCGGGGAAAAGCGCCGCACCCGTCTTCGCACCAAAACCGGTACCGTCCG‑ACATCCCCCCCCC‑AACGGCCATA‑ACGGTTTTCACCCCcccgcccacactcaccccaac                                                    >  2:799458/1‑224 (MQ=255)

                                                                                                                                                                                       tGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGGACAGGGAAGCACCAAATCGCC‑CAGCCCCGCCATGCCCATAAAGGGGGCAGGGGCGGCACCCCGCACCGCACCAAAA‑CGCGACATTTC‑AACCAGCCCCCGGGGGATCAGCGCCGCACGCGCATTCCCACCAAAACCGATACCGGCGG‑ACAACCCCCCACC‑CACGGACAAA‑ACGCTTTTTACCCCgccgccccactccacccccaag                                                   >  2:956506/1‑223 (MQ=255)

                                                                                                                                                                                           aaacCGGCGG‑TTCCGCGCCTGGTTGCCGGTACAGGTAAGCCCCAGCTCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAAt                                              <  1:955100/239‑1 (MQ=255)

                                                                                                                                                                                                 gcgg‑TTACGCGACTGGTCGTCGGTACAGGTAAGCACCAGACCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGa                                         <  1:696631/238‑1 (MQ=255)

                                                                                                                                                                                                  ggg‑TTACGCCCCTGGTTGTCGGTACAGGTAAGCACCAGACCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGAt                                        <  1:1386684/237‑1 (MQ=255)

                                                                                                                                                                                                  cgg‑TTACGCGACTGGTTGTCGGGACAGGGAAGCACCAAAACGCC‑AAGCCCCGCCATGCCCATAAAGGGGGCCGGGGCGGCACCCAGCGCCGCCCCCAAA‑CGCCACATTTC‑AGCCAGCCCACGGGGGATCAGCGCCGGACGCGCATTCGCACCAAAAACGAAACCGCCGG‑ACATCCCCGCACC‑AAAGGCAATA‑ACGTTTTTCACCCCCCCCCCCAACCGCACGCCAATCGAAACCGGGAt                                        >  2:342987/1‑239 (MQ=255)

                                                                                                                                                                                                         cgcgACTGGTTGGCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTa                                 <  1:334791/239‑1 (MQ=255)

                                                                                                                                                                                                             aCTGGTTGTCGGTACAGGTAAGCACCAGACCCCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑Gatctccagcagctgctgcac                              <  1:463937/238‑20 (MQ=255)

                                                                                                                                                                                                               tgtttGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAAcgcg                           <  1:666231/236‑1 (MQ=255)

                                                                                                                                                                                                                ggTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCgg                          <  1:563182/239‑1 (MQ=255)

                                                                                                                                                                                                                          aCAGGGACGCCCCAGACCCCC‑AAGCCCCGCCAGGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGc                <  2:1108464/239‑1 (MQ=255)

                                                                                                                                                                                                                           cAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGCGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AAAGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAACATTTTGcc               >  2:225571/1‑239 (MQ=255)

                                                                                                                                                                                                                              gTAAGCCCCAGACCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGc             <  1:559307/238‑1 (MQ=255)

                                                                                                                                                                                                                              gTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATaa                                                                                                                                                                                                                            <  2:273188/37‑1 (MQ=255)

                                                                                                                                                                                                                               caaGCTCCAGCCCGCA‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCg              <  1:994600/235‑1 (MQ=255)

                                                                                                                                                                                                                                   caccacACCCCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca       <  2:1135680/239‑1 (MQ=255)

                                                                                                                                                                                                                                   cACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGTCACCCCACGGGGGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCCCGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTga        >  2:897295/1‑237 (MQ=255)

                                                                                                                                                                                                                                   cACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCCGTgca       >  1:1109619/1‑239 (MQ=255)

                                                                                                                                                                                                                                   cACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca       >  1:1195779/1‑239 (MQ=255)

                                                                                                                                                                                                                                   cACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAACCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca       >  1:368926/1‑239 (MQ=255)

                                                                                                                                                                                                                                    cccAGATCCCC‑ACGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgca       <  1:1067617/237‑1 (MQ=255)

                                                                                                                                                                                                                                     cccgacagCC‑ATGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcagc     <  1:494667/232‑1 (MQ=255)

                                                                                                                                                                                                                                      cAGACCCCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcagca    <  1:799458/239‑1 (MQ=255)

                                                                                                                                                                                                                                        gATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAgg                                                                                                                                                                                                                  >  1:273188/1‑37 (MQ=255)

                                                                                                                                                                                                                                        gATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcagcacc  >  2:893916/1‑237 (MQ=255)

                                                                                                                                                                                                                                        gATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcagca    >  2:1152932/1‑237 (MQ=255)

                                                                                                                                                                                                                                        gATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTgcagca    >  1:509444/1‑237 (MQ=255)

TCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGG‑TTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCC‑AAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGA‑CGCGACATTTC‑AGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGG‑ACATCCCCGCACC‑AATGGCAATA‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAAT‑GAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGC  >  NC_000913/3782685‑3783154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`