Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,158	0	C	T	20.8%	103.3 / 28.5	81	intergenic (+329/+407)	ycdU/serX	putative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base T (7/9);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

TATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTG  >  NC_000913/1096862‑1097453                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                              tatGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑a                                                                                                                                                                                                                                                                                                         >  2:1090094/1‑301 (MQ=255)  atGTTCTGTACTCTCTTGCGCATTTTGGTGAGTAACTGAGTGACCAGACAGTTGATGTGCATGATTTCCCCACGCCCTACCAGGCGTGAGCGTGGGCACGCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCACTTCGCTCGCCCTGCGCGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑aa                                                                                                                                                                                                                                                                                                        <  2:645896/301‑1 (MQ=17)       cTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATCCCC‑G‑aaaaaaa                                                                                                                                                                                                                                                                                                   >  2:455600/1‑301 (MQ=255)         gTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCCGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAg                                                                                                                                                                                                                                                                                                 >  2:130470/1‑301 (MQ=255)         gTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCGTTAAATATGGCGGTGAGGGGGGGATTGAGTCGCTTCGCTCGCCCGGCGCGCAGCCCGCACCCTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCCGTTGTCGACCCCCGGTCGGGGCCTCTCCTCCCCCCGGTGT‑GTT‑CAATacacccaaacaagg                                                                                                                                                                                                                                                                                                   >  1:332657/1‑291 (MQ=21)           tgctcttccgatcttttGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGc                                                                                                                                                                                                                                                                                                <  2:302074/287‑1 (MQ=255)             cgtctTGCGCATTTGGATTGTGGACGGAGGAACCAGACGGGTGCTGGGCACGATTGCCCCTGGCCCTAACAGACGTGGGCGGGGGCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGGTCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCGTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCg                                                                                                                                                                                                                                                                                             <  1:455119/299‑1 (MQ=17)                    cgcATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACttt                                                                                                                                                                                                                                                                                       >  2:265118/1‑300 (MQ=255)                         tttGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGTCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGcacatcggaagagcgg                                                                                                                                                                                                                                                                                 >  1:301740/1‑286 (MQ=255)                              ttgttgACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCCACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCCTCCCCCCGGTGT‑GTG‑CCATACAC‑G‑AAACAAAAGCCCGTACTTTCGTACGAgta                                                                                                                                                                                                                                                                             >  2:606180/1‑298 (MQ=255)                                   cgatctgTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTcttctt                                                                                                                                                                                                                                                                       <  2:607092/295‑1 (MQ=255)                                     ctggttAACAAGACAGTTGCTGGTGACGATTTCCCCGCGCCCGAACAGACGTGGGCGGGGGCTCCCGGGGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTGCTTTAAAGATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTa                                                                                                                                                                                                                                                                     <  1:264826/296‑1 (MQ=14)                                         gTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGGTCACGGTCTGTCCAACTGGCTGTCGCCAGTTTTCGACCCCCCGTCCGGGCTTCTCAGTCCCCCGGGGT‑GTG‑CACTTCCC‑G‑GACAAAAAGcccagccttccttaccctctctttttagatt                                                                                                                                                                                                                                                                  >  1:606443/1‑272 (MQ=17)                                             ccAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCACTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑Tgg                                                                                                                                                                                                                                                             >  2:358651/1‑300 (MQ=255)                                             ccAGACAGTTGATGTGCACGATTTCCCCGCGCCCTAACGGACGTGGGCGGGGGCACCCGGTGTGTGAAATATAAGAAAAAAAAGCCCGTCCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGc                                                                                                                                                                                                                                                            <  1:1088765/301‑1 (MQ=255)                                              cAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCg                                                                                                                                                                                                                                                           >  2:110157/1‑301 (MQ=255)                                              cAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAG‑CCGCCCCCTGCGTTCCCGGGCCGGTCAACTTGCTTTCCCCCAGTGTCGGGACCCGGTCGGGCCGTCACTCCCCCACCGTGT‑GTT‑CAATTCAT‑G‑AAAAAAACTCCCGTCATTTCCTAcagatctttctcttaatacgtggg                                                                                                                                                                                                                                                           >  2:62793/1‑278 (MQ=12)                                               aGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTTCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCgg                                                                                                                                                                                                                                                          >  2:529477/1‑301 (MQ=255)                                                              aCGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGa                                                                                                                                                                                                                                           >  2:669157/1‑301 (MQ=255)                                                                gggtgcccccgcccacgtctgttagggcgaggggCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATCCAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGAc                                                                                                                                                                                                                                          >  2:1063088/31‑300 (MQ=11)                                                                   ccgcccgcgtcccaaatgggtggggCGGGGGGGCCCGGTGTGGGAAATATACGAAAAAAAAGTCTTGACTTTCGTACGACACCGTGTTTTAAAATGGCGGTCAGGGGGGGATTGATCCGCTGCGCACTCTCTGCGGGCAGCCGGTTCTCGGCGTTGTCGGTCGGCCGAGCTGGCTGTCGCCAGTTGTCGACCCCCGGTGGGGGCTTCTGAGCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACtcgc                                                                                                                                                                                                                                      <  1:130342/279‑1 (MQ=2)                                                                    cccccGCGCCCCAACAGCGGGGGGCGGGGGCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACtcgct                                                                                                                                                                                                                                     <  2:308/300‑1 (MQ=255)                                                                           tgcttagcAGACGAGGCCGCGGTCAACGGGTGAGCGAAATCTGCGAAGAATAGGGCCGTACTTGCGAACGGGGTCTGATCTAAATATGGCGGTGAGGGGGTGATTGCTTTGGTTCGCTCGGCCTGCGGGCAGTCCGCTGACAGCGTTCACGGTCTGTCCAACTGGCTGTGGCCCGCTGTCGACCCTCGCTGGGGGCTGCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTtcgctc                                                                                                                                                                                                                               <  1:301505/293‑1 (MQ=2)                                                                               tAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GCG‑CTATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TTGCGGTGAGGGGGGGATTGACTCGCTtcgctcagat                                                                                                                                                                                                                           >  2:301839/1‑296 (MQ=255)                                                                               tAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGAGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTcgcccgcc                                                                                                                                                                                                                            >  2:775658/1‑299 (MQ=255)                                                                                    tacggtggtgggggcaacatgttggtaaaaTAAAAGAAAGAAAAGCGCGTGCTTTTGTATGAGCGCGTCATTAAATCTGGCGGTGAGTGGGGGAATGACTCGCTTCGCTCGCCCGTCGGGCAGCCCGCTCCTTGCGTTCAAGGTCTGTCGAACTGGCTGTCGCCTGTTCCCTGCCCCCGGTCGGGTCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg                                                                                                                                                                                                                      <  1:201907/272‑1 (MQ=11)                                                                                           ccGGGGTCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGGGCTCTTCTTTAAATATGGCGGTGATGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGTCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc                                                                                                                                                                                                               <  2:746165/299‑1 (MQ=21)                                                                                              ggggCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGGTCGCCCTGCGGGCAGCCCGcc                                                                                                                                                                                                            >  2:394648/1‑299 (MQ=25)                                                                                                         tgtgtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCTCCGGTCGGGGCTTCTCATCCCCCTGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGGACTTTCGTACGATCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTTGCTGCGACCGCCCCGCGGGCACCCCGcctcctccgctca                                                                                                                                                                                                 >  2:1073098/1‑288 (MQ=11)                                                                                                             tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTc                                                                                                                                                                                           <  2:823213/301‑1 (MQ=32)                                                                                                             tacaaattaCGAAAAAAAGGGCCGTAATGTCGCTCAAGCTCTTTTTTGAGTATGGCGTTGGGGGGAGGTTTGAGTAACTTCGTCGGCGCTGGGGGCGGCCCGCTAACTGCATGCACGGTCTGCCCATCTGGCTGTCGCCTGATGCCGGCCCACAGTCGGGGGTGCTCGGCCCCCCGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGt                                                                                                                                                                                            <  1:668444/293‑1 (MQ=2)                                                                                                               aaaagaaCGAAAAAAAAGCCCGTAAGTTCGTACGGGCTCTTCTTTAAATGGGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCGATCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCGGCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTg                                                                                                                                                                                         <  1:1140712/295‑1 (MQ=11)                                                                                                               aaaTATACGAAAAAAAAGCCCGTGCGTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGGCCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTg                                                                                                                                                                                         <  2:968676/300‑1 (MQ=17)                                                                                                                       gAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCaa                                                                                                                                                                                    >  2:704641/1‑298 (MQ=31)                                                                                                                         acaaaagGGCCGTACTTGCTTACGAGTTCTTCGTTAAAGTGGGCGGGGAGGGGGGGGTTGACTCGCTTCGGTCGCCCAGCGGGCAGCCCGTTATCTGCGTTCGCGGCCTGCCCCAATAGCTATCGCCAGTCGCCACCCCCCGCTCGGGGCTTCTCGGCCCCCCGGTGT‑GTG‑CAATAAAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGc                                                                                                                                                                               <  1:394233/299‑1 (MQ=11)                                                                                                                         aaaaaaaGCCCGTAATTTGGATGGAGCTCTGTGTTAAATAGGGCGGGGAGAAGGGTGTTGACTCGCTTCGGTCGGCCTGCGGGCGGCCCGCTCACGGCGTTCACGGTCTGTCCAACGGGCTGCCGCGAGTTGTCGACCCCCGGTGGGGGCTTCTCAGCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGc                                                                                                                                                                               <  1:774770/301‑1 (MQ=11)                                                                                                                                gCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGTGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGCTTGCCCAACTGGACGTCGc                                                                                                                                                                         >  2:706603/1‑300 (MQ=9)                                                                                                                                                        cttTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTTTCCAACTGGCTGGCGCCAGTTGTCGAACCCCGG‑TCgggg                                                                                                                                                 >  2:1122128/1‑299 (MQ=14)                                                                                                                                                          ttAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCCCATCCCCCCGGTGT‑GTG‑CAATACAC‑T‑CAAAAAAAGCCCGTACTTTCGTACCAGCTCTTCTTTAAATA‑TGGCGGTGCGGGGGGGATTGACTGGCTTCGCTCGCCCTGCGGGCATCCCCCTCACTCG‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAgaatgggtgagcacacgtctgact                                                                                                                                               >  2:125862/1‑277 (MQ=11)                                                                                                                                                                                       tcgcttcgctCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTGCGCTCGCCCCGCGGGGAGCCCGCGCACAGC‑GTCTAGTGTCGTTCACACTCGTTGTACCCAGTTTTCTAACCTCGG‑TCCGGGCCTCCCACTCCCTCGgtctgttcctataccgaa                                                                                                                >  1:964569/1‑289 (MQ=11)                                                                                                                                                                                         gcttcgctCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGA‑CCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCCCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGGCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGG‑TCGGGGCTTCTCATCACCCCGGTGTGTGCACTATACGaaaaa                                                                                                             >  2:371173/1‑301 (MQ=2)                                                                                                                                                                                                   cccTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCgg                                                                                                    >  2:461202/1‑300 (MQ=14)                                                                                                                                                                                                     cTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCCCTGACCCCGCCAGCAGCACCCTCACTTC‑GTGAAAGTTCTGTCAAACTGGCTGTCGCCAGTTTTCGAACCCCGG‑TCGGGGCTTACCATCCCCCCGGTAGGTTCAATCTACGAACACAAAGCACGTaa                                                                                                  >  2:680342/1‑300 (MQ=11)                                                                                                                                                                                                      tGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCAAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAACCCGTa                                                                                                   >  2:13324/1‑299 (MQ=14)                                                                                                                                                                                                       gCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACtt                                                                                                >  2:696594/1‑301 (MQ=31)                                                                                                                                                                                                             cgcccgctctcggcgtgcacgctagtccgtaCTGGCGGTCGCCTGTCCCCCACCCCCGCTCGCGGCGTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTAc                                                                                          <  2:16216/271‑1 (MQ=11)                                                                                                                                                                                                              gCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GCTCACGGCCTGGTCGTCGGGCAGGCGGCAGTGTTCGAGCGCCCG‑GCGGGGCTGTTCAGCCCACCGGTGTGTGGAatcgaagcatataaaacccgttatttctctcg                                                                                         >  1:367739/1‑271 (MQ=11)                                                                                                                                                                                                                        tGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTA‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTGCCGGCGCCCTGCCGTCTGCCGGCTCAGTGCAGGTC‑CGGTTTGTCCACCTGTCTGTCTCACAATGGCGAACGACGG‑TCCGTGATGCTCATCCCCCCGATGCGTGAAACTTATGACCAAACAGCTCGCATATTCATACGAGccattatt                                                                               >  1:234298/1‑294 (MQ=2)                                                                                                                                                                                                                            ttaccGGTCTTACCACCTGGCTGTCGCCAGTTCTCGTCCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGATCTCTTCTTTAAATA‑GGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTaaa                                                                           <  2:255641/297‑1 (MQ=9)                                                                                                                                                                                                                                   gccttactaagagtcgtagcaggTTGTCGGCGCCCGGTCGGGGCTTCTCCTCACCCCGGTGT‑GTG‑CAAT‑CACTG‑AAAAAAAGGCTCTTACTTTCGTACGGGCTCTTCGTTAAATA‑TGGCGGTGATGGGGGGGTGGCCTCGCTTCGCTCGCCCGGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGc                                                                     <  1:460721/278‑1 (MQ=11)                                                                                                                                                                                                                                     tGTCCAACTGGCTGTCGCCAGTTGTCGACCCCGGTGCGGGGCTCTTCCTCCCCCCGGTTG‑GTG‑AAATACAC‑GAAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCGGCGGGCAGCCCGCTCACTGC‑GTTCGCGGTCTGTCCACCTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCgg                                                                   <  2:986927/301‑1 (MQ=11)                                                                                                                                                                                                                                     tGTCCAACTGGCTGTCGCCAGTGGTCGACCGCCAGGCGGGGCTTCTCATCCCCCCGGTGA‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTCCTTTCGGACGAGCTCTTCTTTAACTA‑TGGCGGTGCGCGGGGGATTGACGCGCCTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGt                                                                  <  1:695830/301‑1 (MQ=11)                                                                                                                                                                                                                                      gTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAACAAAAAGCCCGTACTTTCGTACGAGCTCTTTCTTAACTATTGCGGTg                                                                 >  2:292576/1‑301 (MQ=14)                                                                                                                                                                                                                                         cAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCCGTTGCCGAACCCCGG‑TCGGGGCTTCCCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTCCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑gg                                                             >  2:184310/1‑301 (MQ=14)                                                                                                                                                                                                                                           actggcgacagCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCACCTGGCCGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑gggg                                                           >  2:22124/11‑301 (MQ=9)                                                                                                                                                                                                                                                ctggcgccagtgtacgcctgccGGTCGGTGCTTTTCCACACCCCGGCGT‑GAC‑CAATACAA‑G‑AACAAAAAGCCCCTGCGTCGGTATGACGTCTTCTTTAAATA‑TGGCGGGGAGGGGGGGATTGCCTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGGGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAt                                                      <  1:1120756/281‑1 (MQ=11)                                                                                                                                                                                                                                                    cGCCAGTTGCCGCCCCCCGTTCGGGCTTTCTCATCCCCCCGGGGT‑TTG‑AAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGAc                                                  <  2:965723/301‑1 (MQ=255)                                                                                                                                                                                                                                                            tgggccgcggatcgtgctttccacaccccccgGGTGT‑‑TG‑CAATACAT‑G‑ACAAAAAAGCCCGTGCTTTCGTACGAGAGCTTCTTTAAAAA‑GGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACGTC‑TTTCACGGTCTGTCCAACGGGCTGGCGCCAGTTGTCGAACCCCGG‑TCGGGGCTCCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACtcgcttcg                                          <  2:1016823/274‑1 (MQ=11)                                                                                                                                                                                                                                                              cGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGGGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACtcgcttcg                                          >  2:781154/1‑299 (MQ=25)                                                                                                                                                                                                                                                               gACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACtcgcttcgctt                                       >  2:665592/1‑300 (MQ=37)                                                                                                                                                                                                                                                                  ccccGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTTCAACTGGCCGTCGCCAGCGGTAGGAGCGTGGTTCGGGGC‑GCTGCTCCTCCCGGTGTGTGCAAACTCCGGAAAAAAAACCCGTACTTTTGTACTCGCGTGTCTTTAAATCTCGCAGTGT‑GGGGAGGATTTGCTCGCTtcgcccgc                                     >  1:71511/1‑300 (MQ=11)                                                                                                                                                                                                                                                                  ccccGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGCGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCAGTCGCCAGTTGTCGAACCCCGG‑TCGGCGCTTCTCATCCCCCCCGTGTCGGCAATAGATGGAAACAAAGACCGTACTTTACTCAGATCTCTTCCTTAAACAAGGCGGTCA‑AGTGGGGAtttcatcacctagcttac                                     >  1:630687/1‑284 (MQ=11)                                                                                                                                                                                                                                                                    ccGGTGGGCGCTTCTCATCCCCCGGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCt                                  <  2:639695/301‑1 (MQ=25)                                                                                                                                                                                                                                                                    ccGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTGACTCGCTTCGCTCGcc                                    >  2:1107975/1‑300 (MQ=18)                                                                                                                                                                                                                                                                        tCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAATCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCg                               <  2:1029862/300‑1 (MQ=37)                                                                                                                                                                                                                                                                            tgCTTCTCAACCCCCCGGTG‑‑GTGTCAATATAA‑A‑AAAAAAAAGCCCGTTCTTTCTTACGAGCTCTTATTTAAATG‑TGGTGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACAGTCTGTCCAACTTGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTTTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa                           <  2:46496/299‑1 (MQ=11)                                                                                                                                                                                                                                                                            ggCTTCTCATCCCCCCGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg                          >  2:72750/1‑301 (MQ=25)                                                                                                                                                                                                                                                                            gactgctcATCCCCCAGTGTT‑GTG‑CAAGAGCA‑A‑AAAAAAAACCCCTTCTTTTCGTACG‑GGTCTTTCTTAAATATTGTCGTGGGGGGGGGGTTTGCCCCGCTTCCCTCGCCTGGCGGCCAGCCCGCTCCCTCC‑TTTAACGTTCTGCCCACCGGCCTGCCCCAATTTGTCGACCCCGGG‑TCGGGG‑TTTTTAACCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg                          <  2:342282/295‑1 (MQ=2)                                                                                                                                                                                                                                                                            acggcagcacacccccGGGGT‑TTG‑GAATAAAC‑G‑ATAAAAAAGCCAGTACTGTAGTACGAGCGCTTCTTTAAATA‑TGGGGTGGAGGGGGCGGTGGACTCGCTGCGCTCGCCCGGCGGCCAGCCCGCTCACTGC‑GTGCACGGTATGTCGAACTGGCGGGCGCCAGTCGCCGGACCCCGG‑TGGGGGCTTCTGAGCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg                              <  1:358261/286‑1 (MQ=11)                                                                                                                                                                                                                                                                              cTTCTCATCCCCCCGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCCGc                         >  2:872664/1‑300 (MQ=21)                                                                                                                                                                                                                                                                               cgtttcttcccgcagatgagag‑CAATACAC‑A‑AAAAAAAAGCCCGTACTTTCGTACGGGCTCTTCTTTAAATA‑TGTCGTTGAGGTGGGGATTGACTCGCTTCGCTCGCCCTGCGGTCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGcc                        <  2:237175/280‑1 (MQ=11)                                                                                                                                                                                                                                                                                      ccccccGGGGTAGTG‑CAATAAAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCCCCCTGCGGGCAGCCCGCTCACTGC‑GTTCCCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc                 <  2:368129/301‑1 (MQ=12)                                                                                                                                                                                                                                                                                       cccccGGTGT‑GTG‑CAATATAC‑G‑AAAAAAAAGCCCGACATTTCGTACGAGCACTTCTATAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCAGCGGGCAGCCCGCTCACTCC‑GTTCACGGTCTGTACAAATGGCTGTCGCCAGTTGTGCAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTg               <  2:239301/301‑1 (MQ=14)                                                                                                                                                                                                                                                                                          ccGGTGT‑GTT‑AAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAAGT‑TGGCGGTGAGGGGGGGATTAACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGt            <  2:723568/301‑1 (MQ=2)                                                                                                                                                                                                                                                                                             tcgt‑gtg‑CAATCCCC‑G‑AAAAACAAGCCCGTACTTACGTACGGGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGACGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc          <  2:419778/298‑1 (MQ=12)                                                                                                                                                                                                                                                                                             atct‑gtg‑CAATATAC‑G‑AAAAAAAGGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTc          <  2:677712/297‑1 (MQ=11)                                                                                                                                                                                                                                                                                                t‑gta‑CAATACAC‑G‑AAAAAAAAGCCTCTACTCTCTTACGAGGTATTAGTTCAATC‑TGGAGTTGAAGGGGTGATTGACTCGCTTCGTTCGGCCTGTGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACg       <  2:520200/300‑1 (MQ=11)                                                                                                                                                                                                                                                                                                 tttg‑CAATATAG‑G‑TAAAAAAAGCCCGTACTTTCGTATGGGCTCTTCTTTAAATA‑TGGCGGTGGGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACgg      <  2:801754/299‑1 (MQ=2)                                                                                                                                                                                                                                                                                                  gtg‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGTCCTGTCCAACTGGCTGGCGCAAGGTGTGGAACCCCGG‑TCGGGGCTTCTCATCACCCCGGTGTGTGCAATATACGAACAAAAAGGCCGTACTTTCATACGAGCACGTCTCTAAATATGGCGGTGA‑GGTGGGGGGTGGAGCGCTTCTCTGGCCCCGCggggcgtccggccatcgtggtacgat      >  1:676998/1‑277 (MQ=1)                                                                                                                                                                                                                                                                                                  gtg‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACggg     >  2:1003540/1‑300 (MQ=11)                                                                                                                                                                                                                                                                                                   tg‑CAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTc    <  2:71571/301‑1 (MQ=11)                                                                                                                                                                                                                                                                                                      cAATATAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGAGCAGCCCGCCGACTGCGTTCCCGTTCtt  >  2:1003374/1‑300 (MQ=2)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                              TATGTTCTGTACTCTCTTGCGCATTTTGATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGT‑GTG‑CAATACAC‑G‑AAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATA‑TGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGC‑GTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGG‑TCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTG  >  NC_000913/1096862‑1097453

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑