Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,227	A→G	G108G (GGT→GGC)	araD ←	L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,227	0	A	G	99.0%	324.5 / NA	100	G108G (GGT→GGC)	araD	L‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (50/48);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAC‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGAC  >  NC_000913/65987‑66466                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  gATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTgcgc                                                                                                                                                                                                                                        >  1:70423/1‑251 (MQ=255)  aTGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGTCGTGAGAATGGCCCAGGACGCCGGGCATTTGCGCTGCAGCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACGTTCCCGGTTTCCCACCCATATTCACCGTTGATTTCTGCGTCGGTCATTTGGCGGGTGCAGGGAAGGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAAGCGACTGTCCCGCCTGCGcc                                                                                                                                                                                                                                       <  2:781459/251‑1 (MQ=255)    gtcGTTATGCACCGCACCTTCGGCATTTTTGCCC‑AGGCAAATGGGCCGGGAGAATGGTCCAGGCCGCCGGGCAGTTGCGCTGCATCGATTCCCTGTTTTTCGAGGGTTTCTACGATGACGTTACCGGTGTCCCACGCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGTTTGCTGGAATCGACTGGCCCTCCTGCGCCCa                                                                                                                                                                                                                                     <  2:359155/249‑1 (MQ=255)      cGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGGGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccccg                                                                                                                                                                                                                                    >  2:575086/1‑248 (MQ=255)       gTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGCATGGACCAGGCCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCCTATTCACCGTTGCTTTCTGCGTCGGTAATTTTGCGGCTGCAGGGAATGGTGCCGGACGAATAGTCGGCGTGGGAGGTGCCGGTTGCTGGCATCGGCAGGCCCGCCCGCGCCCAgct                                                                                                                                                                                                                                  >  2:841327/1‑249 (MQ=255)       gTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccccgct                                                                                                                                                                                                                                  >  2:295607/1‑247 (MQ=255)       gTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGccc                                                                                                                                                                                                                                      >  1:392587/1‑247 (MQ=255)          aTGCACCGCCTCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCAGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAgctagg                                                                                                                                                                                                                               >  2:989699/1‑246 (MQ=255)          aTGCACCGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGACCGCCTGCGCCCAGAtag                                                                                                                                                                                                                                <  1:61680/250‑3 (MQ=255)            gacccgactCTTCGGCATTTTTGCCC‑ATTCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGAtagtg                                                                                                                                                                                                                              <  2:37357/242‑1 (MQ=255)                cGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGCAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTTCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGGCTGGCCCGCCTGCGCCCAGATAgtggcgtg                                                                                                                                                                                                                         >  2:189431/1‑251 (MQ=255)                cGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtg                                                                                                                                                                                                                         >  1:753752/1‑251 (MQ=255)                cGCATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtg                                                                                                                                                                                                                         <  1:684701/251‑1 (MQ=255)                cGCATCTTCGGCATTTTTGCAC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggagtg                                                                                                                                                                                                                         >  2:437666/1‑251 (MQ=255)                  cATCTTCGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAgtggcgtggc                                                                                                                                                                                                                       <  2:689619/251‑1 (MQ=255)                        cGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGt                                                                                                                                                                                                                 >  1:529092/1‑251 (MQ=255)                        cGGCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGt                                                                                                                                                                                                                 <  1:340564/251‑1 (MQ=255)                         ggCATTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTAGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAAATAGTGGCGTGGGGCGAGTg                                                                                                                                                                                                                >  2:838572/1‑251 (MQ=255)                            aTTTTTGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAgggtgg                                                                                                                                                                                                             >  1:553181/1‑246 (MQ=255)                              ttttGCCC‑ATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGTCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTgtgtg                                                                                                                                                                                                            <  2:761169/250‑5 (MQ=255)                                tttgcccATGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGTgggtgtg                                                                                                                                                                                                          >  1:832586/5‑245 (MQ=255)                                       aTGCAAATGGGCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATg                                                                                                                                                                                                   >  1:191033/1‑251 (MQ=255)                                            aatgggCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATgccgcc                                                                                                                                                                                              <  1:478406/248‑1 (MQ=255)                                               gagCCGTGAGAATGGACCAGAACGCCGGGCATTTGCGCTGCATAGATACCCTGTTTTTCGAAGGTTTATACGATGACGTTACCGGTTTCCCACTCATCTTCACCGTTGCTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGGGGGTGGTGCCGGTTGCTGGAATAGACTGGCCCGCCTGCCCCCAGATAGTGGCGGGGCGCGAgggggggtgtgacatgccgccaat                                                                                                                                                                                           >  2:336813/3‑228 (MQ=255)                                                ggCCGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGGGTGTGTGCACAATGCCGCCAATg                                                                                                                                                                                          >  1:306856/1‑251 (MQ=255)                                                   cGTGAGAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGac                                                                                                                                                                                       <  1:934991/251‑2 (MQ=255)                                                        gAATGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGa                                                                                                                                                                                   >  1:813487/1‑250 (MQ=255)                                                         aaTGGACCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGGGGTGCCGGTTGCGGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAAt                                                                                                                                                                                 <  2:729089/251‑1 (MQ=255)                                                           gggcccAGGACGCCGGGCATTTGCGCTGCATCGATCCCCTGTTTTGCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGc                                                                                                                                                                               <  2:197071/247‑1 (MQ=255)                                                              aCCAGGACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTg                                                                                                                                                                            <  1:437666/251‑1 (MQ=255)                                                                caggacgcggggaagttGCGCTGCATCAATACCCTGTTTTTCGAAGGTATCTACGATGACGTTACCTGTTTCACACTCATATGCACCGTTGATTTCTGCGTCGGTCAGTTGGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGTTGGTGCCGGTTGATGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCTTGGCGCGAGTGTGTGTGCACAATGCCTCCAATGGACGGGAACGCCTg                                                                                                                                                                            <  2:536217/234‑1 (MQ=255)                                                                 aggaCGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTAACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCAGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGAt                                                                                                                                                                          >  1:910173/4‑250 (MQ=255)                                                                    gcgcagggCATTTGCGCTGCATCGATACCCTGGTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATGCCCCGTTGATTTCTGCGTCGGGCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATagag                                                                                                                                                                      <  2:235129/246‑1 (MQ=255)                                                                     cGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGGGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTTGGGGGGCACAATGCCGCCAATGGACGGGACTGCCTGATATagc                                                                                                                                                                     >  1:208428/1‑251 (MQ=255)                                                                        cgtgCATTTGCGCTGCATCGATACCCTGTTTTTCGACGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGCGTGTGGGTGCAAAAGGCCGCCCATGGACGGGAATGcatttcagggcagt                                                                                                                                                                  >  2:609605/4‑238 (MQ=255)                                                                            catttcggcttgtcgagtcccttGTTTTTCGAAGGTTTCTACGATGACGTTGCCGTTTTCCCACTCATATTCACCGTTGATTGCGGCGTCGGGCATGTTGCGGGTGCAGGGATTGGTGCCGTAGAAATAGTCGGCGTGGGTGGGTCCGGTGGCGGGAATCGATTGGCCCGCCGGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAAGGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGt                                                                                                                                                              <  2:709655/229‑1 (MQ=255)                                                                            aaTTTGCGCTGCATCGATACCCTGTTTTTCGAAGATTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGt                                                                                                                                                              <  1:157359/250‑1 (MQ=255)                                                                                       cATCGATTCCCTGTTTTTCGAAGGTTTCTACGATGACGTTTCCGGTTTCCCACTCATATTCCCCGTTGATTTCTGCGTCGGTCATTTGGCGGGGGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGTGGCGGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCgta                                                                                                                                                   <  2:537055/251‑2 (MQ=255)                                                                                            aTACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGCGTGTGTGTGCACAATGCCGCCAATGGACGGGGATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCggg                                                                                                                                              >  1:135323/1‑250 (MQ=255)                                                                                            aTACCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCCATGGACGGGAATGCCTGATAGAGCCGTCGGTGAGTTGGCggatcggg                                                                                                                                              >  1:67473/1‑244 (MQ=255)                                                                                              aCCCTGTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCggggg                                                                                                                                            >  1:315787/1‑248 (MQ=255)                                                                                                ccTGTTTTTGGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTAATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGCTGCTGGAATCGACGGGCCCGCCTGCGCCCAGGTAGAGGCGTGGGGCGAGGGGGGGTGCACAAAGACGCCAAAGGACGGGAGTGGATGATAGAGAAGGCGGGGAGTGGGCGGATCggagggg                                                                                                                                          >  2:514757/1‑249 (MQ=255)                                                                                                   gTTTTTCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAggg                                                                                                                                        <  1:709737/250‑1 (MQ=255)                                                                                                      tttCGAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCtt                                                                                                                                     <  1:760431/250‑1 (MQ=255)                                                                                                      ttgcgaagttttCTACGATGACGTTTCCGGTTTCCCACTCATATTCACCGTTTATTTCTGCGTCGGTCATTTGTCGGGTGCGGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTTCCGGTTGCTGGAATCGACTGGCCCGTCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCtt                                                                                                                                     <  1:571396/241‑1 (MQ=255)                                                                                                          gaaGGTTTCTACGATGACTTTACCGGTTTCCCACTCATAGTCACCGTTGATTTCTGCGTCGGTAATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCGGGAATCGACGGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCg                                                                                                                                <  2:338692/250‑1 (MQ=255)                                                                                                             ggTTTCTACGATGACGTTACCGGTTTCCCTCTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCgcac                                                                                                                             <  1:661426/251‑4 (MQ=255)                                                                                                              gTTTCTACGCTGTCGTGTCCGGTTTCCCACTCATATTCACCGTTGTTTGCTGCGGCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAAGAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGATTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCAcc                                                                                                                            <  2:187820/251‑1 (MQ=255)                                                                                                              gTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGAGCACAATGGCGCCAATGGACGGGAATGCCTGATAGAGCCGGCGGAGAGTTGGAGTATCGGGGGAGGGCTTTTTCGCAcc                                                                                                                            >  1:515798/1‑251 (MQ=255)                                                                                                               tttCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCCGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCt                                                                                                                           >  1:45951/1‑251 (MQ=255)                                                                                                                ttCTACGATAACGTTACCGGTTTCCCACTCATATTCACCGTTGTTTTCTGCGTCGGTAATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCATGCGCCCAGATAGCGGCGTGGGGCGAGTGGGTGAGCACAATGACGCCAATGGTGGGGGATGCCTGCTACAGCAGTCGGTGAGATGGCGTATATGAGGAGGGCTTTTTCgcccctt                                                                                                                          >  2:431300/1‑245 (MQ=255)                                                                                                                 tCTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGCCTGGCCCGCCTGCGCCCAGTTAGTGGCGTGGCGCGGGTGTGTGTGACCAATGCCGCCAATGGAAGGGAGTGCCTGGTCGGGCAGTCGGTGGGTTGGCGCATAGGAGGGGGGCGTTTTCgccccttc                                                                                                                         >  2:655068/1‑244 (MQ=255)                                                                                                                  cTACGATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGGTTTTTTGCACCTTc                                                                                                                         >  1:461869/1‑250 (MQ=255)                                                                                                                      gATGACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAAAGCAGTCGGTGAGTTGGCGGATCGGAGGAGGGCTTTTTCGCACCTTC‑AAccc                                                                                                                   >  1:928277/1‑250 (MQ=255)                                                                                                                         gACGTTACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGGGGGCACAATGCCGCCCATGGACGGGAATGCCTGATAGAGCAGTCGGGGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACtt                                                                                                                >  1:517036/1‑251 (MQ=255)                                                                                                                            gttcccGGTTTCCCACTCATTGTCACCGTTAATTTCTTCTTCGGTCATTTTGCGTGTGCAGGGAATGGTGCCTTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGATTCTTCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCAc                                                                                                             <  1:428026/247‑1 (MQ=255)                                                                                                                              tACCGGTTTCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCAcc                                                                                                            <  2:195611/250‑1 (MQ=255)                                                                                                                                ccgttgtCCCACTCATAATCCCCGTTGATTGCTGCGTCGGTCATTTGGGGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGGGGGTGGTGCCGGTTGCGGGAATCGACTGGCCCGCCTGCGCCCAGAGAGTGGCGTGTCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGt                                                                                                         <  2:191033/245‑1 (MQ=255)                                                                                                                                  gttttCCCACTCATATTCCCCTTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGttt                                                                                                       <  1:175925/249‑1 (MQ=255)                                                                                                                                      tCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGa                                                                                                    <  1:427427/250‑1 (MQ=255)                                                                                                                                      tCCCACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGGGGGCTTTTTCGCAACTTC‑AACCACTTCACCGGTTTCGAt                                                                                                   >  1:627667/1‑251 (MQ=255)                                                                                                                                       cccccTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATg                                                                                                  <  1:646951/247‑1 (MQ=255)                                                                                                                                        ccACTCATATTCACCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGGGGGGCGCGAGTGTGTGTGCCCAATGCCGCCCATGGGCGGGAAAGCCTGATAGGGCCGTCGGGGGGGTGGCGTATGGGAGGAGGGCGTTTTTGCCCCTTC‑AAACCCTTCACCggatttggtgc                                                                                                 >  2:534189/1‑242 (MQ=255)                                                                                                                                             cATATTCACCGTGGAGTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTaa                                                                                              <  1:127141/249‑1 (MQ=255)                                                                                                                                                   caccGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGAcc                                                                                      <  1:184322/249‑1 (MQ=255)                                                                                                                                                    accGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGAcc                                                                                      >  1:558049/2‑250 (MQ=255)                                                                                                                                                      cggtGATTTTTGCGGCGGTCATTTGGCGGGGGGAGGGAAGGGTGCCGTAGAAATAGGCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACGGGCCCGCCTGCGCCCAGATAGGGGGGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCata                                                                                   <  2:803028/248‑1 (MQ=255)                                                                                                                                                          gATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGCATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGt                                                                               >  1:491144/1‑251 (MQ=255)                                                                                                                                                           aTTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGGGAGTTGGCGTATTGGAGGAGGGCTTTTTCGAACCTTC‑AACCACTTCACCGGATTCGATGCTAACCACGACCATATCGTc                                                                              >  1:397312/1‑251 (MQ=255)                                                                                                                                                                ggCGTCGGTCATTTTGCGGGTGAAGGGAATGGTGCCGTAGAAATAGTCGGCGGGGGTGGGGCCGGTTGCTGGAGTCGTCTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCgt                                                                         <  2:627816/250‑2 (MQ=255)                                                                                                                                                                  ggTCGGTCAGTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTc                                                                       <  1:968409/250‑1 (MQ=255)                                                                                                                                                                  ggTCGGGCATTTGGCGGTGGCAGGGATTGGGGCCGTAGAAATAGGGGGCGGGGGTGGTGCCGGGGGCTGGAACCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGACTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTGC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTc                                                                       <  2:269596/250‑1 (MQ=255)                                                                                                                                                                   gTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTTCGCCCAGATAGTGGCGTGGCGCGAGTGGGGGTGCACAAAGCCGACAATGGACGGGAATGCCTGGTAGAGCCGTCGGAGAGGTGGCGTAACGGGGGAGGGCTTTTTCGCACCTTT‑AAACCCTTCCCCGGTTTTGGTGCTACccaccaccaacacatccacggtc                                                                       >  2:814474/1‑236 (MQ=255)                                                                                                                                                                        tCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGGAGTCGGGGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTC‑CACCACTTCCCCGGGTTCGATGCTAACCCCGACCATATCGGCCGCGGTccagttg                                                                 >  2:292541/1‑245 (MQ=255)                                                                                                                                                                          aTTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCt                                                              >  1:71857/1‑251 (MQ=255)                                                                                                                                                                               gCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTaa                                                          <  1:170227/250‑1 (MQ=255)                                                                                                                                                                               gCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAAt                                                         <  1:175653/251‑1 (MQ=255)                                                                                                                                                                                      cAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGGCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTC‑AACCACTTCCCCGGTTTCGATGCTAACCACGACCCTATCGTCAGCGGGCCT‑ACTGCTGTAATCGACGcc                                                  >  2:396838/1‑251 (MQ=255)                                                                                                                                                                                       aGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGGGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTC‑AACCAATTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGGCGCCg                                                 >  1:246490/1‑251 (MQ=255)                                                                                                                                                                                        gtgaagggtcccGTAGAAAGAGGCGGCGTGGGTGGTGCCGGTGGTTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCgg                                                <  2:301912/241‑1 (MQ=255)                                                                                                                                                                                         ggAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGGGCGAGTGTGGGTGCACAATGCCGCCAATGGACGGGGATGGCTGATAGAGCAGGCGGTGAGGTGGAGTATCGGAGGAGGGCTTTTTCGCACCTTCCAACC‑CTTCACCGGTTTCGCTGCTACCCACGCCCATATCCTCAGAGGTCATGAT‑GGTTTAATCGACGCCggg                                               >  2:867018/1‑250 (MQ=255)                                                                                                                                                                                         ggAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑CACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAACGGTCATGAT‑GCTGTAATCGACGCCgg                                                >  1:687803/1‑250 (MQ=255)                                                                                                                                                                                             tGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGGGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGt                                           >  1:418529/1‑251 (MQ=255)                                                                                                                                                                                             ggtgaccgtagaatgagctggtgggggtggTGCCGCTTTCTGTCATCGACTGGACTGCCTGCGCCTAGATAGTGGCGTGGCTCCAGTGTGTGGGCACAATTCCGCCAATGGACGGGCATGCCTGATAGAGCAGTCGGTGAGTTTGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑ACCCACTTCACAGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GATGTAATCGACGCCGGAAGGt                                           <  2:210186/227‑1 (MQ=255)                                                                                                                                                                                                  cgtaaaaattAGTCGGGGGGGGGGGTGCCGGTTGCTGGACTCGACTTGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGGGGGTGGGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGAt                                      <  2:707994/242‑1 (MQ=255)                                                                                                                                                                                                  ccGTAGAAATAGTCGGCGTGGGTGGGGCCGGTTGCTGGAAGCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGAt                                      <  1:48962/251‑1 (MQ=255)                                                                                                                                                                                                  ccGTAGAAATAGTCGCCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGACGGTTTGAt                                      >  1:235722/1‑251 (MQ=255)                                                                                                                                                                                                   cGTAGAAATAGTCGGCGGGGGTGGTGCCGGTGGCTGGAACCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGGGTGTGCACAATGCCGCCAATGTACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGTCTTTTGCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATc                                     <  2:48078/251‑1 (MQ=255)                                                                                                                                                                                                     tAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGTAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATcac                                   >  1:982911/1‑251 (MQ=255)                                                                                                                                                                                                            aGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGGCGGTGAGTTGGCGTATCGGGGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGGTGCTAACCACGACCATATCGTCCGCGGTCAAGAT‑GCTGTAATCGACGCCGGGCGGTTTGATCACCAAgccc                            >  2:229186/1‑248 (MQ=255)                                                                                                                                                                                                             gtcgtggggggggtggcgtgTTGCTGGTATCGTCTGGCCCGCCGGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAcgt                           <  2:685297/232‑2 (MQ=255)                                                                                                                                                                                                                ggCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAcgccgc                        >  1:192334/1‑251 (MQ=255)                                                                                                                                                                                                                ggCGTGGGTGGTGCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGAcgccgc                        >  1:422055/1‑251 (MQ=255)                                                                                                                                                                                                                     gtggggggtgggtgtcTGGAGTCGACTGGCCCGCCTGCGCCCAGATATTGGCGTGGCGCGAGTGTGTGGGAACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCGTGTTCGCACCTTC‑AACCTCTTCACCGGTTTCGATGCTAACCACGACCATATCGGCAGCGGTCATGAT‑GCTGTAATCGAAGCCGGAAGGTTTGATCACAAAGACGCCGCGCTcg                   <  2:135638/236‑1 (MQ=255)                                                                                                                                                                                                                       ggggtgCCGGTGGCTGGAATCGTCTGGCCCTCCTGCGCCCAGATAGTGGCGTGGCGCGCGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATTGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTcgcg                 <  1:677789/249‑1 (MQ=255)                                                                                                                                                                                                                         ggtgCCGGTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGTTAGTGGCGTGGCGCGAGTGTGTGTGTACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGAt               <  1:388530/251‑1 (MQ=255)                                                                                                                                                                                                                              cgtgtGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACg          <  1:989699/247‑1 (MQ=255)                                                                                                                                                                                                                              cGGTTGCGGGAATCGACGGGCCCGCCTGCGCCCAGATAGGGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCTGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACg          <  2:45951/251‑1 (MQ=255)                                                                                                                                                                                                                                gTTGCTGGAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGGATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGGTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGACGGAAGGTTTGATCAAAAAAACGCCGCGCTCGCGATCAACGgc        >  1:763802/1‑251 (MQ=255)                                                                                                                                                                                                                                      ggAATCGACTGGCCCGCCTGCGCCCAGATAGTGGCGTGGCGCGAGTGTGTGTGCACAATGCCGCCAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAT‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATTAACGGCGCTGAc  <  1:719547/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  GATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAAACGGGCCGTGGGAATGGACCAGAACGCCGGGCATTTGCGCTGCATCGATACCCTGTTTTTCAAAGGTTTCTACGATGACGTTACCGGTTTCCCACTCATATTCGCCGTTGATTTCTGCGTCGGTCATTTTGCGGGTGCAGGGAATGGTGCCGTAGAAATAGTCGGCGTGGGTGGTGCCGGTTGCTGGAATCGACTGACCCGCCTGCGCCCAGATGGTGGCGTGGCGCGAGTGCGTATGCACAATGCCGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTC‑AACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGAC‑GCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGAC  >  NC_000913/65987‑66466

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