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.9%	324.0 / ‑4.8	91	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 (47/43);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  NC_000913/3560213‑3560699                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     ttACCATCAACGTGTTAGCAACGTTGCGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTACCCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATc                                                                                                                                                                                                                                               >  1:524002/1‑251 (MQ=255)   aCCATCAACGAGTTAGCCACGTTCAGATTCTTGGTGAACATGCGCAACTTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGACGCCGACCTCGCTACACAGCCTCGGGCCATTGGGGCTTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCCCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTCACCGTACTGGCAGGC‑AGGGCCGCACCCCCCTGATGG‑CGCAGGATc                                                                                                                                                                                                                                               >  2:279499/1‑249 (MQ=255)    ccATCAACGTGTTACCTACGTTGAGATTGTTGGTGACATTGCGCAATTTGCTGTGATGGAGCAGTGCGTGCGCTCCCGCTTCCGGCG‑TGGTGCCTATATCAATAAACAGCGTCGACCCTTTGGGGATTT‑GCTCCG‑CCATTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAgg                                                                                                                                                                                                                                            <  1:424944/251‑1 (MQ=255)      aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGCTACACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtt                                                                                                                                                                                                                                          >  1:422003/1‑251 (MQ=255)      aTCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGtt                                                                                                                                                                                                                                          >  1:57498/1‑251 (MQ=255)            gttttAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGGGATTGAGCAGTGCGTGCGTTCCCCCTTCGGGCG‑GGGTGCCGATATCGATAAACGGCGTCGGGCCTTGGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                                    <  1:483357/248‑1 (MQ=255)            gtgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGCGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTc                                                                                                                                                                                                                                    >  1:421297/1‑251 (MQ=255)             tgtTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTcc                                                                                                                                                                                                                                   >  2:28609/1‑251 (MQ=255)                      agTTGAGATTGTTGGCGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTcc                                                                                                                                                                                                                          >  2:114388/2‑250 (MQ=255)                        gggagttgtttgGTGACAATGCGCAAATTGGTGTGATTGAGCAGGGCGGGCGCTCCCGCTCCGGGCG‑GGGTGCCGATATCGAAAAAAAGCGTCGACCCTTGGGGGATTT‑GCCCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                                                                                                                        <  1:412096/242‑1 (MQ=255)                           agaTTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                                                                                                                       >  2:429353/1‑249 (MQ=255)                                gtgggTGACAATGCGAAAGTTGCGGTGATTGAGCAGGGCGTGCGCGACCG‑TTCCGGGGGTGGTGCGGATATCGATAAACAGGGTCGACCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                                                                                                                                <  1:475758/248‑1 (MQ=255)                                      gACAATGCGCAAATTGATGTGATTGAGCCGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGa                                                                                                                                                                                                           >  2:559790/1‑250 (MQ=255)                                                   ttGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGACGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                                                                                                                             >  1:33704/1‑251 (MQ=255)                                                          gATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATACCCAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGCGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                        >  1:84198/1‑249 (MQ=255)                                                           aTTGAGCAGTGCGTGCGCTACCGTTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                                                                                                                        >  1:84690/1‑248 (MQ=255)                                                                  agggcgtgctCTACCCCTGCCGGGG‑TGGTGCCGATACGGATAAAGAGCGTCGAGCCATTGGGGATTT‑GCCCCG‑CCACTTGGCGGGCGATGCGCTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTACCCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGTAGAAATGctc                                                                                                                                                                              <  1:541602/248‑1 (MQ=255)                                                                  aGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                                                                                                                              >  1:467366/1‑251 (MQ=255)                                                                      ggtgcgTTACCGCTCCGGGCG‑GGGTGCCGATATCGATAAACATGGTCGAGC‑ATTGGGGATTTTGCTCCG‑CCACTTTGCGGTCGATGCGCTCTTTTTCTTCGGTCTGGTTGGCCTTGCGATTGTGCCA‑G‑GGGGTGTTAACCGAACGGGAAGGC‑GGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                                                                                                                          <  1:259452/250‑1 (MQ=255)                                                                       gtgcgCTACCGCTTCCGGCG‑TGGCGCCGATATCGATAAATAGCGTCGAGCCATTGGGGACTT‑GCTCCG‑CCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGCGCTAACCGAACTGGACGGC‑AGCGCCGCACCGCCAAGATGG‑CGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCGGGGGGCTGGCGGAGAAATgttcaacc                                                                                                                                                                          >  2:461362/1‑243 (MQ=255)                                                                          cgcTACCGTTCCCGGGG‑TGGGCCCGATATCGAAAAACAGCGTCGACCCATTGGGGTTTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                                                                                                                                      <  1:429353/251‑1 (MQ=255)                                                                            cTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCCTGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctct                                                                                                                                                                    >  2:483977/1‑251 (MQ=255)                                                                                       ggCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                                                                                                                         >  1:275823/1‑251 (MQ=255)                                                                                          gggggggCCGATACCGATAAACAGCTTCGAGCCGTTGGGGATTT‑GTTCCG‑CCACTTTGCGGGGGATCCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGTTCGTGCCA‑C‑GGCGTGTTAACCGACCGGGAAGGC‑AGCGCGGCGCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGTTCGCGGCGAATAGGCTGCGGGCTTACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                                                                                                                       <  1:327408/245‑1 (MQ=255)                                                                                            tggtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                                                                                                                                      >  1:614910/1‑250 (MQ=255)                                                                                              gtgCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCA‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                                                                                                                   >  1:361492/1‑251 (MQ=255)                                                                                                      atCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                                                                                                                                           >  1:392656/1‑251 (MQ=255)                                                                                                         aaaaaaCGGCGTCGGGCCTTTGGGTTTTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                                                                                                                                        <  1:196238/248‑1 (MQ=255)                                                                                                               cAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                                                                                                                                  >  2:198083/1‑251 (MQ=255)                                                                                                                ggtttcgcgCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCTCCCTTTTTCTTCGGTCTGGGTGGCCTTGCAATCGTGCCA‑G‑GGCTTGTTAACCGAGCTGAGAGTC‑AGCGCCGCACGTCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATGGAGGTCGCGGCTAATAGTCTGCGGGCTGACGGATAAATGTTCTGCCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                                 <  1:404894/243‑1 (MQ=255)                                                                                                                 ccgtgcacccatgtgggATTT‑GCCCCGGCCACTTGGGGGGCGTTTCCCTTTTTTTCTCGGGTCTGGGGGGCCTTGCGTTCGTCCGC‑C‑GGTGTGTTACCCGAACGGGAAGGC‑ACCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                                                                                                                                 <  2:277092/237‑1 (MQ=255)                                                                                                                      gacccATTGGGGTTTT‑CCCCCG‑CCACTTTGCGGGCGTTGCGCTCTTTGTCTGCGGTCGGGGTGGCCTTGCGATCGGGCCA‑C‑GGCGTTTTAACCGAACTGGAAGGC‑GGCCCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCTCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGTTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                           <  1:594887/248‑1 (MQ=255)                                                                                                                      gAGCCATTGGGGCTTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGCAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGTGAAATGCTCTACCCGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                           >  2:510036/1‑251 (MQ=255)                                                                                                                      gAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                                                                                                                           >  1:351014/1‑251 (MQ=255)                                                                                                                        gctgttGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                                                                                                                            <  1:443992/244‑1 (MQ=255)                                                                                                                        gCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                                                                                                                            >  2:443992/1‑248 (MQ=255)                                                                                                                        gCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                         >  1:397244/1‑251 (MQ=255)                                                                                                                        gCCATTGGGGATTT‑GCTCCG‑CAACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTCACCGAACTGGACGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCAAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                                                                                                                         >  2:404505/1‑251 (MQ=255)                                                                                                                          cATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCg                                                                                                                       >  1:286891/1‑251 (MQ=255)                                                                                                                               ggagactgaacgcgtccACTGTGGGGGCGTTGCGCTCTTTTTAACGGGGCTGGGTGGCCTCTCGATCGGGCCG‑C‑GGCGTGTTAACCGAACCGGAAGGC‑GGCGCCGAACAGCCATGAGGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAGAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg                                                                                                                    <  2:417948/236‑1 (MQ=255)                                                                                                                                     t‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg                                                                                                            >  1:68779/1‑251 (MQ=255)                                                                                                                                              ccACTTTGCGGGCGATGCGCACTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCCAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGTTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGATGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtttgt                                                                                                     >  2:658581/1‑251 (MQ=255)                                                                                                                                                ccTTTGCGGGCGATGCGCTTTTTTTCTGGGTTCTGGGGGGCCTTGGGCGCGTGCCC‑C‑GGCGTGTTAACCGAACGGGAGGGG‑CGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTgtttgttt                                                                                                   <  2:383726/250‑1 (MQ=255)                                                                                                                                                   gtgggggcgaggcgctagtttgcgtcGGTCTGGGTGGCCTTTCGACCG‑GCCA‑CCGGGGTGTTAACCGAGCGGGAGGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATGCCGTTGTGACGTTGTGTTTGTTTCAt                                                                                                <  2:315431/227‑1 (MQ=255)                                                                                                                                                   gtGCGGGCGATGCTCTTTTTTTCTTGGGGCTGGGTGGCGTGGCGATCGGGCCA‑G‑GGGGTGTTAACCGACCGGGAAGGA‑GGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAt                                                                                                <  2:245777/250‑1 (MQ=255)                                                                                                                                                   gtGCGGGCGATGAGCTCTTTTTCGGGGGTGTGGGTGGCCTGGCGACCG‑GCCAAC‑GGCGTGTTAACCGAACGGGAAGGC‑AGCCCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCa                                                                                                 <  2:68779/249‑1 (MQ=255)                                                                                                                                                    ggCGGGCGATGCGCTTTTTTTCTGCGGTCGGGGTGGCCTTGCGATCGGGCCA‑G‑GGCGTGTTAACGGATCTGGAAGGC‑GGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAtt                                                                                               <  2:528593/250‑1 (MQ=255)                                                                                                                                                       gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                            <  1:626316/251‑1 (MQ=255)                                                                                                                                                       gggCGAGGCGCTCTTTTTCTGCGGTCGGGGTGCCCTCGCGCTCGCTCCA‑C‑GGGGTTTTAACCGCCCGGGAAGCC‑AGCCCCACACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTtat                                                                                            <  2:286891/251‑1 (MQ=255)                                                                                                                                                               cgcTTTTTTTCTCCGGTCTGGGGGGCCGTGCGCTCGTGCCA‑C‑GGCGGGTAAACCGAACGGGAAGGC‑GGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                                    <  2:493744/251‑1 (MQ=255)                                                                                                                                                               cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCt                                                                                    >  1:220922/1‑251 (MQ=255)                                                                                                                                                                   ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCg                                                                                                                                >  2:68900/1‑203 (MQ=255)                                                                                                                                                                    tttttCTGCGGTCTGGGGGGCCTGGCGATCGTGCCA‑C‑GGGGTTTTAACCGAACTGGAGGGC‑GGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGaa                                                                                <  2:6217/250‑1 (MQ=255)                                                                                                                                                                             ggccgggtgggCCTGGCAATCGCGCCC‑C‑GGGGTGTAAACCGAACGGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTGCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgt                                                                      <  1:658581/242‑1 (MQ=255)                                                                                                                                                                               tCTGGATGGCCTAGCGATCGTGCCA‑C‑GGCCTGTTAACCGAACTGGAAGGC‑AGAGCCGCAACGCCCTGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTAATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGGTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAATTCCCTGGAATTAttttcgttt                                                                    >  2:55321/1‑251 (MQ=255)                                                                                                                                                                                gtGGGTGGCCTGGCGATCGGGCCA‑C‑GGCGGGTTACCCGAACTGGAAGGC‑AGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgttt                                                                    <  2:84690/249‑1 (MQ=255)                                                                                                                                                                                cTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTAttttcgtttt                                                                   >  2:62411/1‑251 (MQ=255)                                                                                                                                                                                    gTGGCCTGGCGATCGTGCCA‑C‑GGCGGGTTAACCGAACTGGAAGGC‑AGCGCCGCGCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATGTTCGTTTTcgcg                                                               <  1:559790/251‑1 (MQ=255)                                                                                                                                                                                     tGGCCTCGCGCTCGTGCTA‑C‑GCCGTGTTAAACGACCTGGAAGGC‑AGCGCCGCACTGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                              >  2:319515/1‑251 (MQ=255)                                                                                                                                                                                      ggCCTTGCGATCGTGCCT‑C‑GGCGTGTTCACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGGTCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCACTTCGGTACTGACATAACCCTGCTGTTTAACAAGTTCGATAATACCGCTGTGACGTTGTGTTTGTTTCATTTATACATCCCTGGAATTATTTTCGTTTTCgcgc                                                              >  2:99894/1‑250 (MQ=255)                                                                                                                                                                                      ggCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCgcgc                                                              >  1:336765/1‑250 (MQ=255)                                                                                                                                                                                         cTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTg                                                          >  1:275004/1‑251 (MQ=255)                                                                                                                                                                                            gCGATCGTGCCA‑C‑GGGGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGc                                                       <  1:145256/251‑1 (MQ=255)                                                                                                                                                                                             cgtgcGTGACA‑C‑GCCGTTTTAACGGAAGGGGAAGGC‑AGCGCCGTACCGCCATGAGGG‑CGCAGGATCAGGTCTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTg                                                          <  1:292/243‑1 (MQ=255)                                                                                                                                                                                             cGATCGTGCCA‑C‑GGCGTGTTGACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCATGTTTTGCTCCGCCCGCTCATTGGGGTCGAGGCGACACGTCTGCGGGCTGACGGAGAAGTGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACATTGTGTTTGTTTCATGTATCAATACCTTGCATTAgtttcgttttctgtcatttatat                                                      >  2:292/1‑239 (MQ=255)                                                                                                                                                                                              gATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGAGCCGCACCGCCCTGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                                                                                                                            >  2:626830/1‑164 (MQ=255)                                                                                                                                                                                              gATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGCC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                                                                                                                                            <  1:626830/164‑1 (MQ=255)                                                                                                                                                                                                cctggCCA‑C‑GGGGTGTTAACCGAACTGGAAGGC‑GGCGCCTCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAAt                                                   <  2:392656/247‑1 (MQ=255)                                                                                                                                                                                                   gcccA‑C‑GGGGTGTTAACCGACCGGGAAGGC‑AGCCCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaa                                                <  2:397244/249‑1 (MQ=255)                                                                                                                                                                                                   cgggcccgggggggtaacccGAACTGGAAGGC‑GGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATc                                                  <  2:549702/233‑1 (MQ=255)                                                                                                                                                                                                     ccA‑C‑GGCGTGTTAACCGAACTGGTAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATcaac                                               >  2:557937/1‑250 (MQ=255)                                                                                                                                                                                                              gtgtTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGcc                                        >  2:340562/1‑250 (MQ=255)                                                                                                                                                                                                               ggtaaaCCGAACTGGA‑GGCGAGCCCCGCCCCGCCATGAGGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCAt                                      <  2:614910/247‑1 (MQ=255)                                                                                                                                                                                                                gtTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGGAG‑CAGGATCACGTTTTGCTCCGCCAGCTCCCTGGGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAACTGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAAAAAAAGCCATc                                     >  1:620714/1‑251 (MQ=255)                                                                                                                                                                                                                 ttAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCCGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAACTCCCTGGAATTATTTTCGTTTTCGCGCATTGCGCGACTCTACAAAAGCCAt                                      >  2:587957/1‑249 (MQ=255)                                                                                                                                                                                                                 ttAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCATGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATc                                     <  1:199093/250‑1 (MQ=255)                                                                                                                                                                                                                   aaCCGAACTGGAAGGC‑GGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCATTTCGATAATCCCGTTGTGACGTGGTGTTTGTTTCATTTATAAATCCCTGGAATTATGTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCAATCGCt                                  <  1:461362/251‑1 (MQ=255)                                                                                                                                                                                                                     cccgacctggaggccgccgccgcACCGCCATGATGG‑CGCAGGAGCAGGTTTTGCGCCGCCAGCTCATGGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATGTTCGTTTTCGCGCATTTAGCGAATCACCAAAAGCCATCGCTa                                 <  2:453099/234‑1 (MQ=255)                                                                                                                                                                                                                            ggAAGGC‑AGGGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACg                         <  2:275823/251‑1 (MQ=255)                                                                                                                                                                                                                                aggaaGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAACCCCACGGCt                      <  1:36811/247‑1 (MQ=255)                                                                                                                                                                                                                                    cgcgcCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc                  <  2:253664/250‑1 (MQ=255)                                                                                                                                                                                                                                    aGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc                  >  1:202090/1‑251 (MQ=255)                                                                                                                                                                                                                                    aGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAAcc                  <  1:323361/251‑1 (MQ=255)                                                                                                                                                                                                                                      cgccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCGTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTGCATTTAGAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCg                <  1:425012/251‑1 (MQ=255)                                                                                                                                                                                                                                          gcACCGCCATGCTGG‑CGAAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATCACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGCCGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTAGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCgc            >  2:673872/1‑250 (MQ=255)                                                                                                                                                                                                                                           cACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAt           >  1:603938/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCGGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGGTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGCATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTCAACCCACGGCTAACCCGGCGAtg          >  1:386280/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtg          >  1:122771/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCCGCCATGATGG‑CGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTACATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtg          >  2:485390/1‑251 (MQ=255)                                                                                                                                                                                                                                                ccATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGGGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGAtgtgtg      <  1:62750/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gatgG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAGCCCACGGCTAACCCGGCGATGTGTGCTcc  <  2:467366/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     TTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCG‑TGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTT‑GCTCCG‑CCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCA‑C‑GGCGTGTTAACCGAACTGGAAGGC‑AGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTTATAAATCCCTGGAATTATTTTCGTTTTCGCGCATTGAGCGAATCAACAAAAGCCATCGCTAAACCCACGGCTAACCCGGCGATGTGTGCTCC  >  NC_000913/3560213‑3560699

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