Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,686,538	T→C	18.3%	S879G (AGC→GGC)	bcsC ←	cellulose synthase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,686,538	0	T	C	18.3%	187.4 / 7.5	99	S879G (AGC→GGC)	bcsC	cellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base T (39/37);  new base C (6/15);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 8.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTA  >  NC_000913/3686305‑3686772                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              aCCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGGTATCCTCGCTGTAACTGATGCCGCCGACCCCATCCACCACGTTGAAGCCCATCGGCCCGGGACCGATATCCCAGCTCCAg                                                                                                                                                                                                                                           >  2:420568/1‑238 (MQ=255) aCCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCCCTGTAACTGATGCCGCCGCCCCAATCCCCCACGTTGAGGCCCATCGGCGTGGTACCGACATCCCAGcccccca                                                                                                                                                                                                                                          >  1:495229/1‑236 (MQ=255)   cTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGCCCCGATATCATCTCTGTAACCGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCCCCAGAc                                                                                                                                                                                                                                         >  1:656058/1‑238 (MQ=255)            aGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGGGGTACCGATATCCCCGCCCCAGACGGCATTTcc                                                                                                                                                                                                                                >  1:804835/1‑237 (MQ=255)                cAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCCCCGACCACATCCACCCCGTTGAAGCCCATCGGCCTGGTCCCTCAATCCCAACCCCAACCGCCATTTCGcaac                                                                                                                                                                                                                            >  1:494028/1‑235 (MQ=255)                       aGCCCCACCCCGGCGGCCCGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGc                                                                                                                                                                                                                    <  2:817307/239‑1 (MQ=255)                        gCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGGGGTACCGGTATCCCAGCTCCAAACGGCCATTCCCCAGCCGACCGc                                                                                                                                                                                                                    >  2:308153/1‑238 (MQ=255)                             acacCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGAAATCCCAGCTCCAGACCCCATTTCTCCCGCCGACCGCCAcacc                                                                                                                                                                                                              >  1:514857/1‑236 (MQ=255)                             acacCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGGACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCAc                                                                                                                                                                                                                 >  2:137056/1‑236 (MQ=255)                             acacCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGAGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGAAATCCCAGCTCCAGACGCCATTTCTCCAGCCGACCGCCAc                                                                                                                                                                                                                 >  2:844871/1‑236 (MQ=255)                                  gggggccagtgcGCCACCCCATGTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTGG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCAc                                                                                                                                                                                                         <  2:563885/228‑1 (MQ=255)                                   tCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCATAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACCTCCACCCCGTTGAAACCCATCGGCGTGGTACCGCAATCCCACCCCCACAAGTCCTTTCCCCAGCCGAACCCCCCCCCGGCAcc                                                                                                                                                                                                        >  1:276461/1‑239 (MQ=255)                                   tCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGCCATTTCGCCAGCCGACCGCCACGCTGGCAcc                                                                                                                                                                                                        >  2:133577/1‑239 (MQ=255)                                   tCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGGACCGCTATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACCCCGGCAcc                                                                                                                                                                                                        >  2:789186/1‑239 (MQ=255)                                   tCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGGGGTACCGCAATCCCAACCCCAGACGTCATTTTGCCAGCCGCCCGCCCCCCCGGCAcc                                                                                                                                                                                                        >  1:679798/1‑239 (MQ=255)                                   tCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGGGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCCAATCCACCACGTTGAAGCCCATCGGCGTAGTACCGAATTCCCAGCTCCAAACGTCATTTCCCCACCCGACcccccaccgggcacc                                                                                                                                                                                                        >  1:282587/1‑228 (MQ=255)                                     ggCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCCCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCCATATCCCAGCCCCAGACGCCATTTCGCCAGCCGACCGCCACCCTGGcccccg                                                                                                                                                                                                      >  2:838757/1‑237 (MQ=255)                                                 aCCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCCCGCTGGCACCGGAATCCGCCTgg                                                                                                                                                                                           >  2:161949/1‑238 (MQ=255)                                                  ccccATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGGGGTACCGAAAACCCAGCTCCCAACGTCATTTCGCCACCCCACCCCCACCCGCGCACCGGAAACAGACGGGcc                                                                                                                                                                                         >  1:49319/1‑238 (MQ=255)                                                            tCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCACCCGCCCGCCCCCCTGGCCCCGGAACCCGCCGGGCCGCGGTTGcc                                                                                                                                                                                >  1:677285/1‑238 (MQ=255)                                                               gtgtgtttgcgtgggATTC‑TTTTTC‑CCCCACAAAGGCCAGCAAAGAACTCGGAAGTGG‑CCGCTGGGGGGCG‑TTAGCGGTGTAA‑CCCG‑GCGGCCCGATTTCATCGCTGTAACTGATGCCGCCGCCCCCACCCACCACGTTGACGCCCATCGGCGTTGTACCGATACCCCACCTCCATACTTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                                                                                                            <  1:564898/227‑1 (MQ=255)                                                               cGGTTTTGC‑CGGGGGTTT‑TTTTTGCCCCACCAAGGGCCAGGAAAGAACTGGAGAGGGC‑CTGACGGTGGGCG‑TTAGCGGT‑TAA‑CCCAGGCGGCCCGATATCATCGCTGTACCAGATGCCGCCGACCCCATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATACCCCAGCTCCAGCCGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                                                                                                            <  1:362096/239‑1 (MQ=255)                                                               cGGTATTGCTCGGGGATTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                                                                                                            <  2:855045/239‑1 (MQ=255)                                                               cGGTATTGCTCGGGGAGTT‑TTTTTG‑CCCCCCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTTTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCCCATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                                                                                                            <  2:896915/239‑1 (MQ=255)                                                               cGGTATTGCCCGGGGAGTC‑TTTTTG‑CCCACCAAGGGCCAGAAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                                                                                                            <  2:739288/239‑1 (MQ=255)                                                                cgtattggctggggggTCTTTTTTC‑CCCACCAAAGCCCAG‑AAAGAAACGGGGAGGGGCCCGGCGGTGGGCT‑TTAACGG‑GTAACCCCG‑GGGGCCCGATTACCTCGCTTTAACGGAGCCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGGGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTa                                                                                                                                                                            <  1:275996/229‑1 (MQ=255)                                                                  tATGGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTTTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACa                                                                                                                                                                          <  2:348813/238‑1 (MQ=255)                                                                    ttGCCCGGGGGGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGAGGGG‑CCGACGGGGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTc                                                                                                                                                                       <  2:677285/239‑1 (MQ=255)                                                                    gtGCTGGGGGGGTT‑TTTTTA‑CCCCCCAAGGGCCAGAAAAGAACGGGAGATGGC‑CCGAGGGGGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTACCTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTc                                                                                                                                                                       <  1:575583/238‑1 (MQ=255)                                                                       ctgggggCGTC‑TTTTTC‑CCCACCAAAGCACACAAAAGCACTGGAGATGGG‑CCGCGGGGGGGCT‑TTAACGGTGGAA‑CCCC‑GCGCCCCGTAACCCTCGCTGAAACTGATGCCGCCGACCACATCCCCCACGTTGAAGCCCAGCGGCGTGGAACCGATATCCCAGCTCCAGACTTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTg                                                                                                                                                                    <  1:257746/236‑1 (MQ=255)                                                                       cTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGCTATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCCCCCTGGCCCCGGAATCCCACCTGCTGCGGTTTCCGCTACAGTCCTg                                                                                                                                                                    >  2:411342/1‑239 (MQ=255)                                                                        tCGGGGATTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTg                                                                                                                                                                    <  2:9234/238‑1 (MQ=255)                                                                        tCGGGGAGTC‑TTTTTG‑CCCCCCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTg                                                                                                                                                                    <  2:622400/238‑1 (MQ=255)                                                                        gcgggggggcttttttG‑CCCACCAAGGGCGGGTAAAGATCGGGAGTTGGG‑CCGGCGGGCGGGT‑TTAACGGGGTAA‑CCCG‑GGGGTCCGTTATCTTCGCTGTAATTGCTCCCGCCGACCCCTTCCGCCACGTTGATGCCCACCGCCGGGGTACCGATATCCCAGCTCCGGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTg                                                                                                                                                                    <  1:664246/228‑1 (MQ=255)                                                                        cccgggggttttttttC‑CCCCCCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGTCGGTGGGCG‑TTAACGGTATAA‑CCCA‑GCGGCCCGATATCATCGCGGTAACGGAGCCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATACCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTg                                                                                                                                                                    <  1:137056/228‑1 (MQ=255)                                                                        ccGGGGGGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTg                                                                                                                                                                    <  2:97381/237‑1 (MQ=255)                                                                         gggggAGTT‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTa                                                                                                                                                                  <  2:77148/238‑1 (MQ=255)                                                                            ggAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTTGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACCGGCTGCGGTTGCCGCGACAGGCCTGTAATg                                                                                                                                                               >  2:455485/1‑239 (MQ=255)                                                                            ggAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCCCCAGCCGACCGCCCCCCTGGCACCGGAAACCGCCTGGCTGCGGTTGCCGCTACAGGCCTGAAATg                                                                                                                                                               >  2:377637/1‑239 (MQ=255)                                                                                tC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATTTa                                                                                                                                                             >  2:575583/1‑237 (MQ=255)                                                                                          cccACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGGACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGCCTGGCTTCGGTTGCCCCTACAGTCCTGTAATGTACAGGTGcccc                                                                                                                                                   >  2:236445/1‑239 (MQ=255)                                                                                            caccaAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACCCTTGCCCCGGAATCCAACTGGCCCGGGTTGCCGCCACAGTCCTGTAATGGACAGGGGCCCCAg                                                                                                                                                 >  1:890416/1‑239 (MQ=255)                                                                                                    gCCCGCAAAGAA‑TGGGGCTGGG‑CCCGGTGTGGGGGTTTAACGGTATCA‑CCCC‑GCGGCCCGATAACACCGCGGGAACTGACCCCCCCGCCCACCGCCCCCACGTTGAACCCCATCGGCGTGGTACCGATAACCCAGCCCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑c                                                                                                                                           <  1:299902/236‑1 (MQ=255)                                                                                                        ccAAAGAACGGGGGGTGGG‑CCGGCGGGGGGGC‑TTAACGGTGTAC‑CCCA‑GCGGCCAGATCTCCCCGCTGTAACTGGTGCCGCCGACCACACCCCCCACGTTGAAGCCCAGCGGGGTGGTACCGATATCCCAGCGCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCAt                                                                                                                                    <  1:844871/238‑1 (MQ=255)                                                                                                             aagcgggagggggg‑CCCGCGGTGGGCG‑TTAACGGGGTAA‑CCCG‑GCGGCCCGATTACAGCGCTGTAACTGCGGCCCCCGACCACATCCCCCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCAAGACGTCATTCCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGcc                                                                                                                               <  1:121961/228‑1 (MQ=255)                                                                                                               ctgggggggggg‑CAGGCGGTGGGCG‑TTAACGGTGTAG‑CCCA‑GCGGCCCGATATCATCCCTGAAATGGAGGCCGCCGACCCCATCCCCCACGTTGAAGCCCTTCGGGGGGGTACCGATATCCCAGCCCCAACCGTCATTCCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGcc                                                                                                                               <  1:161949/228‑1 (MQ=255)                                                                                                                 ctgggtgggggcccggggTGGGCT‑TTAACGGTTTAA‑CCCG‑GCGGCCCGATTGCCGCGCTGTAGCTGATGCCCCCGACCCCATCCACCACTTTGAACCCCAGCGGCGTGGTACCGATTTCCCAGCCCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATc                                                                                                                            <  1:883121/223‑1 (MQ=255)                                                                                                                 cggggggtgggcgcagggggggCG‑TAACCGGTGAAC‑CCCA‑GCGGCCCTAACTCACCGCTGTAACTGATGCCGCCGCCCCCCTCCCCCACGTGGAAGCCCATCGGCGGGGTACCGATACCCCAGCTCCAGAGGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATc                                                                                                                            <  1:592336/220‑1 (MQ=255)                                                                                                                   gagaTGGG‑CCGACGGTGGGCG‑TTAACGGTTTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGCCCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCCCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCgg                                                                                                                          <  2:420604/238‑1 (MQ=255)                                                                                                                   gagaTGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCCCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGc                                                                                                                         >  2:274759/1‑239 (MQ=255)                                                                                                                       tGGG‑CCGACGGTGGGCT‑TTAACGGTTTAA‑CCCA‑GCGGCCCGATTCTCTCGTTGTAACTGATGCCGCCGACCACATCCCCCCCGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGCCGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCAtt                                                                                                                      <  2:268355/238‑1 (MQ=255)                                                                                                                        ggg‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTc                                                                                                                                                                       >  2:623831/1‑189 (MQ=255)                                                                                                                        ggg‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTc                                                                                                                                                                       <  1:623831/189‑1 (MQ=255)                                                                                                                          g‑ccgccGGTGGGCT‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGt                                                                                                                   <  2:804835/238‑1 (MQ=255)                                                                                                                          g‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTCCAGTCCTGTAATGTACAGGGGGCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTg                                                                                                                  >  2:741375/1‑239 (MQ=255)                                                                                                                                       tg‑TTAAGGGTTTAC‑CCCA‑CCGGCCCTATTGCAGCCGTGAACCGGACTCCCCCGACCACGTCCCCCCTGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCCCCAGGCGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGcc                                                                                                       <  1:579342/237‑1 (MQ=255)                                                                                                                                       cG‑TTAACGGTGTAA‑CCCA‑GCGGCCAGATATCAGCGCTGTAACTGATCCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGGGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGcc                                                                                                       <  2:625760/238‑1 (MQ=255)                                                                                                                                            aaCGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGt                                                                                                  >  2:33749/1‑239 (MQ=255)                                                                                                                                                 tgtAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCCGACGTCATTTCGCCAGCCGACCGCCACTCTGGCACCGGAATCCGCATGGCTGCGGTTGCCCCTCCAGTCCTGTAAGGTACAAGGGCCCCCGTTGT‑CACCCCATTTGCCATCGGCATTAATGGAAAAAATACCgccctccaat                                                                                             >  1:837041/1‑230 (MQ=255)                                                                                                                                                  gtAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTc                                                                                                                                                                       >  2:570645/1‑165 (MQ=255)                                                                                                                                                  gtAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCAGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAAACCGACTGGCTGCGGTTGCCGCTACCGTCCTGTAATGTACAGGGGCCCCAGTTGT‑CATCCCATTTGGCATCGGCATTTGGGGAGAAACAGCCGACGTTCaaat                                                                                            >  1:350494/1‑239 (MQ=255)                                                                                                                                                   tAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGtt                                                                                                                                               >  2:321196/1‑188 (MQ=255)                                                                                                                                                             cGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCCTTAGTGGGGAAACTGCCGCCGTTCATATTGCCGAAAt                                                                                   >  2:368410/1‑239 (MQ=255)                                                                                                                                                             cGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCACACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACGGGCTGCGGTTGCCGCTACCGTCCTGTAATGTACAGGGGCCCCCGTTGT‑CCTCCCACTGACCATCGGCAGTAGGGGGGAAACTGCCGCCGGTCAAATTTAAGaaac                                                                                   >  1:499487/1‑238 (MQ=255)                                                                                                                                                             cGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAAACGTCATTTCGCCAGCCGACCGCCACCCTGGCCCCGGCATCCGACTGGCGGCGGTTTCCCCCCCCGTCCCGTAACGGACAGG‑GCCCCCGGTGTCCCTCCCATTTGCCCACGGCATTAGCGGGGAAACTGCCCACGTTcgaattaccgcaac                                                                                   >  1:2095/1‑226 (MQ=255)                                                                                                                                                                     attcatcCCTGTAACTGATGCCGCCGACCCCACCCACCACGTTGAAGCCCATCGGCGTGGTACCGATACCCCAGCTCCAGCCGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑g                                                                          <  2:494028/236‑1 (MQ=255)                                                                                                                                                                               gTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATc                                                                >  2:614404/1‑238 (MQ=255)                                                                                                                                                                                gaaCTGATGCCGCCGACCACAACCCCCCCGTTGAAGCCAGTCGGCGTGGTACGGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCg                                                              <  2:875741/238‑1 (MQ=255)                                                                                                                                                                                  acgacacccgccgCCCCCACCCCCCCCGTTGAGCCCCAGCGGCGGGCAACCGATAGCCCAGCCCACGAGCTCATTTCGTCAGCCGACCGCCACGCTGGCCCCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCAGTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGcc                                                            <  1:350105/232‑1 (MQ=255)                                                                                                                                                                                  aCTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTCCAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGAGGTTCCTATTTACGAAATAA‑CTCC‑GAAAAAACAAcccc                                                             >  1:63517/1‑236 (MQ=255)                                                                                                                                                                                     ttaccccccgCCCCCATCCACCACGTTGAAGCCCAACGGCGTGGTCCCGATACCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGt                                                         <  2:514857/232‑1 (MQ=255)                                                                                                                                                                                     gATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTTACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGt                                                         >  2:262757/1‑239 (MQ=255)                                                                                                                                                                                      aTGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTCC‑GACAGacactccaccagcc                                                        >  2:96/1‑226 (MQ=255)                                                                                                                                                                                        gccgccGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAAAATCCGCCCGTCAg                                                      >  2:522514/1‑239 (MQ=255)                                                                                                                                                                                          ccccgacccacgcccccacGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGaa                                                    <  2:56104/224‑1 (MQ=255)                                                                                                                                                                                              gACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGGCCCCCAGTTGT‑CATCCCCTTTGCCATCGGCATTATGGGGGAAAAGGCCGCCGTCCATATTGACGAAA‑CA‑CCGCGGAAAGAAAAACCCCCCgcacacaaaag                                                >  1:69798/1‑229 (MQ=255)                                                                                                                                                                                              accccccccccccctgTTGACGCCCGTGGGGGTGGCACCAATACCCCAGCCCCCGCGGTCTCTTCGCCGGCCGCCCGCCGCGCTGGCACCGGAATCCGACGGGCTGCGGTTGCCGCTACAGTCCTTTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACg                                                <  1:236445/224‑1 (MQ=255)                                                                                                                                                                                                    ccccccccccttgcaccccacgggCGTGTTACCTATGCCCCAGCCCCCGCGGTTAGTCCCCCAGCCGCCCGCCACGTTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTC‑TATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATc                                         <  1:377637/217‑1 (MQ=255)                                                                                                                                                                                                    aTCCACCACGTTGAAGCCCACCGGCGTGGTACCGAAAACCCGGCTCCAGACGTCTTTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGGTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCAt                                          <  2:2095/239‑1 (MQ=255)                                                                                                                                                                                                         ccccGTTGAAGCCCATCGGCGTGGTACCGATACCCCAGCTCCAGACGTCATTTCCCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCAcc                                     <  2:679798/236‑1 (MQ=255)                                                                                                                                                                                                         ccccGTTAAAGCCCCTCGGCGTGGTCCCGATACCCAAGCCCCAGGCCTCATTTCGCCAGCCGCCCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCAcc                                     <  2:111479/236‑1 (MQ=255)                                                                                                                                                                                                         ccacTTTGAAGCCCAACGGCGGGGGACCGATCGCCCAGCCCCGGACCTCATTCCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCAGCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCAcc                                     <  2:781133/239‑1 (MQ=255)                                                                                                                                                                                                         cagcGTTGAAGCCCATCGGCGTGGTACCGATAGCCCAGCTCCAGACGTCAGTTGGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCAcc                                     <  2:453969/236‑1 (MQ=255)                                                                                                                                                                                                           ccgttaatcccccccGGTGTGGTCCCGATTCCCCACCCCCCGCCGTAATTTCGCAGGCCGACCGCCGCGCTTGCACCGAAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  1:411342/225‑1 (MQ=255)                                                                                                                                                                                                           ccgtgaaaccccATCGGAGGGGTACCGATACCCCAGCCCCAGACGTCATTCCGCCAGCCGCCCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  2:758474/230‑1 (MQ=255)                                                                                                                                                                                                           ccGTTGAATCCCATCGGCGTGTGACCGATATCCCAGCCCCAGCCGTCATTTCGCCCGCCGACCGCCACCCTGGCACCGGAATCCGACTGGCTGCCGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  2:186512/238‑1 (MQ=255)                                                                                                                                                                                                           aggaacaggCCCATCGGCGTGGTACCGATTTCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  2:867198/231‑1 (MQ=255)                                                                                                                                                                                                           acgggtaacccaatgggtgTGGAACGGAGATCCCGGTCCCAGACGTGTTTTCGCCAGCCACCCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACAACCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  2:129648/221‑1 (MQ=255)                                                                                                                                                                                                           acGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  2:63517/239‑1 (MQ=255)                                                                                                                                                                                                           acGTTGAAGCCCATCGGCGTAGTCCCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  2:282587/239‑1 (MQ=255)                                                                                                                                                                                                            tggtggtgtgggtCGGGGTGGCCCCGATCTCCCGGCCCCAGCCGTCTTTCCGCCAGCCGACCGCCACGCGGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTg                                   <  1:294038/226‑1 (MQ=255)                                                                                                                                                                                                            cGTTGAAGCCCATCGGCGTGGTCCCGATATCCCAGCTCCAGACCTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTGc                                  <  2:890416/239‑1 (MQ=255)                                                                                                                                                                                                              ttcaaGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCCACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGGTTTAGTGGGAAAACTGCCAACGGTCAAATTAGCGAAATAACCTGC‑AAAAGAA‑AACCACCCCACCAAAAACGGAGCAACCACCCGCaa                                >  1:390501/4‑239 (MQ=255)                                                                                                                                                                                                                      ggtCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACAACCGCCCGTCAGAATACGGCGCATCCACCAGCAACATggaa                         >  2:294038/3‑236 (MQ=255)                                                                                                                                                                                                                       aTCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGCAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGGGCATCCACCTGCAACATGGTatt                       >  2:178494/1‑237 (MQ=255)                                                                                                                                                                                                                            cgttgccccGATAACCCACCGCCAGTCGGCATTTCGCCCCCCCTCCGCCCCCCTGGCACCGGAATCCGCCTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTgc                   <  1:741375/231‑1 (MQ=255)                                                                                                                                                                                                                              tGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTATTGGGGAAACTGCCGACGTTCCTATTGACGAAATCA‑CTCC‑GAAAGAACATCCCCCCGTCAGCATACCGGCCCTCCACCTCCAACATgggcgctgtttct                >  1:867198/1‑228 (MQ=255)                                                                                                                                                                                                                                    ccatacccccccTCCGGCCGTCTTTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGGCTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGa          <  1:274759/228‑1 (MQ=255)                                                                                                                                                                                                                                       taccccAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATc        <  2:693554/235‑1 (MQ=255)                                                                                                                                                                                                                                           ccAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTG‑AACATGGGAGGTTGGGCTTTCAGATCGGAGTa  >  2:196834/1‑239 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              ACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTC‑TTTTTG‑CCCACCAAAGGCCAGCAAAGAACTGGAGATGGG‑CCGACGGTGGGCG‑TTAACGGTGTAA‑CCCA‑GCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTGCCGCTACAGTCCTGTAATGTACAGGTGCCCCAGTTGT‑CATCCCATTTGCCATCGGCATTAGTGGAGAAACTGCCGACGTTCATATTGACGAAATCA‑CTGC‑GAAAGAACATCCGCCCGTCAGAATACGGCGCATCCACCTGCAACATGGTAGTGTGCGCTTTCAGATCGGAGTA  >  NC_000913/3686305‑3686772

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