BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,137,349	T→C	10.4%	E103G (GAG→GGG)	yijF ←	DUF1287 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,349	0	T	C	10.4%	291.8 / 3.0	117	E103G (GAG→GGG)	yijF	DUF1287 family protein
Reads supporting (aligned to +/- strand): ref base T (53/47); new base C (5/9); total (58/59)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGT  >  NC_000913/4137111‑4137578

gcaaaTTTTATGTATCACCAGCGGCGTACCGTCCCGGGCAAAGCCATCTGAAACCACCCCAATATGTCCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCCCCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACt                                                                                                                                                                                                                                        <  1:15049/236‑1 (MQ=255)

gcCAAGGTTAGGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACt                                                                                                                                                                                                                                        <  1:594968/239‑1 (MQ=255)

   aaTGTTATGGAACACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAg                                                                                                                                                                                                                                     <  1:173070/239‑1 (MQ=255)

      gTTATGTATCCCCACGGGCGTACCGTCGCGGGCGAAGCCAGCTGAAACCCCCCCAATATGTGCAAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTGGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCaa                                                                                                                                                                                                                                  <  2:427147/239‑1 (MQ=255)

          tGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCCGG‑TTTCCCGATTAGGCACCCGACGGGGATCGATGTTGCT‑TGCCGGGCGCTTTA‑ACTTCC‑ATTTTTTCGGGGACTCACGAAAAtt                                                                                                                                                                                                                              >  1:879060/1‑239 (MQ=255)

          tGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCCAAACCAGG‑GTTCCAGGTTTGGAGCCCCGCGCGGATCGCGTTTTCT‑GTCCGGGGGCTTTT‑AAACCC‑ACTTTTTCGGGTACCCAGCaaaaat                                                                                                                                                                                                                              >  2:514399/1‑237 (MQ=255)

           gTATCCCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGGGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAAttc                                                                                                                                                                                                                             <  1:556764/239‑1 (MQ=255)

               ccccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATTTCGCCCGCTTGATAGTCCCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCg                                                                                                                                                                                                                         <  1:495510/237‑1 (MQ=255)

               cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCCGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCGTTA‑ACTTCC‑ATTTTTTCGGGTACTCAGCAAAAttcttc                                                                                                                                                                                                                          >  1:732435/1‑238 (MQ=255)

                 ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCCCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCCGG‑GTTTCCGCTTTGGGACCCCCCGGGGGTCGAT‑TTGTTCTTCCCGGCGCCTTT‑AATTCC‑CTTTTTTCGGGGACCTACCAAAAttttttccc                                                                                                                                                                                                                       >  2:641040/1‑235 (MQ=255)

                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCCTGGCGGCCTAACCAGG‑CTTTCAGGTTCTGCGCCCGACGGTGATCGGTGGTGAC‑GTCCCGGCGCTTTACATTTGC‑TTTTTTG‑GGGTAGTCCGCAAAAATTTACGcccac                                                                                                                                                                                                                    >  2:195928/1‑236 (MQ=255)

                        cgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCCGAATAGGCACCCGACGGGGACTGAGGTTGCT‑GCCCCGGCGCTTTA‑AATTCC‑CATTTTTCGGGGAATCCGCCAAATTCTTCGccctctttt                                                                                                                                                                                                                >  2:247498/1‑232 (MQ=255)

                        cctACCCTCGGGGGCAGAGCCATCTTAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCt                                                                                                                                                                                                                 <  2:415477/236‑1 (MQ=255)

                                         aGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGc                                                                                                                                                                                              >  1:257105/1‑239 (MQ=255)

                                                   aaCCCCCCCAAGAGGTGCCAGCCCGTTATCCAGTCGCCAGGAGCCAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCt                                                                                                                                                                                    <  1:166085/239‑1 (MQ=255)

                                                        ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAACAGGTTTTCCAG‑TTAGGGACCCGGCGGTGATCGCTGGTGCT‑GTCCGGGCGCTTTT‑ACTTCA‑CTTTTTGCGGGTACTCAGGAAAATTCTTCGCCTAATTCTCTTGAACCAG‑GTTCTGCAAATCCACCTTTTgc                                                                                                                                                                               >  2:361775/1‑238 (MQ=255)

                                                         cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTg                                                                                                                                                                                >  1:94582/1‑237 (MQ=255)

                                                               ttggccccgCCCGTTTTCCAGTCGCCAGGAGAAAATATCGCCCGGTTGATAGTCACGGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCa                                                                                                                                                                        <  2:127335/231‑1 (MQ=255)

                                                                   gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTCTCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGGACCCGACGGTGATCGATGTTGCT‑GTCCCGGCGCCTTA‑ACTTCC‑ATTTTTGCGGGTACTCCGCAAAATTTTTCGCCATATCCTCGTGAACCAA‑TTTCTGCAAATCGACCTTCCggacgcacaaagc                                                                                                                                                                    >  2:443062/1‑228 (MQ=255)

                                                                   gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGGACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTTGCTGCGCAAta                                                                                                                                                                     >  1:127335/1‑237 (MQ=255)

                                                                   gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGAACCCCACGGTGATCGAGGTTGCT‑GTCCGGGCGCTTTA‑TCTTCC‑CTTTTTGCGGGGACCTGGCCAAATTCTTTGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCCTCTGGCTGCGCAAtc                                                                                                                                                                     >  2:772777/1‑237 (MQ=255)

                                                                   gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCCGCTAAACCAGG‑GTTCCAGATTAGGCTCCCGACGCGGATCGTTGTTCCT‑GTCCTGGGGCTTTT‑CCTTCC‑CTTTTTTCGGGTACCCACCAAAAATCTTTGCCATACTTTTTTGAACCAG‑TTTTTGCACAtacacactctcgcgtcgcaattg                                                                                                                                                                    >  2:373625/1‑217 (MQ=255)

                                                                    ccacccctgtttCCATTCCCCGGGAGACAATTTCGCCCGTTTGATAGTCCCGGGGGTTTTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATgcg                                                                                                                                                                   <  2:94582/228‑1 (MQ=255)

                                                                    ccAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAACAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTT‑ACTTCC‑ATTTTTGCGGGGACTTAGCAAAATTTTTCGCCCTATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATgcg                                                                                                                                                                   >  1:922489/1‑239 (MQ=255)

                                                                     cAGCCCGTTATCCGGGCGCGAGGGACCAATATCGCCCGGTTGATGGTCACGGGGGTTCTTGCTGGTGGGGCGGGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGGTTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATgcg                                                                                                                                                                   <  2:879060/238‑1 (MQ=255)

                                                                       gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCCGG‑TTTTCCGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTTCCGGTACTCAGCAAAATTCTTCGCCATATCTTCCTGAACCCG‑TTTTTGCAAAACCACCTTCTCGGTGCCCAAATCGCgg                                                                                                                                                                >  2:386133/1‑239 (MQ=255)

                                                                       gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTTGGCACCCGACGGTGATCGATGTTGGT‑GTTCGGGCGCTTTA‑ACTTCC‑ATTTTTTGGGGGACTCCTCCAAATTCTTCGCCATTACTTCGCGAACCCG‑TTTCTGCAAAATAACCTTCGTGCCCCGCAATGCGAgg                                                                                                                                                                >  2:12762/1‑239 (MQ=255)

                                                                         ccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGAt                                                                                                                                                              <  1:854689/239‑1 (MQ=255)

                                                                            ttATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCCAAAACAGG‑GTTCCAGATTAGGCACCCGACGGGGAACGATGGTTCT‑GGCCGGGGGCCTTA‑ACTTCCAATTTTTGC‑GGTAATCACCACAAATCTTCGCCAGATATTCGTTAACCAG‑TTTATGCAAAACGACCTTTTGGCTGCCCAATTCGCGGaaaa                                                                                                                                                            >  2:734699/1‑235 (MQ=255)

                                                                                cccgTCGCCAGGGGAAAATTTCGCCCGCTGGATGGCCCCTGGGGTTTTTGGTGGTGGGGCTTGTTTTACCATGGGGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACcaca                                                                                                                                                       <  2:1048369/236‑1 (MQ=255)

                                                                                    tccccctgtgaaattttcccccGTTTGATACGCCATGGGGTTCGTGGTGGTGGGGGCGGTTTTTACATGGCGGACAAACCCGG‑TTTCCAGGTTAGGCACCCGGCGGTGATCAATGTTGCC‑GTCCGGGCGCTTTA‑CCTTCC‑ATTTTGGCGGGTACGCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCGCCTCATCgg                                                                                                                                                   <  2:676670/220‑1 (MQ=255)

                                                                                    tCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCCTTA‑ACTTCC‑ATTTTTGCGGGGACTCAGCAAAATTCTTCGCCCTATCTTCGTGAACCAG‑TTTCTTGAAAACGACCTTCTGGCCGCGCAAAGCCCGGATacccaccacgg                                                                                                                                                   >  2:336852/1‑234 (MQ=255)

                                                                                     cgccggagcaCAATATCGCCCGCTGGATAGTCCCTGGGTTTCTTGCTGGGGGGGCGTTTTTTATCATGGGGGCAAAACCAGG‑TTTCCAAATTAGGCACCCGACGGTGATCGATTTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGCCCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGa                                                                                                                                                  <  1:652347/230‑1 (MQ=255)

                                                                                     cGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCgg                                                                                                                                                   <  1:809997/238‑1 (MQ=255)

                                                                                      gCCAGGGGGCAATTTCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGGTGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCCCCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGaa                                                                                                                                                 <  1:443062/239‑1 (MQ=255)

                                                                                        cAGGGGACAATTCCCCCCGCTTGATTGTCACTGGGGTTGTTTCTGGTGGGGCGTTTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCCCCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACa                                                                                                                                               <  2:581506/239‑1 (MQ=255)

                                                                                             aaCAATCTGGCCCGCATGTTAGTGCCGGGGGTTATGGTTGGTGGGCCGGTTTTTGTCATGGCGGCTAACCCTGG‑TTTCCTAGTTAGGAACCCGACGTGGACAGATGTTGCG‑GTCCTGTCGATTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTGTGGCTGCGCAATGCGCGGATCCCCACATCGGAACATaca                                                                                                                                           <  2:757484/237‑1 (MQ=255)

                                                                                                agatgggaataagagccagttccTGGGTTTCTTGCTGGTGGGGCGTGTTTTAGCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACcg                                                                                                                                        <  1:666149/216‑1 (MQ=255)

                                                                                                aaTATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTgcg                                                                                                                                                                          <  1:932457/204‑1 (MQ=255)

                                                                                                aaTATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCCACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGGAATCAGCAAAATTCTTCGCCCTATCTTCGTGGAACAGTTTTCTCCAAA‑CGACCTTCTGGCCGCGCAAAGCGGGGCACCCCCCATCGGAAAAAACAcccc                                                                                                                                       >  2:464141/1‑237 (MQ=255)

                                                                                                  tttctcccgctgtatagatccTGTGGTTCCTGCTGGGGGGGTGTGTTTTATCGTGGCGTGCAACCCAGG‑TCTCCAGTTTAGGCATCCGACGGTGAACGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGTGTACTCAGCAAAATTCTTCGCCATAGCTTCTTGAACCAG‑GTTCTTCAAATCGACCCTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGt                                                                                                                                     <  2:392224/219‑1 (MQ=255)

                                                                                                       cccGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGACGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTg                                                                                                                                <  1:506986/239‑1 (MQ=255)

                                                                                                          gCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGCCGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ACTTTTGCGGGGACTCAGCAAAATTTTTCGCCCCATCTTCGGGGACCCG‑TTTCTGCGAATCCACCTTCCGCCTCCGCCCACTCCGGACCACCACATCCGAACaacaacccccgtcttgggg                                                                                                                             >  2:606386/1‑221 (MQ=255)

                                                                                                            ttGATAGTCCCTGGTGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTg                                                                                                                                <  1:341490/234‑1 (MQ=255)

                                                                                                            ttGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑GTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGGACATAAACCCCGTTCTTTCGGaa                                                                                                                           >  1:721420/1‑239 (MQ=255)

                                                                                                                     cttGGGTTTTTTGTTGGTGGGGCTTTTTTTACCAGGGCGGTTAAACCAGG‑TTCCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGc                                                                                                                  <  2:257105/237‑1 (MQ=255)

                                                                                                                     ccgggggTTCTTGCTGGTGGGGCGTTTTTTATCAGGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCCCCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGc                                                                                                                  <  1:642536/236‑1 (MQ=255)

                                                                                                                     aCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATgcg                                                                                                                                                                   >  1:880716/1‑190 (MQ=255)

                                                                                                                     aCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATgcg                                                                                                                                                                   <  2:880716/190‑1 (MQ=255)

                                                                                                                      gtGGGTTTCTTGCGGGGGGGGCGTTTTTTATCAGGGCGGTAAACCCAGT‑TTTCCAGATTAGGCCCCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGAAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGc                                                                                                                  <  1:734699/237‑1 (MQ=255)

                                                                                                                           gTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCtt                                                                                                                                                                                   <  1:1037402/168‑1 (MQ=255)

                                                                                                                             tCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGGGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCCCATCGGAACATAAACCCCGTTCTTGCGGAACCTCACCCCCCggcaaa                                                                                                          >  1:144463/1‑235 (MQ=255)

                                                                                                                              cTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑CTTCCAGGTTAGGCACCCGCAGGGGGTCGGTGTTGCT‑GTTCGGGGGCTTTA‑ACTTCC‑ATTTTTGCGGGGACTCAGGCAAATTCTTTGCCATATCCTCGTGCGCCAG‑TTTCTTCAAATCGAACTTCTTGCTGCGCAATGCGCGGATCACCACCTCGGAACAAACACCGCGGTCTTggggggggctcttgtacaccaccg                                                                                                         >  2:666149/1‑219 (MQ=255)

                                                                                                                               ttttgtgtggGGGGTTTTTTATCAGGGCGGTAACACTGGG‑TTGCCAGTGTAGACACCCGACGGTGATCGATGTTGGT‑GTCCGGGCGATTTA‑CTTTCC‑ATTTTTGGGGGTACTCGGCAAAATCCTTCCCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCCCAATGCGCGGATCACCACATCGGAACATACACCGCTTTCTTGCGGAACATCACCGCCAGGAtaagt                                                                                                        <  2:163520/233‑1 (MQ=255)

                                                                                                                               tggctgtgggggCGTTTTTTACTGTGGCGGATAACCCAGG‑TTGCCAGATTAGGCCCCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGAtaagt                                                                                                        <  1:464141/231‑1 (MQ=255)

                                                                                                                                 gCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTT‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTTTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTTCGGCACATTCCCCCCAGCAtaagaaa                                                                                                      >  1:163520/1‑239 (MQ=255)

                                                                                                                                        gggCGTGTTTTATCATGGCGGCTAAAACAGG‑TTTTCAGATTAGGCACCCGACGGGGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑CTTTTTGCGGGTACTCAGCACAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCCAATCGACCTTCTGGCTGCGCACTGCGCGGATCACCACATCGGGACATACACCCCGTTCTTGCGGAACAACACCCCCCGGGtaagaaaagagcca                                                                                               >  2:935514/1‑232 (MQ=255)

                                                                                                                                             tgtTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTTGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCCCCTTCTTGCGGAACATCACCGCCCGGGTAAGTAAACTGCACATAc                                                                                           >  1:675553/1‑238 (MQ=255)

                                                                                                                                                tttagcagtgcgtgtAAACCAGG‑TTTCCAGTTTAGGCCCCCGACGGTAACCGATGTTGCT‑TTCGGGTCGCTTAA‑ACTTCC‑ATTTTTGCGGGTACTGAGCAAAATGCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGAATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCg                                                                                       <  1:245488/225‑1 (MQ=255)

                                                                                                                                                tttATCAGGGCGGCTAAACCGTG‑TTTCCAGATTAGGCACCCGACGGTGATCGATTTTGCT‑GTCCGGGCGCTTTA‑GCTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACCCCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCg                                                                                       <  1:577853/239‑1 (MQ=255)

                                                                                                                                                tttAGCAGGGCGGCTAAACCAGG‑TTTCCATATTAGGCCCCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTGGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCg                                                                                       <  1:641040/239‑1 (MQ=255)

                                                                                                                                                  ttgtcgggccgttAAACCGGG‑TTTCCATTTTGGCCCCCCGAGGGTTATCGGTGTTGCT‑TGCCGGGCGTTTTA‑ATTTAC‑ATTTTTGCGGTGACTCAGCAAAATTCTTCCCCATATCTTCGTGAACAAG‑TTTCTGCAAACCGCCCTTCTGGCTGCGCAATGCGCGGATCCCCACATCGGAACATACACCGCGTTCTGGCGGAAAATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGa                                                                                     <  2:153859/227‑1 (MQ=255)

                                                                                                                                                  tatctgggggttaaacacggg‑TTTCCAGTTTGGGCCCCCGTCGGTGCTGGAGGTTGCT‑GTCGGGGCGTTTTA‑ACTTCC‑ATTTTTGCGGGTACGCAGCAAAATTCTCCGCCATATCTTCGGGAACCAG‑TTTCTGCAAATCGCCCTTCTGGCTGCGCAATGCGCGGATCACAACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGa                                                                                     <  2:766310/220‑1 (MQ=255)

                                                                                                                                                  tatctggggggcTAACCCGGG‑TTTCCAGTTTAGGCACCCTACGGTGATCGATGTTGCG‑TTCCGGGCTTTTTA‑ACTTCC‑ATTTTTGCGGTTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGCCCTTCTGGCTGCGCAATGCGCGGATCCCCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGa                                                                                     <  2:144463/230‑1 (MQ=255)

                                                                                                                                                  gaacatggggtatAAACCAGT‑TTCCCAGATTGGGCCCCCGAGGGTGATCGTTGTTGCT‑GTCCGGGTGCTTTA‑ACTGCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGa                                                                                     <  1:373625/227‑1 (MQ=255)

                                                                                                                                                  atcctggggcgtaaaacaggg‑TTTCCGGATGGGGCCCCGGGCGGTGATCGATGTTGCT‑GCCCGGGTGCTTTA‑ACTTCC‑ATTTTTGCGGGTACCCAGCAAATGTCTCCGCCCTATCTTCGTGAACCAG‑TTCCTGCAAATCGACCTTCTGGCTGCGAAATGCGCGGATCCCCACATCGGAACATACACCGCTTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGa                                                                                     <  2:732435/220‑1 (MQ=255)

                                                                                                                                                   aTCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCCGGATAAGAAAGCTTCACCTAcccccgg                                                                                     >  1:425528/1‑235 (MQ=255)

                                                                                                                                                       tGGCGG‑TAAACCCGGTTTTCCAGATTGGGCCCCGGAGGGTGATCGATGTTGCG‑GTCGGGGCTCTTTA‑ACTTCC‑ATTTTGGCGGGTACCCAGAAAAATTCTCCGCCATATCTTCGTGAACCAG‑TTGCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTa                                                                                <  2:676379/239‑1 (MQ=255)

                                                                                                                                                        ggcggcTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATAAACCGCGTTCTTGCGGGACATCACCCCCCGGATAAGAAAACTGCCCATACGCCGGGTCGTaa                                                                               >  1:417464/1‑239 (MQ=255)

                                                                                                                                                        ggcggcTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATAAACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGACCATACGCCGGATCGTaa                                                                               >  1:1033089/1‑239 (MQ=255)

                                                                                                                                                         gcggcTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATTCGCGGATAACCACATCGGAACATACACCGCGTTCTTGCGGAACATCCACGCCAGGATAAGAAAGCTGCACCTACGCCGGATCGTaaa                                                                              >  1:81927/1‑239 (MQ=255)

                                                                                                                                                              aaaacgcgg‑TTCCCAGTTTAGCCCCCCGTGTGTTGTCGATGTTCCT‑GGCCGGGCGCTTTA‑ACTGCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCAGATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGc                                                                         <  2:294419/232‑1 (MQ=255)

                                                                                                                                                                      g‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGGTCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGCAAAGGAAGATGCACAAACGCCGAACGGTAAAAAAGCGGGCGGcc                                                                 >  2:528020/1‑239 (MQ=255)

                                                                                                                                                                             ggtttgggcgccgggcGGTGGGCGATGTTGCT‑GTCCGGGCGCTTTA‑TCGCCT‑ATTTTTGCGGGTATGCAGGAAAATTGTCCGCCATAGCGTCGGGAACCAG‑TTTCGGAAAATGAACCTTCAGGGTGCGCAATGCGCGGACCACCCCAGCGGAACATACGCCGTGTTCTGGCGGAACTACACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTg                                                           <  2:702944/224‑1 (MQ=255)

                                                                                                                                                                                ttAGGCCCCCGCGGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTCCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAAGCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTg                                                           <  1:247498/236‑1 (MQ=255)

                                                                                                                                                                                 tAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGc                                                                                                                  <  1:371992/180‑1 (MQ=255)

                                                                                                                                                                                 tAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGc                                                                                                                  >  2:371992/1‑180 (MQ=255)

                                                                                                                                                                                         cGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCCGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATATGCTGTCTGGCAcc                                               >  1:1078089/1‑239 (MQ=255)

                                                                                                                                                                                                     aTGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGACGGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGc                                   >  1:1002605/1‑239 (MQ=255)

                                                                                                                                                                                                      tttttccttcccggggcttttA‑CCTCCC‑ATTTTGGGGGGTCCCCAAAAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAACCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGc                                   <  2:1033089/221‑1 (MQ=255)

                                                                                                                                                                                                      tggtGCT‑GTCCGGGCCCTTTA‑ACTTCC‑ATTTTTGGGGGTACCCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGCCCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGc                                   <  1:125220/235‑1 (MQ=255)

                                                                                                                                                                                                      tGTTGGT‑GCCCGGGCGCTTAA‑ACTTCC‑ATTTTTGCGGGTTCTAAGCAAAATTCTTCGCCATATCCTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCCCCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGc                                   <  1:1086015/238‑1 (MQ=255)

                                                                                                                                                                                                      tGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTGGCGGAACATCACCGCCCGGCTAAGTAAGCTGAACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGATGTCGGGCACCGTCTGCGATGGCa                                  >  1:226352/1‑239 (MQ=255)

                                                                                                                                                                                                         tGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGTTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCACTTTGCTGTCTGGCACCGTCTGCGATGGCCAgg                               >  1:162790/1‑239 (MQ=255)

                                                                                                                                                                                                         tGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTCCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTGACTGACTGGCACCGTCTGCGATGGCAAgg                               >  1:195065/1‑239 (MQ=255)

                                                                                                                                                                                                          gCT‑GGCCGGGCGTTTTA‑ACTTCA‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGCCCTTCTGGCTGCGCAATGCGCGGATCACCACACCGGAACATACACCGCGTTCTTGCGGAACACCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGt                              <  2:122653/239‑1 (MQ=255)

                                                                                                                                                                                                           cT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGTTGGCAAGGtt                             >  1:61770/1‑239 (MQ=255)

                                                                                                                                                                                                               tCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGAACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAGGGTTGGt                          >  1:322363/1‑239 (MQ=255)

                                                                                                                                                                                                                 cGGGCTCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATc                                                                                                                                                                                         <  2:443996/78‑1 (MQ=255)

                                                                                                                                                                                                                 cGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATc                                                                                                                                                                                         >  1:443996/1‑78 (MQ=255)

                                                                                                                                                                                                                 cGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTAt                        >  1:1038341/1‑239 (MQ=255)

                                                                                                                                                                                                                 cGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGAGGGCAAGGTTGGTAt                        >  1:373305/1‑239 (MQ=255)

                                                                                                                                                                                                                    gcgcTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCCTCTGCGATGGCAAGGTTGGTATTg                      >  2:96801/1‑238 (MQ=255)

                                                                                                                                                                                                                      gcTTTA‑ACTTCC‑TTTTTTGCGGGTACCCAGCAAAATTCTTCGCCATATCTCCGTGAACCAG‑TTCCTGCAAATCGGCCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGcc                   <  1:225099/239‑1 (MQ=255)

                                                                                                                                                                                                                       gttttattttCC‑TTTTTGGGGGGTCCCAAGAAAATTTTTTCCCCATATTTTCGGGAACCAG‑TTTCTGAAAATGCGCTTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATCCCCCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGcc                   <  2:675553/231‑1 (MQ=255)

                                                                                                                                                                                                                         ttA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCggg                >  1:702133/1‑239 (MQ=255)

                                                                                                                                                                                                                           a‑aCTCCC‑ATTTTTGCGGTTATCCAGCAAAATTCTTGCCCATATCTTCGTGAACCAG‑TTCCTGCAATTCGCCCTACTGGCTGCGCAATGCGCGGATCACCACATCGGAACATCCACCGCTTTCTTGCGAAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGAt              <  1:361775/239‑1 (MQ=255)

                                                                                                                                                                                                                               ttCC‑ATTTTTCGGGGTATCCAGCAAAATTCTTCGCCATATCTTGGTGAACCAG‑TTTCTGCAAATCGACCTTCGGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTGGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTg           <  1:772777/239‑1 (MQ=255)

                                                                                                                                                                                                                                    aTTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGTTGGCAAGGTTGGAATTGGCCTGGATTTGTAcc       >  1:145636/1‑239 (MQ=255)

                                                                                                                                                                                                                                         cgCGGGTACTCAGCTAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGCACATAACCGCCAGGATAAGAAAGCTGCCCATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGCCTCGGAGGGCACGGTTGGTATTGGGCGGGGATTGAGCCACGGt  >  2:611421/2‑239 (MQ=255)

GCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGG‑TTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCT‑GTCCGGGCGCTTTA‑ACTTCC‑ATTTTTGCGGGTACTCAGCAAAATTCTTCGCCATATCTTCGTGAACCAG‑TTTCTGCAAATCGACCTTCTGGCTGCGCAATGCGCGGATCACCACATCGGAACATACACCGCGTTCTTGCGGAACATCACCGCCAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGT  >  NC_000913/4137111‑4137578

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