Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,336,594	0	A	C	27.0%	138.0 / 11.6	75	G17G (GGT→GGG)	yfaA	DUF2138 domain‑containing protein YfaA
Reads supporting (aligned to +/- strand):  ref base A (37/17);  new base C (1/19);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.24e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCAT  >  NC_000913/2336347‑2336842                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          ccAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCTGCGG‑AAAG‑CAGTCCTTTTGCGCCAAAACCt                                                                                                                                                                                                                                                        >  2:568099/1‑251 (MQ=255)        aTCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGAc                                                                                                                                                                                                                                                 <  2:243742/251‑1 (MQ=255)           ggCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCAt                                                                                                                                                                                                                                              <  2:68600/251‑1 (MQ=255)                gATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCCACAAGACCATAGAAc                                                                                                                                                                                                                                         >  2:13973/1‑251 (MQ=255)                 aTTCTGATAATAGAAACCGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGAAGCCAGCGG‑AAAG‑CCGTGCNATTGCGCCAAACCCTCCAAGACCATaacacc                                                                                                                                                                                                                                        >  2:360355/1‑246 (MQ=255)                          ataGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCtt                                                                                                                                                                                                                               >  1:326794/1‑251 (MQ=255)                            aaaaaaCGAAATCTTCGCTCAGAACATCGTGGAGAAAGGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTC‑GATGTGGC‑GCGTTTGAGAAAGG‑GGCCCGGCT‑CATTGAGATCGTTTGGCAGATTGTT‑GGTCCCCACCCGGGGCAACAACGTTTTTTTTGGCCCACTGCCAGCATATTGAAACGCCCAGACCCCAGCGG‑AAGG‑CAGTCTTATTGGCCCAAAACCTACAAGCCCATATAACCGCCAGCCTTTc                                                                                                                                                                                                                             <  2:536592/249‑1 (MQ=255)                             aaaaaCGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCg                                                                                                                                                                                                                            <  1:185974/250‑1 (MQ=255)                               aaaCGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGcc                                                                                                                                                                                                                          >  2:69165/1‑251 (MQ=255)                                 aCGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCGGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGGCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCtt                                                                                                                                                                                                                        <  1:13973/251‑1 (MQ=255)                                  cGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAAGGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCCCCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCttt                                                                                                                                                                                                                       <  1:292691/251‑1 (MQ=255)                                   gAAATATTCGCTCAGAACATCGTGGAGAAAAGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCTAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCAGAGAACCGCCAGCCTTTCGCCtttt                                                                                                                                                                                                                      >  1:232184/1‑251 (MQ=255)                                   aaaaTCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACACCGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCAAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCtttt                                                                                                                                                                                                                      <  1:685728/250‑1 (MQ=255)                                       cctgccctaaGAACAGCGTGGAGAAAGGGAATGGTGAGGAGATCTTTGGGAAGGTGGTCCAGA‑GAGTC‑GCTGTTGAGAAAG‑‑GGGCCGGGTGCTTTGATATCGTTGTGCAGATTGTG‑GGGCCCCACCCGCGGCGACAAGGTTTTGTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAACCCTACAGGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTc                                                                                                                                                                                                                  <  2:166084/242‑1 (MQ=255)                                             cTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAAACCATAGACCCGCCAGCCTTTCGCCTTTTTTTCACCACt                                                                                                                                                                                                            >  2:447427/1‑251 (MQ=255)                                                    cagcGTGGAAAAACGGAATGGTGAGGAGATGTTGGGGAAGCTGGCTCAGA‑GAGTC‑GATGGCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCAATTTGCGGATTGTT‑GTGCCCCACCAGCGGCGACAACGTTTTTTCTGGCCCACGGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CGGTGCTATTGCGCCAAACCATCAAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGcc                                                                                                                                                                                                     <  2:4641/248‑1 (MQ=255)                                                    cATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGcc                                                                                                                                                                                                     >  2:233070/1‑251 (MQ=255)                                                             gaAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCtt                                                                                                                                                                                            >  1:518024/1‑251 (MQ=255)                                                             aaaaCGGAATGGTGAGGAGATCTTTGGGAAGCGGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCtt                                                                                                                                                                                            <  1:439953/250‑1 (MQ=255)                                                              aaaCGGAATGGTGAGGAGATCTTTGGGAAGCGGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTGGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtg                                                                                                                                                                                           <  1:187953/251‑1 (MQ=255)                                                                  ggAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTgtgt                                                                                                                                                                                        >  2:640251/1‑250 (MQ=255)                                                                   gAATGGGGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCGTTCGAGAAAGG‑GCCCCGGTT‑CTTTGTGATCTATTTGCAGATTGTG‑GTGCACCACCGGCGGCAACACCGTTTTTTCGGTCCCACTGCCAGCATATTGCAACGCCCAGCCGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAACCCTAAAAGACCATAAAACCGCCAGCCTGTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGtat                                                                                                                                                                                      <  2:551266/251‑1 (MQ=255)                                                                       ggTGGGGAGATCTTGGGGAAGCTGGCTAAGA‑‑AGTCCGGTGTGGA‑TAAGACGGCCCGGCT‑CATTGAGATCGATTTGCAG‑TTGTTTGTGCACCACCAGCGCCGAAAACGTTTTTTCTGGCCCACTCCCACCATATTGCACCCCCCAGACGCCAGCGG‑AAAG‑CAGTTCTTTTCCGCCAAACCCCCCAAGACCAGAGACCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATagcc                                                                                                                                                                                  <  2:455987/251‑1 (MQ=255)                                                                           gggagATCTTGGGGAAGCGGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GTTCCGGCT‑CATTGAGATCTATTTGCAGATTGTT‑GTGCACCCCCAGCGGCGACAACGTTTTTTCTGCCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGT‑AAAG‑CAGTGCTATTTCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATagccagcc                                                                                                                                                                              <  1:230751/250‑1 (MQ=255)                                                                                 tCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt                                                                                                                                                                        >  1:550723/1‑251 (MQ=255)                                                                                         aaGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GAGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑cc                                                                                                                                                               >  2:714565/1‑250 (MQ=255)                                                                                          aGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGTT‑CATTGAGATCGTTTGGCAGATTTTT‑GGGCACCCCCAGCGGCGACAACGTTTTTTCTGCCCCACTGCCAGCATATTGCACCGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCAAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑cac                                                                                                                                                              <  1:599720/251‑1 (MQ=255)                                                                                                  aaga‑gagTC‑GCGGTCGAGAAGGA‑GGTCGGGCT‑CATGGAGATCGATTTGAAGATTGTT‑GTGAACCACCAGCGGCGACAACGTTTTTTGTGGCCCACGGCCATCATATTGAAACGCCTAGACGCAAGCGG‑AAAG‑CAGTGCTATTGCGCCAAACCCTACAAGCCCATAGATCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAAc                                                                                                                                                      <  1:72716/250‑1 (MQ=255)                                                                                                     a‑gagTC‑GCTGTCGAAAAAGA‑GGTCCGGCC‑CATTGAGACGGATTTGCAGATTGTT‑GGCCCCCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACCCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACttt                                                                                                                                                   <  1:168969/251‑1 (MQ=255)                                                                                                                   gagaAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACTACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCCAACTTTACCGTGCCCTa                                                                                                                                        >  2:440461/1‑250 (MQ=255)                                                                                                                    agaAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTTTG‑GTGCCCCACCAGCGGCGACATCGTTTTTTCTGGCCCACTGCCAGCATATGCCAACGCCCAGACGCCAGCGA‑AAAG‑CAGTGCTATTGCGCCAAAACCAACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAAt                                                                                                                                      <  1:35182/251‑1 (MQ=255)                                                                                                                       aaGA‑GGTCCGGTT‑CATTGAGTTCGATTTGCAGATTTTT‑GTGCCCCCCCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGA‑AAAG‑CAGTGCTATTGCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACg                                                                                                                                   <  1:181212/251‑1 (MQ=255)                                                                                                                       aaGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACg                                                                                                                                   >  1:652234/1‑251 (MQ=255)                                                                                                                        aga‑gGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCCTTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACg                                                                                                                                   >  1:379925/1‑250 (MQ=255)                                                                                                                          a‑gGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGAc                                                                                                                                 >  1:161451/1‑250 (MQ=255)                                                                                                                                  gCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCttttttt                                                                                                                                                                                                                   >  2:516232/1‑162 (MQ=255)                                                                                                                                      cATTGAGATCGATTTGCAGATTGTT‑GTGCCCCACCAGCGGCGACATCTTTTTTTCTGGCCCCCTGCCAGAATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCAAAAACCTACAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAg                                                                                                                       <  2:518024/250‑1 (MQ=255)                                                                                                                                      cATTGAGATCGATTTGCAGATTGTT‑GTGCCCCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCa                                                                                                                        <  1:603137/249‑1 (MQ=255)                                                                                                                                       aTTGAGATCGTTTTGCAGTTGGTT‑GGGCCCCCCCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGCCGCCACCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGAc                                                                                                                     <  2:44979/251‑1 (MQ=255)                                                                                                                                         tGAGAGCGATTTGCAGATTGTT‑GTGCCCCACCAGCGGCGACAACGTTTTTTCTGGCCCCCTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACt                                                                                                                    <  1:706954/250‑1 (MQ=255)                                                                                                                                                  tttGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGACAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTg                                                                                                          >  1:600954/1‑251 (MQ=255)                                                                                                                                                          ttgtt‑gTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttc                                                                                                  >  1:607846/1‑251 (MQ=255)                                                                                                                                                            gtt‑gtGCACCCCCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttca                                                                                                 <  1:370933/250‑1 (MQ=255)                                                                                                                                                               tggccaccccCAGGGGCGACAACGTTTTTTCTGGCCCCCGGCCAGCATATTGCAACGCCAAGACGCCAGCGG‑AAAG‑CAGTCCTATTGCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACTACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttcaat                                                                                               <  1:440461/246‑1 (MQ=255)                                                                                                                                                                gtgCACCCCCAGCGGCGACAGCGTTTTTCCTGGCCCACGGCCAGCATATTGCAACGCCCAGACCCCAGCGG‑AAAG‑AAGTGCTATTGCGCCAAACCCTACAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttcaattc                                                                                             <  1:73/251‑1 (MQ=255)                                                                                                                                                                     ccaccaGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCACCCTCTCGCATTTTTTTAACCAATCCTTGCCACCTCCCTTGTGTCTAGCCAGCCATTTTTTACGTT‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGCCTTTGCAGCCTGGACTTTTTAATTCAAAc                                                                                         >  1:410275/1‑250 (MQ=255)                                                                                                                                                                         caGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAggg                                                                                    <  1:69165/251‑1 (MQ=255)                                                                                                                                                                             ggcgACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCAAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAggg                                                                                    >  1:620982/1‑247 (MQ=255)                                                                                                                                                                              gcgACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTTTTGCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATa                                                                               <  1:714565/251‑1 (MQ=255)                                                                                                                                                                               cgccaccGTTTTTTCTGGCCCACTGCCAGCATATTGCACCGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAg                                                                              <  1:447427/245‑1 (MQ=255)                                                                                                                                                                                  aaactttttttCTGGCCCACTCCCAGCAGATTGCAACCCCCAGACGCCAGCGG‑AAAG‑CAGTGTTAATGCGCCAAAACCTAAAAGACCAAAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTc                                                                           <  1:481149/246‑1 (MQ=255)                                                                                                                                                                                       ttttttCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCttt                                                                      >  1:270091/1‑251 (MQ=255)                                                                                                                                                                                          tttgtgccccACTGCCACCATAGTGCAACGCCAAGACGCCCGCGG‑AAAG‑CAGTGCTATTGCGCCAACCCCTACAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                    <  1:169303/243‑1 (MQ=255)                                                                                                                                                                                          ttccggccccACTGCCATAATATTGCAACGNCCAGACGCC‑GCCG‑AAAGGCAGTGTTATTGCCCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTg                                                                    <  2:326794/243‑1 (MQ=255)                                                                                                                                                                                               ggCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATg                                                               >  1:88392/1‑250 (MQ=255)                                                                                                                                                                                                gcccattcccAGCATATTTCCACCCCCAGACGC‑GGCGGGAAAG‑CAATGTTTTTGCGCCAAACCCGTCAACACCATAAAACCGCCCGCCTTTCGCATTTTTTTCACCTCTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CTCAGCCAAACTTTCCCGTGCCCTAATACGACAAAAGCCCAGACTTCGCAGCCTTGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATg                                                               <  1:1691/241‑1 (MQ=255)                                                                                                                                                                                                  ccACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACttt                                                                                                    >  1:450296/1‑210 (MQ=255)                                                                                                                                                                                                       gCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTa                                                       >  1:538536/1‑250 (MQ=255)                                                                                                                                                                                                                 tGCAACGCCCAGACGCCAGCGG‑AAAG‑CATTGCTATTGCGCCAAACCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttct                                             <  1:568099/250‑1 (MQ=255)                                                                                                                                                                                                                    aaCGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAg                                                                              >  1:332609/1‑214 (MQ=255)                                                                                                                                                                                                                     aCGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCCCTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAAAAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttcttc                                        >  2:90958/1‑251 (MQ=255)                                                                                                                                                                                                                       gCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTg                                      >  1:595628/1‑251 (MQ=255)                                                                                                                                                                                                                          cAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATCCGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCt                                   >  2:547078/1‑251 (MQ=255)                                                                                                                                                                                                                           aGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtc                                  >  1:552176/1‑251 (MQ=255)                                                                                                                                                                                                                            gACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtcg                                 >  1:543686/1‑251 (MQ=255)                                                                                                                                                                                                                              cGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGAAAAATTATTCTTCTTCTGGCtcgtc                               >  1:519715/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggg‑AAAG‑CAGTGCTATTGCGCCAAACCCTACAGGACCATAAAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAAc                         <  1:565970/250‑1 (MQ=255)                                                                                                                                                                                                                                      ggataG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATCGAACCGCCAGCCTTTCACCTTTTTTTCACCACTCATTCCCACATTCCTTGTGTCTAGCCAGCCATTTTTTACGGG‑CACAGCCACACTTTACCGTGCCCTAATACGACAAAAGCCCAAACTTTGCAGCCTGGACTTTTCAATTAAAACAAGGGAGATAGCTCCCTTTTGGCAGGAAGAAGTAAAATTCTTCTTCTTATGGCCCGTCTTCAACg                        >  2:148610/5‑251 (MQ=255)                                                                                                                                                                                                                                      g‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGCAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGt                       >  1:601343/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggcaatcAGTGCTGTTGCGCCAAATCCTCCAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGCGTCTAGACAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAACACGACTAACGCACCGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCACCCTTTGGGCATGACGAAGTAAAGTGATTCTTCTTCTGGCTCGTCGTCAACg                        >  2:278005/7‑251 (MQ=255)                                                                                                                                                                                                                                         acg‑CAGTGCTATTGCGCCAAAACCTACAAGACAATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCAATCATTGCAACATTCCTAGTGTATAGCCCGCCATTTTTTACGGGACCCA‑CCAACCTTTACCGTGCCCTCATACGACAAAACCCCAGACTTTGCCGCCTGGACTTTTCAATTCAAACAAGGGAGCTAGCTCCCTTTTGGCATGCAGCAGTCAAATTATTCTTCTTCTGGCACGTCGTCAACGTcc                     >  2:281294/3‑251 (MQ=255)                                                                                                                                                                                                                                                  cTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTT‑CACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCgg             >  1:633369/1‑250 (MQ=255)                                                                                                                                                                                                                                                    atgcctccAAAACCTACAATCCCATAGATCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCg              <  2:491102/242‑1 (MQ=255)                                                                                                                                                                                                                                                    aTTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACAACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAACAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCgg             >  1:65632/1‑249 (MQ=255)                                                                                                                                                                                                                                                      tGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCg         >  1:587140/1‑251 (MQ=255)                                                                                                                                                                                                                                                          ccaaaaccaaaaAGACCAAAGAACCGCCAGCCTTTCGCCGTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGAttt     <  2:161451/251‑1 (MQ=255)                                                                                                                                                                                                                                                             aaaCCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGCTCGATTTCAt  >  2:530687/1‑251 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          CCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGA‑GAGTC‑GCTGTCGAGAAAGA‑GGTCCGGCT‑CATTGAGATCGATTTGCAGATTGTT‑GTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAAG‑CAGTGCTATTGCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGG‑CACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGGAGCGATTTCAT  >  NC_000913/2336347‑2336842

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