Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,883	1	.	C	95.7%	222.8 / 1.7	71	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (39/28);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCC‑T‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815636‑3816132                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           cgcgCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACCCCTGCTTGCTAATCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑c‑cccctc                                                                                                                                                                                                                                                          <  2:190955/251‑4 (MQ=255)  gcgcCAGGTATTAAACAGCACGAGGTTGAACAACTAGGGGCTTTTGCGCCCAGATTTCAGCGTAACCTCGCCAAACTTTCACCCCTGCTTGATAAGCTAAAATTCACTAAACTGGCGCTGATATGCTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGAAGACTCACAAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTTCGTCGCTACATTGGATTCGATT‑c‑cccctcg                                                                                                                                                                                                                                                         >  2:525249/1‑247 (MQ=255)        ggTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCag                                                                                                                                                                                                                                                  <  2:48257/251‑1 (MQ=255)              aaaCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCca                                                                                                                                                                                                                                             <  1:249039/250‑1 (MQ=255)                       ggCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATg                                                                                                                                                                                                                                     >  1:501719/1‑249 (MQ=255)                        gCGTTGAAGAAATAGGGGCTTTAGCGCCCGGATTTCAGCGTAAACGCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGt                                                                                                                                                                                                                                  <  1:189314/251‑1 (MQ=255)                                           ttttGCGCCCGGATTTCATCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CTGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                                                                                                                                                                                                                 <  2:360990/249‑1 (MQ=255)                                           ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                                                                                                                                                                                                                 >  1:360990/1‑249 (MQ=255)                                                       aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTccg                                                                                                                                                                                                   >  2:3171/1‑249 (MQ=255)                                                         ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                                                                                                                                                 >  1:218842/1‑251 (MQ=255)                                                          tCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                                                                                                                                                                                 <  2:426899/250‑1 (MQ=255)                                                           cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCCTGGGGGAACGGCTCGCCTTCTGCCGt                                                                                                                                                                                               >  2:324218/1‑251 (MQ=255)                                                                  aaCTCGCCAAACTTTAACACCTGCTTGCTAAGCTCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTACCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCCGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGACGTCCCCTGc                                                                                                                                                                                        >  1:45800/1‑251 (MQ=255)                                                                  aaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGc                                                                                                                                                                                        >  1:581531/1‑251 (MQ=255)                                                                   ccTCGCCAAACTTTAACACCTGCTTGCTCAGCGCAAATTCAATAAACTGGCGCTGATATCGTTTCATGCCTTCGCTCCTCATCTTACCTTTCCACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa                                                                                                                                                                                       <  2:162476/250‑1 (MQ=255)                                                                       gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAACAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                                   >  2:479500/1‑251 (MQ=255)                                                                          aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCATCATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat                                                                                                                                                                                <  2:488957/251‑1 (MQ=255)                                                                            aCTTTAACTCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatga                                                                                                                                                                              <  1:49593/251‑1 (MQ=255)                                                                            aCTTTAACACCTGCTTGCTACGAGCAACTTCACTCTCCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACACACAAAAGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCAAACAAGATGAGTAGCTCTTAATGGGGGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg                                                                                                                                                                               >  2:69692/1‑250 (MQ=255)                                                                                        gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCg                                                                                                                                                                      <  2:104401/247‑1 (MQ=255)                                                                                        gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCg                                                                                                                                                                      >  1:104401/1‑247 (MQ=255)                                                                                         cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTc                                                                                                                                                                 >  2:68221/1‑251 (MQ=255)                                                                                           tgctAAGCGCAAATGCAATAAACTGGCGCTGATATGGTGTCATGCCTTCGCTCCTCATCTTACGTGTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCGCTGCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATGCGCCCGCCTTCg                                                                                                                                                                <  2:96653/250‑1 (MQ=255)                                                                                            gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCG‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGt                                                                                                                                                              <  1:363406/251‑1 (MQ=255)                                                                                            gcgaaGCGCAAATTCAATAAACTGGCGCTGATATGGTCTCATGCCTTCGCTCCTCATCTTACGTTTCTACAGACAACAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCTCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCG‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCTCGTCTTctt                                                                                                                                                               <  2:1919/247‑3 (MQ=255)                                                                                               aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                                                                           <  2:275641/251‑1 (MQ=255)                                                                                                 gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                                                                          >  1:18230/1‑250 (MQ=255)                                                                                                       aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGt                                                                                                                                                   >  1:555126/1‑251 (MQ=255)                                                                                                          cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTc                                                                                                                                                 >  1:417200/1‑250 (MQ=255)                                                                                                               aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                                                                           >  1:224723/1‑251 (MQ=255)                                                                                                                   ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtct                                                                                                                                       >  2:402048/1‑251 (MQ=255)                                                                                                                      gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTg                                                                                                                                    <  2:319558/251‑1 (MQ=255)                                                                                                                       cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCC‑TTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGc                                                                                                                                   >  1:216577/1‑251 (MQ=255)                                                                                                                              ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCag                                                                                                                              >  1:13005/1‑249 (MQ=255)                                                                                                                              ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCag                                                                                                                              >  2:461199/1‑249 (MQ=255)                                                                                                                                 ttCATGCCTTCGCTCCTAATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTACGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCAGAGTCt                                                                                                                         >  2:205613/1‑251 (MQ=255)                                                                                                                                   cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCtt                                                                                                                        <  1:209728/250‑1 (MQ=255)                                                                                                                                   cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                                                                       >  1:318129/1‑251 (MQ=255)                                                                                                                                         ttctctcctcATCGTACTTTTATACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATACATAGCTCTTCATGGGTGCACGGCTCGCCTTCTACCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCATCACTCCGTTCCTGTAGGTCTCTGCGGCTGCGTCTTCaaaga                                                                                                                  >  2:329716/5‑247 (MQ=255)                                                                                                                                           cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAtt                                                                                                               <  2:104669/251‑1 (MQ=255)                                                                                                                                           cGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTAt                                                                                                                <  1:347981/250‑1 (MQ=255)                                                                                                                                            gCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTc                                                                                                              <  1:578240/251‑1 (MQ=255)                                                                                                                                              tcctcATCTTACTTTTCTCCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAATATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGCTT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGGGTCTTAAACGTATTcc                                                                                                             >  2:285073/1‑250 (MQ=255)                                                                                                                                                 tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTcc                                                                                                             >  1:557025/1‑247 (MQ=255)                                                                                                                                                        aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAATCATTTTGCCAGCGCCGCCGTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCATAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCGGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                   <  2:196932/250‑1 (MQ=255)                                                                                                                                                           tttCTACAGACCAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGACAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACACGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTACCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCCCGATCGCAAACCGc                                                                                               >  2:512771/1‑251 (MQ=255)                                                                                                                                                            ttCTACAGACAAAAACAAGGCGACTAATCAGTCGCCTTAAAAATCAGTTTGACAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑CACCC‑T‑CGGGCCTGAGCCAACAAGATGAGTAGCTCTTCCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCt                                                                                              >  2:198754/1‑251 (MQ=255)                                                                                                                                                               tACAGACTAAAAAAAGTAGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAACACCGCATTCTGCGCCGCTACAATGGAATCGATTTC‑CCC‑T‑CGGGCAAGACACAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTAGCCTTCTGCCGTCCCCTGACCTTCAATGATGCTCCCGTCTTCGGTCATCACTCCGTTCCTGTCGGTCTCTGCGGTAGAGCCTTCAACGCATTCCAAATCGCAATACGCTTCg                                                                                           >  2:390610/1‑251 (MQ=255)                                                                                                                                                                aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCg                                                                                           <  1:458869/250‑1 (MQ=255)                                                                                                                                                                 cagacaaaacacacacGACTCATCAGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAAGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGCATTCCAGATCGCAAACCGCTTCGcc                                                                                         >  2:382096/16‑251 (MQ=255)                                                                                                                                                                        aaaaaaGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAACCCGCTTCGCCGTTcac                                                                                   >  2:297926/1‑250 (MQ=255)                                                                                                                                                                          aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAt                                                                                >  1:63154/1‑251 (MQ=255)                                                                                                                                                                              ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCg                                                                            >  1:59565/1‑251 (MQ=255)                                                                                                                                                                              ggCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCg                                                                            >  2:465339/1‑251 (MQ=255)                                                                                                                                                                                 gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                                                                  >  2:156125/1‑178 (MQ=255)                                                                                                                                                                                 gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGtt                                                                                                                                                  <  1:156125/178‑1 (MQ=255)                                                                                                                                                                                 gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACa                                                                         <  1:6251/251‑1 (MQ=255)                                                                                                                                                                                    tcatcaGTCGCCTTAATAATGACTTTGCCAGCGACGCCTTTTGCGTCGCTACAATGGATTAGATT‑C‑CCCAT‑CGGGCCCGAGCCAACCAGATGAGTAGCTTTTCATGGGAGAACGGCCCGCCGTCAGCCGACCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATCTCGGTCTCTGCGGCCGCGTCTTCAACGTATTCCCGCGCTCACACCCCTTCGCCATTCACACTTCCggcaat                                                                       >  2:19584/1‑245 (MQ=255)                                                                                                                                                                                               cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCAAGATCGCAAACCGCTTCGCCGTTCACAATTCCGCCAGCAACTTCGGCTAc                                                           >  1:370347/1‑251 (MQ=255)                                                                                                                                                                                                  aaaaaTCAGTTGGCCAGCGCCGCCTTCTGCGTCGCAACAATTGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt                                                        <  1:3673/251‑1 (MQ=255)                                                                                                                                                                                                        cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCttt                                                  >  2:443206/1‑251 (MQ=255)                                                                                                                                                                                                          ggtggCCAGCGCCGTCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCa                                                <  1:197045/247‑1 (MQ=255)                                                                                                                                                                                                               ccAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCCTCgg                                            >  1:498672/1‑250 (MQ=255)                                                                                                                                                                                                                    gccgccTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCACAGCCAACAAGATGAGTAGCTATTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGt                                      >  2:513750/1‑251 (MQ=255)                                                                                                                                                                                                                    gccgccTTCTGCGGCGCTACACTGGATGCGATT‑C‑CCCCT‑CGGGCCCGAGCCAACAAGGTGAGTAGATCTTCATGGGTGAACGGCTCGCCTTCTGCCGTACCCTGCACTTCACTGATGCGCCCGTCTTCGGTCATCCCTCCGTTCATGTCGGACTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCCAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCAGCTACCATCCCTTTCATCGGCTTggg                                      >  2:256984/1‑250 (MQ=255)                                                                                                                                                                                                                       gccTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGTTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtttt                                   >  1:448268/1‑251 (MQ=255)                                                                                                                                                                                                                           tCTGCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCTCTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCa                                 <  1:267420/249‑1 (MQ=255)                                                                                                                                                                                                                              gCGTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCt                              <  2:213647/249‑1 (MQ=255)                                                                                                                                                                                                                                gTCGCTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                          <  1:461199/251‑1 (MQ=255)                                                                                                                                                                                                                                    cTACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTTTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCCCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGttt                      >  2:974/1‑251 (MQ=255)                                                                                                                                                                                                                                     tACAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTTTGCAGTCCCTTGCACTTAAATAATGCGCCCGTCTTCGGCAATCACTACGTCCATGTCGGTCTCTGCGGAAGAGCTTTCAACGATTTCCAGATCGCAAACGGCTTGGCCGTTCACAATTCCGACAGAAATTGCGGCTACCATCCCTTCAATCGAATTGGTTTTAAGCTTCCCGtttt                     >  1:140649/1‑251 (MQ=255)                                                                                                                                                                                                                                       cAATGGATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcc                   >  1:163094/1‑251 (MQ=255)                                                                                                                                                                                                                                           ggATTCGATT‑C‑CCCCT‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccacca               <  1:177496/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ttcc‑ccc‑T‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt       >  1:144700/4‑251 (MQ=255)                                                                                                                                                                                                                                                    tcc‑ccc‑T‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt       >  1:194048/3‑250 (MQ=255)                                                                                                                                                                                                                                                    tac‑ccc‑T‑CGGGCCAGAGCCAAAAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGACGTTTTCCCCCAGCTTCTGt       >  2:209260/3‑250 (MQ=255)                                                                                                                                                                                                                                                      c‑ccc‑T‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgc    >  1:385941/1‑251 (MQ=255)                                                                                                                                                                                                                                                        ccc‑T‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATGCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcg   <  1:587809/251‑1 (MQ=255)                                                                                                                                                                                                                                                         cc‑T‑CGGGCCAGAGACAACAAGCTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACCATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCTACCAGCTTCTGTAgcgc  >  1:412555/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                           CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATT‑C‑CCC‑T‑CGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815636‑3816132

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