Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,893	A→C	A219A (GCT→GCG)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,893	0	A	C	98.9%	331.3 / NA	94	A219A (GCT→GCG)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (53/40);  minor base G (1/0);  total (54/40)
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

TAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTAC  >  NC_000913/3815648‑3816138                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      gaaaCAGCCCGGCGTTGAAGAAATTGGGGCTTTTTCGCCCGGATTTCAGCGTAAACTTGCCAAACTTTAACACCTGCTTGCTAAGCTCAAATTCAATAAACTGGCGCTGATATAGGTTCATGCCTTCGCTCCTCATCTTACTTTTCTACATACAAAAAAAAGTCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑TGGCCAGCGcc                                                                                                                                                                                                                                                    <  2:599708/249‑1 (MQ=255)   aaCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCcaac                                                                                                                                                                                                                                                 >  1:580372/1‑251 (MQ=255)       gCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAgg                                                                                                                                                                                                                                             <  1:507852/251‑2 (MQ=255)             cGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                                        >  1:136522/1‑250 (MQ=255)             cGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGt                                                                                                                                                                                                                                        >  1:353214/1‑250 (MQ=255)                    gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCACATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGGGAGAAGCTCTTc                                                                                                                                                                                                                                >  2:470245/1‑251 (MQ=255)                    gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTc                                                                                                                                                                                                                                >  1:416062/1‑251 (MQ=255)                     aaaTAGGTGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGTTTTCATGCCTTCGCTCCTCATCTTACTTATCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GTGCCAGCGCCAACAAGGTGAGTAGCTCTTCa                                                                                                                                                                                                                               <  1:418348/251‑1 (MQ=255)                           gggCTTTTGCGCCCTGATTTCAACTTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATAGGGTTTCATGCCTTCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACTCAGCAGTCGCCTTAAAAATCAGTTGGCCAGCGCTGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCGC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTg                                                                                                                                                                                                                         <  2:668967/251‑1 (MQ=255)                               ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACg                                                                                                                                                                                                                     >  1:761394/1‑251 (MQ=255)                                  tGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATTCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGc                                                                                                                                                                                                                   <  1:423312/250‑1 (MQ=255)                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                                                                                                                                                                                            >  1:402768/1‑251 (MQ=255)                                            tttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGACCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTCTGc                                                                                                                                                                                                        >  1:757509/1‑251 (MQ=255)                                                          cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                                                                                                                                                                                          <  1:758762/251‑1 (MQ=255)                                                                     aacaaatgcttgcgAAGAGCAACTTCAAAAAACTGGCGCTGACATGGTTTCATGACATTGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTAATCAGTCGCCTTAAAACTCAGTTTACCAGCGCCGCCTTCTGCGTCGCAACAATGGATTCGAT‑T‑C‑C‑CAC‑TC‑GGGCCAGGGCAAACAAGGTGAGTAGATCTTCATGGGAGTCCGGCTCGCCTTCTGCCGTCCCCTGCACCTCAATGATgcgc                                                                                                                                                                               >  2:709674/7‑251 (MQ=255)                                                                        acCTGCTTGCTAAGCGCACATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCCCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAACGATGCGCCCg                                                                                                                                                                            >  2:623177/1‑251 (MQ=255)                                                                            gctcgctAAGCGCACATTCGGTATGCTGGCGCTGATATGGGTTCGAGTCTACGCTCCTCCTCTTGCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGACTTAAAAACCATTGTGCCAGCGTCGCCTTCGGCGTCGCGACAATGGATTCGAT‑C‑C‑C‑C‑CCTC‑TGGCCAGCGCCATCAGGGTGAGTAGGTCTTCATGGGTGAACGGCGCGCCTTCTGCCGGCCCCTGCCCTTCCATTATGCGCCCGTCtt                                                                                                                                                                        <  2:99806/247‑1 (MQ=255)                                                                                     gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGGGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGGTGCGCCCGGCTTCGGtcatcc                                                                                                                                                               >  2:317667/1‑250 (MQ=255)                                                                                       gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGCGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACt                                                                                                                                                             >  2:151293/1‑251 (MQ=255)                                                                                        cAAATTCAATAAACTGGCGCTGAGATGGTTTCATGCCTTCGCTCCTCATCTTTCTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATGC‑C‑C‑C‑C‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTa                                                                                                                                                            <  2:754788/251‑1 (MQ=255)                                                                                         aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTAc                                                                                                                                                           >  1:742893/1‑251 (MQ=255)                                                                                          aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACg                                                                                                                                                          <  1:577556/251‑1 (MQ=255)                                                                                               aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                                                     <  1:686795/251‑1 (MQ=255)                                                                                               aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGGTCAt                                                                                                                                                     >  1:567086/1‑251 (MQ=255)                                                                                               aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCATAAAAATCATTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCACTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                                                     <  2:602887/251‑1 (MQ=255)                                                                                                  aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTc                                                                                                                                                  >  1:193991/1‑251 (MQ=255)                                                                                                   aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCg                                                                                                                                                 >  1:622438/1‑251 (MQ=255)                                                                                                     cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTCAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                                                                                                                               >  2:238483/1‑251 (MQ=255)                                                                                                         cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTctct                                                                                                                                           >  1:50185/1‑251 (MQ=255)                                                                                                           cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                                         >  2:210802/1‑251 (MQ=255)                                                                                                           cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                                         >  1:528305/1‑251 (MQ=255)                                                                                                            tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                                        >  1:408344/1‑251 (MQ=255)                                                                                                             tatataggttCAGGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCTACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGTTTCGAT‑T‑C‑C‑C‑CCTC‑GGTCCAGCGCCAACAAGGTGAGTAGCTCTTCATGTGTGAACGGCTCTCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCTTCTTAGGGCATCACTTCCTTCATGTCGATCTCTGCg                                                                                                                                        <  2:649173/242‑1 (MQ=255)                                                                                                             gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCg                                                                                                                                        <  1:145440/250‑1 (MQ=255)                                                                                                                atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGct                                                                                                                                     <  1:532581/250‑2 (MQ=255)                                                                                                                atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGct                                                                                                                                     <  1:66617/250‑2 (MQ=255)                                                                                                                atGGTTGCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGctg                                                                                                                                    <  1:276153/251‑3 (MQ=255)                                                                                                                       cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTc                                                                                                                             <  2:50185/251‑1 (MQ=255)                                                                                                                        aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTc                                                                                                                             >  2:36407/1‑250 (MQ=255)                                                                                                                                 ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGCCTTCGGTCATCACTACGTTCATGTCGGACCCTGCGGCTGAGTCTTCCACGTATTcc                                                                                                                   >  2:454324/1‑251 (MQ=255)                                                                                                                                   cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                                                   >  1:652634/1‑249 (MQ=255)                                                                                                                                   cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAg                                                                                                                 <  2:289564/251‑1 (MQ=255)                                                                                                                                    ctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                                                   >  1:425536/1‑248 (MQ=255)                                                                                                                                     tcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                                                   >  1:529217/1‑247 (MQ=255)                                                                                                                                        tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCAAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGc                                                                                                            >  2:242161/1‑251 (MQ=255)                                                                                                                                            acatatCTACTGACAAAAAAAAGGCGACTCATCAGGCGACTTAAAAATCACTTCGCCAGCGCCGCCTTCTGCGTCGGTACAATGGATTCGAT‑TAC‑C‑C‑C‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCATCGGTCATCACTACGTTCATGTCGGTCTCTGAGGCTGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                         <  2:730254/245‑1 (MQ=255)                                                                                                                                            aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGACCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaat                                                                                                        >  1:390482/1‑250 (MQ=255)                                                                                                                                                     cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGGTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaattgcttcacc                                                                                               <  2:742893/251‑11 (MQ=255)                                                                                                                                                        aCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGtt                                                                                            <  1:317667/251‑1 (MQ=255)                                                                                                                                                           aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGtttac                                                                                         <  2:528305/251‑4 (MQ=255)                                                                                                                                                                aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAtac                                                                                    >  2:135374/1‑249 (MQ=255)                                                                                                                                                                 aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAAtacc                                                                                   >  1:187261/1‑248 (MQ=255)                                                                                                                                                                     gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                >  2:83748/1‑250 (MQ=255)                                                                                                                                                                     gACTCATCAGTCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATGCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACa                                                                               <  1:606690/251‑1 (MQ=255)                                                                                                                                                                       cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACCATACCGAc                                                                                >  2:516169/1‑248 (MQ=255)                                                                                                                                                                       cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGAc                                                                                >  1:503134/1‑248 (MQ=255)                                                                                                                                                                          atcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAAc                                                                          >  1:745094/1‑251 (MQ=255)                                                                                                                                                                            caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTTCAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTg                                                                        <  1:671926/251‑1 (MQ=255)                                                                                                                                                                             aGTCGCCTTAAAAATCAGGTTGCCCGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGc                                                                       <  2:523267/251‑1 (MQ=255)                                                                                                                                                                             aGTCGCCTTAAAAATAAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑CAC‑A‑C‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGc                                                                       >  1:640445/1‑251 (MQ=255)                                                                                                                                                                                 gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCt                                                                   <  2:416062/251‑1 (MQ=255)                                                                                                                                                                                 gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCt                                                                   >  1:188276/1‑251 (MQ=255)                                                                                                                                                                                      aagaaTCAGTTTGCCAGCGCCGCCTGCTGCGGCGCTACAATGGATTCGTT‑T‑C‑C‑C‑CCGC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTGCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCAt                                                              <  2:585197/248‑1 (MQ=255)                                                                                                                                                                                       aaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGATACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑TGGCCAGCGCCAACACGGGGAGTAGCTCTTCACGGGAGAACGGCTCGCCTTCTTCCGTCCCCTGCACTTCAATGTTTCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACCTATTCCCGAACGCAAATTTCTTCACCGTTTACAATACCGTCAAAGACTGCGGCCACCATc                                                             >  2:220580/1‑251 (MQ=255)                                                                                                                                                                                             ggTTTGCCAGCGCCGCCTTCTGCTTCGCTACNAGGGATTCGAT‑T‑C‑C‑C‑CCGC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCAAGGGTGAACGGGTCGCATTCGGCCGCCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCATTACGTTCATGTCGGTCGCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTc                                                       <  2:334687/250‑1 (MQ=255)                                                                                                                                                                                             aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑C‑CCCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTc                                                       <  2:622438/251‑1 (MQ=255)                                                                                                                                                                                             aGTTTGCCAGCGCCGCCTTCAGCGTCGCTACAAGGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTGTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAAGGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTGACAAGACCGACAGAAACTGCGGCTACCATCCCTTTc                                                       <  2:243553/251‑1 (MQ=255)                                                                                                                                                                                              gTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTc                                                       >  2:604604/1‑250 (MQ=255)                                                                                                                                                                                                ttaccAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑CAC‑C‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGAGAACGGCTCGCCTTCTGCCGTCCCCTGCCCCTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAACTGCGGCTACCATGCCttttat                                                     >  2:756779/4‑247 (MQ=255)                                                                                                                                                                                                     cgcgcCTCCGTCTGCGCCGCTACAAAGGATTCGATGC‑C‑C‑C‑C‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCTTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACTTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCAGCgggtt                                               <  2:580372/250‑4 (MQ=255)                                                                                                                                                                                                     aGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCgggtt                                               <  2:529217/251‑4 (MQ=255)                                                                                                                                                                                                     aGCGCCGCCTTCTGCGTCGCTACAAGGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATTTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCgggtt                                               <  2:481548/251‑4 (MQ=255)                                                                                                                                                                                                        gccgccTTCTGCGTCTCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                            >  1:615787/1‑251 (MQ=255)                                                                                                                                                                                                         ccgccTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGGTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGggttggtt                                           <  1:42565/251‑1 (MQ=255)                                                                                                                                                                                                          cgccTTCTGCGTCGCTACAATGGATTCGCT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCCGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGCCTTAAACGAATTCCCGATCGCAAATTCCTTCACCGTTTACAATACCGACAGAAACTCCGGCCACCCTCCCTTTCATCGGGTTGGtct                                          >  2:291257/1‑249 (MQ=255)                                                                                                                                                                                                          cgccTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGttt                                          >  1:408197/1‑251 (MQ=255)                                                                                                                                                                                                                 tGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCATCCGTCCCCTGCCATTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACTGCGGCTCCAATCCCTTTCATAGGGTTGGTTTTCAGCtt                                   >  2:43602/1‑251 (MQ=255)                                                                                                                                                                                                                  gcgccGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCt                                    <  1:494258/245‑1 (MQ=255)                                                                                                                                                                                                                   cGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGc                                 <  1:400353/251‑1 (MQ=255)                                                                                                                                                                                                                       gCTACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAACTGCGTCTACCATCACTTTCATCGGGTTGGTTTTCAGCTTGCCGtt                             >  1:481560/1‑251 (MQ=255)                                                                                                                                                                                                                         tACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGtttt                           >  2:675160/1‑251 (MQ=255)                                                                                                                                                                                                                         tACAATGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGtttt                           >  1:148489/1‑251 (MQ=255)                                                                                                                                                                                                                           aaaTGGATTCGTG‑T‑C‑C‑C‑CCAC‑GGGCCAGCGCAAACAAGGTTGGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCAACACTACGTTCATGGCGGTCTCGGCGGCTGAGTCTTCAACGTATTCCAGAGCGCAAATTTCTTCACCGTTTACAATACCGACAGAAACTGCGGCTTCCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTcc                         <  2:719131/250‑1 (MQ=255)                                                                                                                                                                                                                              tGGATTCGAT‑T‑C‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTccacc                      >  1:195433/1‑251 (MQ=255)                                                                                                                                                                                                                                gattcgat‑t‑c‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAg                    <  1:210802/251‑1 (MQ=255)                                                                                                                                                                                                                                    cgat‑t‑c‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTc                <  1:675160/251‑1 (MQ=255)                                                                                                                                                                                                                                    cgat‑t‑c‑C‑C‑CCTC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTc                >  1:686177/1‑251 (MQ=255)                                                                                                                                                                                                                                         tcc‑c‑c‑c‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTAGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGATTTCAGCTTGACGTTTTCCACCAGCTTATGt             >  1:304856/3‑250 (MQ=255)                                                                                                                                                                                                                                         tcc‑c‑c‑c‑TC‑GGGCCAGCGCAAACACGGAGAGTAGGTCTTCATGGGGGAACGTCTCGCCTTCTACCGGCCGCTTCACTTCAATGATGCGCCCGTCTGCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCGGAGTCTTCCACGTATTCCAGATCGCAGATTGCTTCCCGGTTTACATTCCCGACCGAAGATGCGGCTACCATCCCCTTCATCGGGTTGGTGTTCACGTTGCCGTTTTCAACCAGCTTCTgc             >  2:402186/3‑249 (MQ=255)                                                                                                                                                                                                                                           c‑c‑c‑c‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGGCCCCTGCACTTCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgc          <  2:652634/251‑1 (MQ=255)                                                                                                                                                                                                                                             c‑c‑c‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcg         >  1:734745/1‑251 (MQ=255)                                                                                                                                                                                                                                               c‑c‑TC‑GGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTATTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTTCAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGtctcgc        >  2:578509/1‑246 (MQ=255)                                                                                                                                                                                                                                                       ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCCTCGGTAGCGCATCtg   >  2:255344/1‑250 (MQ=255)                                                                                                                                                                                                                                                        gCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtgc  >  1:15889/1‑249 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      TAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑T‑C‑C‑C‑C‑TC‑GGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTAC  >  NC_000913/3815648‑3816138

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