Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,455	+G	intergenic (+589/+167)	rtcR → / ← glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,455	1	.	G	98.8%	294.5 / ‑4.5	81	intergenic (+589/+167)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (33/47);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCTTTTACCATCAACGTGTT‑AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAT  >  NC_000913/3560206‑3560690                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            ccttcttTTACCATAAACGTGTT‑AGCAACGTTGAGAGTGTTGGTGACAATGCGCAAATTGTTGTGAGTGAGCAGTGCGGGCGCTACCGC‑TCCCGGCGTGGGG‑CCGATATC‑GATAAACAGCGTCGAGCCTTTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGAGGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cg                                                                                                                                                                                                                                             <  1:203943/250‑2 (MQ=255)  cttcttTTACCATCAACGTGTT‑AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cgc                                                                                                                                                                                                                                            >  1:429612/1‑249 (MQ=255)  cttcttTTACCACCAAC‑TGTTTAGAAACGTTGAGATTTTGGGTGACAACGTGAATATTGCTGTGATGGAGCAGGGGGTGC‑CCAGCGCCTTCCGCTGGGGGG‑ACGATATC‑GATAACCAGCTGCGATCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTTTTTTCCTTGGGTCTGGGGGGCCTTGCGATCGTTCCACGGCTT‑GTTAACCAAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑c                                                                                                                                                                                                                                              <  2:383029/249‑1 (MQ=255)         tACCATCAACGTGTT‑AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCAACGCCATGATGG‑CGCAGGATc                                                                                                                                                                                                                                      >  1:30814/1‑250 (MQ=255)         aaccaccAACGTGTA‑AGCAACGTGGATATTGTGGGTGACAAGGCGAAAATGGCTGTGATGGAGCAGGGCGTGCGCGACCGC‑TTCGGGCGTGGTG‑CCGATACC‑GATAACCAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCAAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCa                                                                                                                                                                                                                                     <  2:71916/250‑1 (MQ=255)           ccATCAACGTGTG‑AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑AATAAACAGCGCCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAgg                                                                                                                                                                                                                                   <  1:168819/251‑1 (MQ=255)                    tgtT‑AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                          >  1:176734/1‑251 (MQ=255)                             aCGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTCCCGT‑TTCCGGGGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                  <  1:324337/251‑1 (MQ=255)                              cGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGATACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                                                                                                                  >  1:302991/1‑250 (MQ=255)                                  gagaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGAAGCGCCGCCACGCCATGTTGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTgt                                                                                                                                                                                                             >  1:303015/1‑250 (MQ=255)                                   agtttgtggggaacaaggctcAAATTGCTGGGAGTGAGCGGGTCGTGCGCTCCCGT‑TCCCGGGGTGGGG‑CGGA‑ATCCGAAAAACGGCGGCGAGCCATTGGGGT‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                                                                                                                            <  2:136386/231‑1 (MQ=255)                                      ttgttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                                                                                                                         >  1:448873/1‑251 (MQ=255)                                       tgttgGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAggcgg                                                                                                                                                                                                        >  2:234244/1‑248 (MQ=255)                                        gttgGTAACAATGCGCAAATGGCTGTGATTGAGCGGTTCGTTCGCTACCGC‑TCCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                                                                                                                        <  1:445713/250‑1 (MQ=255)                                           ggTGACAATGCGAAAATTGGTGTGATGGAGCAGTGCGTGCGCTCCCCC‑TTCCGGCGTGGTG‑CGAATAGC‑GATAAACAGCGTCGAGCCATGGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                                                                                                                                    <  2:147911/251‑1 (MQ=255)                                            gTGACAATGCGCAATGTGTTGTGATTGAGCATGGCTTGCGCTCCCGC‑TTCCGGTGGGTTC‑CCGATATC‑GATAAACAGCGTCTGGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACAATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                                                                                                                                   <  1:129072/251‑1 (MQ=255)                                                aaaTGCGCAATTTGCTGTGATTGGGCAGTGCTTGCGCTACCGT‑TCCCGGCGTGGTG‑CCTATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CCCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                                                                                                                               <  1:195134/250‑1 (MQ=255)                                                 aaTGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                                                                                                                >  1:3819/1‑249 (MQ=255)                                                          atttgtgtgATGGATCAGGGCGTGCGCTCCCGC‑TTC‑GGGGTGGGGGCCGATAGC‑GATAAACAGCGTCGGGCCATGGGGGA‑TT‑TG‑CCCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                                                                                                                                     <  2:173127/246‑1 (MQ=255)                                                            tgctggaattgtggaGTGCGTGCGCTCCCGC‑TTCCGGTGTGGTG‑CCAATATC‑GATAAACATGTTGGAGCCATTGGTGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACTGCGT‑GTTAACCGACCTGGAAGGCAGCGCCGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                    <  1:362445/236‑1 (MQ=255)                                                               tgtgATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGCTATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCTGCCACTTTGAGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                                                                                                                                >  1:205878/1‑251 (MQ=255)                                                                 ggatggatcagggcgtgccCTCCCGC‑TTCGGGGGTGGTG‑CCGATATC‑AATAAACGGCGTGGGGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTGTTCTTCGGTCTGGGTGGCCTTGCGACTGTGCCACGGGGG‑GTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                                                                                                                              <  1:463262/239‑1 (MQ=255)                                                                  gATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                               >  2:154072/1‑249 (MQ=255)                                                                        gcagtgggtggccTACCGC‑TTCGGGGGTGGGG‑CCGATATC‑GAAAAACAGGGTCGAGCCTTGGGGGT‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                                                                                                                                                                       <  1:128248/239‑1 (MQ=255)                                                                                    cttccGC‑TTCAGGCGTGGTG‑CCGATATC‑GATAAATAGCGTCGAGCAATTGGGGA‑TG‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCACCGCACCGCCACGATGGAG‑CAGGAGCAGGTTTGGCTCCGCCAGCACATTGAGGTCGCGGCGAATAGTCTGCGAGCTGACGGAGCAATGctcgacaaacttt                                                                                                                                                           >  2:7853/4‑241 (MQ=255)                                                                                      cccGC‑TTCGGGCG‑GGTGGCCGATATC‑GATAAACGGCGTCGGGCCATTGGGTA‑TT‑GG‑CTCCGCCCCTTTGCGGACGATGCGCTCTTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                                                                                                                          <  2:396430/249‑1 (MQ=255)                                                                                              cggggggggtg‑CAGATACC‑GATAAACAGCGGCGA‑CCAATGGGGG‑TT‑TGTCTCCCCCACTTTGCGGACGAGGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGT‑GTTAACCGAACTGGAAGGCATCGCCGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCCTCGGTACTg                                                                                                                                                  <  2:272174/244‑1 (MQ=255)                                                                                               cccggcggggtgcagttacggATAAAAAGCGTCGGGCATTTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                                                                                                    <  1:132688/230‑1 (MQ=255)                                                                                               cGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCAATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                                                                                                                  >  2:322316/1‑250 (MQ=255)                                                                                                ggCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTg                                                                                                                                                                                                                              >  1:117622/1‑173 (MQ=255)                                                                                                ggCGTGGTG‑CCGATATC‑GATAAACAGCGGCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACAAGCTCTTCGGTACTGAc                                                                                                                                                >  2:406196/1‑251 (MQ=255)                                                                                                    tggtg‑CCGATTTC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑GG‑CCCCGCCCTTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                                                                                                                         <  1:137968/238‑1 (MQ=255)                                                                                                    tggtg‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCAGCTCTTc                                                                                                                                                         >  2:137968/1‑238 (MQ=255)                                                                                                    tggtg‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                                                                                                             >  1:291649/1‑250 (MQ=255)                                                                                                     ggggtgcgattatcaaaaaccggcttCGACCCATTGGGTA‑TT‑GG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCCGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                                                                                                             <  2:207400/226‑1 (MQ=255)                                                                                                       tg‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGGCTGCGGGCTGACGGCGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                                                                                                                         >  2:239142/1‑251 (MQ=255)                                                                                                              tatC‑GATAAACAGCGTCGAGCCATTGGGGC‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCTCGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGGACTGACATAACCCTGCTGt                                                                                                                                   >  2:516867/1‑251 (MQ=255)                                                                                                                     tAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                                                                                                             >  2:272931/1‑251 (MQ=255)                                                                                                                      aaccgccggcgcgcttttGGGGA‑TT‑TG‑CCCCGCCACTTTGCGGCCGATGCCCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTGACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                                                                                                                            <  1:406196/235‑1 (MQ=255)                                                                                                                      aaacagggtcttgCCATTGGGTA‑TT‑TG‑CTCCTCCCCTTGGCGTGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGGCCCCCGGCTT‑TTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                                                                                                                            <  2:43999/239‑1 (MQ=255)                                                                                                                        aCAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGt                                                                                                                          >  2:153858/1‑251 (MQ=255)                                                                                                                           gcgccGAGCCATTGGGGA‑TT‑TG‑CGCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                                                                                                       <  1:322316/247‑1 (MQ=255)                                                                                                                               gggggccttGGGGAGTT‑GG‑CTCCGCCGCTTGGCGGACGAGGGGCTCTTTTTTTTGGGTCGGGGTGGCCTTGCGACCGGGCCACGGCGT‑GTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata                                                                                                                    <  2:385997/244‑1 (MQ=255)                                                                                                                               ccgctctgttgggA‑TTGTG‑CTCCGCCACTTGGCGGCCGCTGCGCTTTTTTTCGTCGTTCGGGGTGGCGTTGCGATCGTGCCACGGCGT‑GTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGATAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                                                                                                       <  2:129572/238‑1 (MQ=255)                                                                                                                                 aGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                  >  1:250525/1‑250 (MQ=255)                                                                                                                                    ggcctgtgggtattggccccGCCCCTTGGCGGCCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGTGCCCCGGCGT‑GTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                                                                                                                 <  2:447831/233‑1 (MQ=255)                                                                                                                                      cctttggggtttgcccccGCCCCTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGT‑GTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc                                                                                                               <  2:3819/235‑1 (MQ=255)                                                                                                                                       tGGGGA‑TT‑TC‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                                                                                                            <  1:490658/250‑1 (MQ=255)                                                                                                                                         gggA‑TT‑TG‑CTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg                                                                                                         >  1:231322/1‑251 (MQ=255)                                                                                                                                          tggggattggcgcagccatgttGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGGGT‑GTTAACCGAACTGGAAGGCAGCGCGGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                                                                                                            <  1:488123/230‑1 (MQ=255)                                                                                                                                           ggggtttagtctcGCCACTTTGCGGCCGGTGCGCGCTTTTTCTTCGTTCGGGGTCCCCTTCCGATCGTGCC‑CGGGGTTCTAAACCGATCTGGACGGCAGCTCCTCACCGCCATGATGG‑CGCAGGATCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGTCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAGCCAGTTCGATAATACCGTtgt                                                                                                          <  2:451410/239‑1 (MQ=255)                                                                                                                                                    cTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                                                                                                              >  2:367328/1‑190 (MQ=255)                                                                                                                                                    cTCCGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                                                                                                                              <  1:367328/190‑1 (MQ=255)                                                                                                                                                      ccGCCACTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt                                                                                               >  1:111648/1‑251 (MQ=255)                                                                                                                                                            cTTTGCGGACGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgtt                                                                                           <  1:150554/249‑1 (MQ=255)                                                                                                                                                              ttGCGGACGATGCGCTCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGGGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCATGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTc                                                                                         <  2:250525/249‑1 (MQ=255)                                                                                                                                                                      gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtatatat                                                                               >  1:462495/1‑248 (MQ=255)                                                                                                                                                                      gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAAt                                                                               >  1:410649/1‑251 (MQ=255)                                                                                                                                                                              ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                              >  2:163182/1‑228 (MQ=255)                                                                                                                                                                              ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt                                                                                              <  1:163182/228‑1 (MQ=255)                                                                                                                                                                                      tcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttt                                                                 <  1:402077/249‑1 (MQ=255)                                                                                                                                                                                        ggTCTGGGTGGCCTTGCGATCTTGCCACGGTGT‑GTTAACCGAGCTTGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGTTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCTGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGTATTAttttcgt                                                             <  1:284161/251‑1 (MQ=255)                                                                                                                                                                                         gtcggggTGGCCTGGCGATCGTGCCACGGGGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtt                                                            <  2:281560/247‑1 (MQ=255)                                                                                                                                                                                             gggTGGCCTTGCGATCGTGCCGCGGCGT‑TTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttttcg                                                        <  2:168870/251‑1 (MQ=255)                                                                                                                                                                                             gggTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttagttttcg                                                        >  2:315705/1‑251 (MQ=255)                                                                                                                                                                                              ggTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgc                                                       >  2:377376/1‑251 (MQ=255)                                                                                                                                                                                               gTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGCTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTcgcg                                                      >  1:340655/1‑251 (MQ=255)                                                                                                                                                                                                tGGCCTTGCGATCGTGCCACGGCAT‑GTTAACCGAACTGGAAATCAGCGCCGCACCGCCATGATGG‑CGCCGGATCAGGCTTTGCTCCGCAAGCTCATTGAAGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTCGATAATCCCGTTGTGACGCTGTGTTTGTTTCATTTATTAATCCCTGGAATTCTTTTCGTTTTCgcgc                                                     >  2:548906/1‑251 (MQ=255)                                                                                                                                                                                                  gCCTTGCGTTCGTCCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAt                                                   <  1:150825/251‑1 (MQ=255)                                                                                                                                                                                                       gCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAGGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGc                                              <  2:176734/251‑1 (MQ=255)                                                                                                                                                                                                             gTGCCACGGCGT‑GCTAACTGAACTGGAAGGCAGCGCGGTATCGCCATGATGG‑CGCAGGATCACGTGTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCGCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAAGACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTTGCGCATTGAGCGAATca                                        <  1:102989/251‑1 (MQ=255)                                                                                                                                                                                                              tggccacgggggtttAACCAAATTGGAGGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                                          <  2:30355/236‑1 (MQ=255)                                                                                                                                                                                                                ccACGGCGT‑GTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCATTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaaca                                     <  2:471452/251‑1 (MQ=255)                                                                                                                                                                                                                        gggtgaaccgaacgtgaggtcggggccgcACCGCCATGATGG‑CGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTGTCGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGGTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATGATGGGCGTTTTCTCGCGTTGAGCGAATCAACAAAAGCCa                              <  2:52979/227‑1 (MQ=255)                                                                                                                                                                                                                        gggtaaaCCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCa                              <  2:429612/246‑1 (MQ=255)                                                                                                                                                                                                                           ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCg                           >  1:382349/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCAc                 >  2:398303/1‑251 (MQ=255)                                                                                                                                                                                                                                      ggAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCa                  <  1:173155/249‑1 (MQ=255)                                                                                                                                                                                                                                        aaGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                  <  1:271913/215‑1 (MQ=255)                                                                                                                                                                                                                                        aaGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCAtt                                                  >  2:271913/1‑215 (MQ=255)                                                                                                                                                                                                                                          ggCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                             >  2:385413/1‑234 (MQ=255)                                                                                                                                                                                                                                          ggCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                             <  1:385413/234‑1 (MQ=255)                                                                                                                                                                                                                                                ggcgcACCGCCGTGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAAGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATTACCCTGCTTTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAGCAACAAAAGCCATCGCTAAACCCACGGCTAACCCgg      <  2:406531/249‑1 (MQ=255)                                                                                                                                                                                                                                                    cACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAt  <  2:231322/251‑1 (MQ=255)                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            TCTTCTTTTACCATCAACGTGTT‑AGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGC‑TTCCGGCGTGGTG‑CCGATATC‑GATAAACAGCGTCGAGCCATTGGGGA‑TT‑TG‑CTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGT‑GTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAT  >  NC_000913/3560206‑3560690

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