Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,847	C→T	G497G (GGG→GGA)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,847	0	C	T	98.7%	203.6 / ‑4.9	76	G497G (GGG→GGA)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (35/40);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGC  >  NC_000913/66600‑67090                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          cgagcttgtggtccggaaaaGACGCGCTATCGGCGCAACATGCGTAGAAGGCCGGGTTCGCTTTGCTTATCCGGCCTTCAAAAACGGACCGCGTAGGCC‑TGATAAGACGCGTCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCATTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCg                                                                                                                                                                                                                                                        <  2:6456/233‑1 (MQ=16)      gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATa                                                                                                                                                                                                                                                   <  1:360681/251‑1 (MQ=18)        gttggCCTGATAAGACGCGCTAGCGGCGCATCAGGCGTTGTATGCCGGTTGCGCTTTGCTTATACGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGGGTCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCTGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                                                                                                                                                                                                  <  2:408807/247‑1 (MQ=18)            gCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTc                                                                                                                                                                                                                                             <  1:145343/251‑1 (MQ=18)                 aTAAGACGCGCCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACTCCTGATGCGACTACTGAGCGACGAAATCCGTAATACACTGCGTTcc                                                                                                                                                                                                                                        <  2:640660/251‑1 (MQ=18)                        gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTTTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcag                                                                                                                                                                                                                                 <  2:429084/251‑1 (MQ=18)                                    cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGGCTTTaa                                                                                                                                                                                                                      >  1:555911/1‑250 (MQ=18)                                     aTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAc                                                                                                                                                                                                                    <  1:224452/251‑1 (MQ=18)                                       cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAACTCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGc                                                                                                                                                                                                                  >  1:630660/1‑251 (MQ=18)                                          gCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTgg                                                                                                                                                                                                               >  1:355268/1‑251 (MQ=18)                                           cGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGc                                                                                                                                                                                                              >  1:328889/1‑251 (MQ=18)                                                 attccGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                                                                                        <  2:234386/248‑2 (MQ=18)                                                 aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgg                                                                                                                                                                                                        <  1:467564/251‑2 (MQ=18)                                                   tccgtattcgcTTTGCTTATCCGGCCTACAAAATCTGAACGCGTAGTCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTTAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGTCTACTTAGCGACGAAATCCGTAATACTCTTCTTTCCAGCGCAGCGCTTCTTTAAACGCTGGCAGGCggg                                                                                                                                                                                                       <  1:609674/242‑3 (MQ=18)                                                     cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATAAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggt                                                                                                                                                                                                     >  2:512529/1‑246 (MQ=18)                                                      ggATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgggtatc                                                                                                                                                                                                   <  1:305698/251‑7 (MQ=18)                                                        aTGCGCGTTGCTTATCCGGCCTAAAAAATCGGAACGCGTATGCC‑GGTTAAGACGCGCATGCGTCGCTTCTGGCGTTTAATGCCGGATTCTCTTTGCTTTTCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCTTCATCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGTCAGGCGGGTA‑TCGt                                                                                                                                                                                                <  2:447001/251‑1 (MQ=19)                                                          gcgcTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCG‑TCGGCCTGATAAGACGCGTCAGAGTCACATCAGGCGTTACACACCAGATGCGACTACTTAGCGACGAAATCCGAAATACACTACAGTCCAGCGCAGCGCGTCCTTAAACGCaggcaggcaggttgcgtt                                                                                                                                                                                                >  2:116528/1‑240 (MQ=17)                                                             gttgtcttagcggggcTACAAAATCGGAACGCGTAGGCCGTG‑TAAGACGCGCCCGCGTCGCATCGGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAACGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCGTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTGAAACGCTGGCAGGCGGGTA‑TCGTTatca                                                                                                                                                                                           <  2:247892/236‑1 (MQ=18)                                                                    tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGTG‑CAGGCCTGGTAAGACGCGTCCGCGTCGCATAAGGCGCTACACACCTGACGCGACGACTTCGCGTCGAAATCCGGAATACGCTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggactcgctacaacaccac                                                                                                                                                                                     >  2:169781/1‑231 (MQ=12)                                                                    tATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTACTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CAc                                                                                                                                                                                     <  1:234613/249‑1 (MQ=21)                                                                     aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgggggggatagttatccatccccgg                                                                                                                                                                                    >  2:24034/1‑228 (MQ=18)                                                                     aTCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGt                                                                                                                                                                                  <  1:24034/251‑1 (MQ=21)                                                                        cGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTTTCAAT‑CACCGTggt                                                                                                                                                                               >  1:295393/1‑249 (MQ=21)                                                                                       gaacgcgTAGGCC‑TGATAAGACGCGCCAGCGTCGCGTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAACATCGCAGCGTG‑TAGGCCTGATAAGACACGTCAGCGGCGCATCAGACGTTACACACCAGATGCGCCTACTTAGCAACGACATCCGTTATACACTTCGTTCCAGCACAGCGCGTCTTTAACAGCTGGCCGGCGGGGA‑TCGGTATCAAT‑CAccacgagtgtcatgttgtgagt                                                                                                                                                                >  2:432537/4‑231 (MQ=12)                                                                                         acgcgTAGACC‑TGATAAGACGCGCCAGCGTCGCATGTGGCGTTGGATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGATAT‑CCTGCTAAGATGCGTCAGCGACGAATCATGCGTTACACACCTGATGCGAAAAATTATGGACGAACTCCGTCATACACTTCAGTCCAGCGCAGCGCGTCTTTAACCGCTTGCAGGAGGGTG‑GTGTTATACAT‑CACCGTggtgtctatgcagtgacatc                                                                                                                                                              >  2:73815/2‑232 (MQ=9)                                                                                             gTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGc                                                                                                                                                         >  1:280445/1‑251 (MQ=25)                                                                                                      gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGgc                                                                                                                                                 <  1:485226/251‑1 (MQ=25)                                                                                                              gcgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCg                                                                                                                                         <  1:203741/251‑1 (MQ=25)                                                                                                                        cGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CAACCTGATTTTCATGTCGTA‑CATCTCGGCGacccggtgcatacacttgtggttc                                                                                                                                >  2:417907/1‑227 (MQ=21)                                                                                                                         gaaTCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGGCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTAGCAAT‑CACCGGGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAg                                                                                                                              <  2:295393/249‑1 (MQ=25)                                                                                                                          cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgc                                                                                                                             <  1:512529/251‑2 (MQ=37)                                                                                                                              aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCCTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTTGGGTTCAgcgcct                                                                                                                         >  1:173148/1‑246 (MQ=37)                                                                                                                               ggCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGGGATTTCAATGTCGGG‑CATCTCGGCGAACAGGCGCATATCGTTGAGGGTCAGCGCATg                                                                                                                        >  1:637094/1‑251 (MQ=25)                                                                                                                               ggCGTTGAATGCCGGATGCGCTTTGCTTAGCCGGCCTACAAAATCGCAGCTTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATg                                                                                                                        <  1:699048/251‑1 (MQ=255)                                                                                                                                   ttGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCG‑GAGATCGTTATCAAT‑CACCGTGATTTCAATGTCGGG‑AATCTTAGCGAACCGGCGCataacggttaggatcagcgccaggctg                                                                                                                    >  1:380731/1‑227 (MQ=21)                                                                                                                                           cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGGTCCGCGCATCGC‑T‑GAAGACggg                                                                                                          >  1:428156/1‑250 (MQ=37)                                                                                                                                            ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGGACCGGCGCATATCGTTGGGGGTCCGCGCATGGC‑T‑GAAGACggg                                                                                                          >  1:71588/1‑249 (MQ=25)                                                                                                                                                  gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTgc                                                                                                   <  1:31632/251‑1 (MQ=38)                                                                                                                                                        cTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATCAAT‑CACCGTGATTTCCATGTCGTG‑CATTTCGGCGAAATGGCGCATATCGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGccc                                                                                             >  1:234604/1‑250 (MQ=37)                                                                                                                                                         ttATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGaactg                                                                                                                                                   <  1:80130/198‑4 (MQ=37)                                                                                                                                                           aTCCGGCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGGTTTTAAT‑CAACGTGTTTTCAATTTC‑TGTCATCTGGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCCTGGC‑TTGCAG‑CGGTCTGGTGCGCGccaccc                                                                                          >  2:298902/1‑250 (MQ=19)                                                                                                                                                              cGGCCTACAAAATCGCAACGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGA‑TCGTTATTAATACACGGT‑ATTTCAATGTCGTG‑CATCTTGGCGTACTCGCGCATATCGGTGGGGGTCAGCGCATGGC‑T‑GGAGACGGGATTGGGCGGGCCACCcgcg                                                                                       >  2:347379/1‑251 (MQ=255)                                                                                                                                                                gCCTACAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACAAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAG‑TCGGTATCAAA‑TAACGAGATTTCATTATCGTT‑CATTTCGGCGAACTGGCCCCTATAGGTGAGGGTTACGGCATGGC‑T‑GAAGACgggagaggtggtgccacccgca                                                                                       >  1:370894/1‑229 (MQ=255)                                                                                                                                                                    aCAAAATCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTAGTTCCAGCGCAGCGCGTCTTTAAACGCTGGAAGGCGGGTA‑TCGTTATCCCA‑CACCGTGTTTTCCATATAGTG‑CATATCGGCGAATTGGCGCATATCGTTTGGGGTCAGCGCCTGTCAT‑GA‑GACGGTATGGGGGGGGCCACCCGCGAGGATc                                                                                 >  2:534223/1‑251 (MQ=255)                                                                                                                                                                       aaagcGCAGCGTG‑TAGGCCAGATAAGACGCGGCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATTGTGCGCGCCACCAGCGAGGATCCAc                                                                              <  2:301142/247‑1 (MQ=255)                                                                                                                                                                       aaaacGCAGCGTG‑TAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCAc                                                                              <  2:328889/247‑1 (MQ=255)                                                                                                                                                                       aaaTCGCAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGGATGGTGCGCGCCACCAGCGAGGATCCAc                                                                              >  1:484964/1‑251 (MQ=255)                                                                                                                                                                             cAGCGTG‑TAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg                                                                        <  2:355268/251‑1 (MQ=255)                                                                                                                                                                              aGCGTG‑TAGGCCTGACAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCTGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGAACCACGCTTCgg                                                                       >  1:432032/1‑251 (MQ=255)                                                                                                                                                                                 gtg‑tAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GGAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGGTTCGGa                                                                      >  1:478823/1‑249 (MQ=255)                                                                                                                                                                                 gtg‑tAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGAATGGTGCGCGCCACCAGCCAGGATCCCCGCTTCGGa                                                                      >  1:66449/1‑249 (MQ=255)                                                                                                                                                                                      aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGCTTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCCTATCGTTGTGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGGGCGCCCCCAGCGGGGAGCCCCGCTTCGGAAGCAgg                                                                >  2:381247/1‑250 (MQ=255)                                                                                                                                                                                       ggccttaTAAGACGCGGCAGCGTCGCATCAGTCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGGCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTGCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGtt                                                               <  2:238413/245‑1 (MQ=255)                                                                                                                                                                                       ggCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcg                                                                                                                            <  1:586858/192‑3 (MQ=255)                                                                                                                                                                                           tgATAAAACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAA‑TCGTTATCAAT‑CACCGTGGTTTCAATGCAGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCCAGGA‑T‑TAAGACGGTAAGTTTCGCGCCCCCATCGAGGGTCCACGCTTCGGGAGCAGTTGGc                                                            >  1:127941/2‑250 (MQ=255)                                                                                                                                                                                              tAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGGACTGGCGCATAACGTTGAGGGTCAGCGCATGGC‑T‑GAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGGAGCAGTTGGCAGa                                                         >  1:296606/1‑250 (MQ=255)                                                                                                                                                                                               aaGACGCGTCAGCGTTGTATCAGGCGTTACACCCCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTGCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTGCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGAt                                                        <  2:663465/250‑1 (MQ=255)                                                                                                                                                                                               aaGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGAt                                                        <  1:381247/250‑1 (MQ=255)                                                                                                                                                                                                aGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATcc                                                      <  1:671835/251‑1 (MQ=255)                                                                                                                                                                                                aGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACGTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATcc                                                      <  2:194814/251‑1 (MQ=255)                                                                                                                                                                                                  aCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCgg                                                    <  1:264377/251‑1 (MQ=255)                                                                                                                                                                                                        cAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTgcg                                              <  1:290081/250‑1 (MQ=255)                                                                                                                                                                                                           cGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGGG‑CATCTCGGCGAACTGGCGCCTATCGTTGAGGTTCAACGCATGGC‑T‑GAAGACGGGATGGGGCGCGCCACCAGGGCAGATCCACGCCTTGGAAGcgggttggagatacgggttcgcttt                                           >  2:454739/1‑227 (MQ=255)                                                                                                                                                                                                            gTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTc                                         <  1:181078/251‑1 (MQ=255)                                                                                                                                                                                                                              cacctgATGCGACTACTTAGCGACGAAATCCGTAATACACGTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTT‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCatc                       <  2:285378/246‑3 (MQ=255)                                                                                                                                                                                                                               acctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGCG‑CATCTCGGCGGACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GGAGGCCGTATTGGGCCCGCCACCAGAGAGGATCCACGCCTCCGGAGCCGTTGGGAGGTCCGGGTTG‑GCTTTccccaccgccatcgcccccc                      >  2:625289/5‑233 (MQ=255)                                                                                                                                                                                                                                cctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCgg                     >  1:318767/4‑251 (MQ=255)                                                                                                                                                                                                                                cctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCgg                     <  1:347379/248‑1 (MQ=255)                                                                                                                                                                                                                                   gATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCa                   <  2:637094/250‑1 (MQ=255)                                                                                                                                                                                                                                       cgactactgatctacgaAATCCGTAATACACATCGTTCCAGCGCTGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCGCTGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc              <  2:622112/238‑1 (MQ=255)                                                                                                                                                                                                                                       cgacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc              <  1:514594/248‑1 (MQ=255)                                                                                                                                                                                                                                        gacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCCTATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGACACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCCCATTCGCCACCGGCAGTTTCg             >  1:129301/3‑251 (MQ=255)                                                                                                                                                                                                                                        gacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCATTTGGCAGATCCGG‑TTgc                                               <  1:559675/215‑1 (MQ=255)                                                                                                                                                                                                                                         acTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAA‑CACCGTGATTTCAATGTCGTG‑CATCCCGGCGAACTTGCGCCTATCGTTGAGGTTCAGCGCATGGC‑T‑GCAGACGGTATGGGGCGCGCCACCCGCGAGGCACGCCGCTCCGGAAGCAATTGGCCGGTCCGG‑TTGCGATTTCCACCACGCCTTCGCCACCGGCAggctctg            >  2:22628/2‑245 (MQ=255)                                                                                                                                                                                                                                                   gacgaAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTA‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGc  <  1:702729/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          GCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCC‑TGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTG‑TAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTG‑TCGTTATCAAT‑CACCGTGATTTCAATGTCGTG‑CATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGC‑T‑GAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGG‑TTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGC  >  NC_000913/66600‑67090

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: ‑