BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,278,433	T→G	10.7%	T67P (ACC→CCC)	Y75_RS21180 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,278,433	0	T	G	10.7%	567.3 / 6.8	219	T67P (ACC→CCC)	Y75_RS21180	membrane protein
Reads supporting (aligned to +/- strand): ref base T (84/102); new base G (24/8); total (109/110)
Fisher's exact test for biased strand distribution p-value = 2.04e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACA  >  NC_007779/4278196‑4278668

cGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGGTG‑TTGTGCCCG‑ATAACGGTg                                                                                                                                                                                                                                                >  1:2116746/1‑239 (MQ=255)

      cgcCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGTTTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTcc                                                                                                                                                                                                                                          >  1:760962/1‑239 (MQ=255)

      ccccGTCTACCGCGCCGCCGAGGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATAAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTc                                                                                                                                                                                                                                           <  2:525921/236‑1 (MQ=255)

                     cgccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTTTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCa                                                                                                                                                                                                                           >  1:250106/1‑239 (MQ=255)

                       ccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCCACCCGATG‑TTGTGCCCA‑ATAGGGGTGATGTCCGCCATTGTGATGCCaa                                                                                                                                                                                                                          >  1:1339128/1‑238 (MQ=255)

                              acgttGATGATGACATTCCTTTCTAATTTAT‑ATTAAAAGCAATAATTTATATTGTTTT‑AAGGTGCATTATTAACCCATCTGCCAGCCATAAAAAATAAGCCTCCG‑TTATGGAGGCTTCTACTTATCAGTTCAAAATCATTTTTCATTGTGTCGG‑TGGCTTATTCACCT‑CTTCTATCA‑TCTCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCGAGCAGccc                                                                                                                                                                                                                   <  1:365893/235‑1 (MQ=255)

                                tgttatgatggCATTCCGTTCTGATTTGTCATAAAAA‑CAATAGGTAATATTGTTTT‑AAGGTGGATTATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGGCCGCCATTGTGATGCCAAGCAACCCCAt                                                                                                                                                                                                                <  1:904974/235‑1 (MQ=255)

                                  tgatgatgaCATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGGTTAATCT‑TTTCTATCA‑TCAACCACCAGATG‑TTTTGCCAG‑AGAGAGGTGCTGTTCTCCATTTTGATGCCTAACAACCCCaatt                                                                                                                                                                                                              >  1:107333/1‑236 (MQ=255)

                                         gaCATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCCTCT‑CTTCTATCA‑TCACCCACCCGGTG‑TTGTGGCAG‑ATAGCGGGGATTTTCGCCATTGTGGAGCCCAACAACCCCCATGCG‑AACCa                                                                                                                                                                                                      >  2:415998/1‑239 (MQ=255)

                                           tgttCGGTTTTGATTTCTCATTAAAA‑CAATAGATTATATTTTTTA‑AAGGTGGATGATAAAAGCTTCTCCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAAc                                                                                                                                                                                                    <  2:650751/237‑1 (MQ=255)

                                           atttCGGTACTGTTTTGCAATAAAAA‑AAATAGTTTATATTTTTTA‑AAGGGGGAAGTTTAAAGTCTCTCCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTACCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGAGGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCaa                                                                                                                                                                                                     <  2:1651588/236‑1 (MQ=255)

                                                  ttCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTTATCT‑CTTCTATCA‑ACACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCCACCCCCTTGCG‑AACCAACA‑GGCGGa                                                                                                                                                                                            >  1:1859244/1‑239 (MQ=255)

                                                  ttCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTATTATTTATCT‑CTTCTATCA‑TCACCGACCCGATG‑TTTTTTCCG‑ATAGCGGTGATGTCCGCCTTTTTGATGCCCAGCAACCCCATTTCG‑AACCAAAA‑GTCGGa                                                                                                                                                                                            >  1:670829/1‑239 (MQ=255)

                                                  ttCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGGTTGAATTATCT‑CTTTTATAA‑ATAACGGCCCCATG‑TTTGGTCAG‑AAAAGGGGGATTGTGGCCATTGTGATGCCCAGCAAACCCCTTTTG‑AACCAAAA‑GGCGGa                                                                                                                                                                                            >  2:323847/1‑239 (MQ=255)

                                                       aTTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTTTTTCATTCACCCCTTCG‑TG‑TTTTGGCCG‑ATAGCGGGGATGTCCGCCATTGTTTTGCCAATACACCCCCTTGCC‑AACCAACA‑AGCGGGaacaa                                                                                                                                                                                       >  1:1607668/1‑239 (MQ=255)

                                                       aTTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGGCTGATTCATCT‑CTTTTATCA‑TCACCCACTCCCTG‑TGTTGCCCG‑ATAGCCGTGATGTCCCCCCTTTTTATCCCACGCAAACCCCTTCCG‑AACAAACA‑GGCGGAaacaa                                                                                                                                                                                       >  2:711803/1‑239 (MQ=255)

                                                           gTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGcc                                                                                                                                                                                   <  1:1225209/239‑1 (MQ=255)

                                                             aaaaaaaa‑CATTAGGTTAT‑TTGTTTTAAAGGTGGATGATAACACCTTCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGTTTCTACGTATCAGGTCAAAATCATTGGCCTTTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCa                                                                                                                                                                                 <  2:1067674/235‑1 (MQ=255)

                                                                tAAAA‑CAATAGGTTATATTGTTTT‑ATGGTGGTTGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GAc                                                                                                                                                                             <  1:2029292/239‑1 (MQ=255)

                                                                 aaaa‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTTATTCTTCT‑CTTTTATCC‑TCACCGACTCCGTG‑TTGTGCCAG‑AAAACGGTGATGTCCCCCATTTTGATGCCCAGCAACCCCCATGCG‑CACCAACA‑AGCGGGAACAAACCCCA‑AAc                                                                                                                                                                             >  2:1852654/1‑238 (MQ=255)

                                                                         tAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGctgctg                                                                                                                                                                     >  1:1329886/1‑239 (MQ=255)

                                                                           ggTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCaa                                                                                                                                                                                                     <  1:408275/207‑1 (MQ=255)

                                                                            gTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCg                                                                                                                                                                           <  1:585635/230‑1 (MQ=255)

                                                                             ttATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGAt                                                                                                                                                                   <  1:1915624/237‑1 (MQ=255)

                                                                               atatTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTTATCT‑CTTCTATCA‑TCACCGACTTGATG‑TTGTGCCAG‑ATAGCGGGGATGTCCGCCATTGTGATGCCAAACAAACCCATTGCG‑AACCAAAA‑GGCGGAAAAAACGCCCA‑AACCgctc                                                                                                                                                                        >  1:52977/1‑229 (MQ=255)

                                                                                 atTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑GGGGTAGGTCCTCT‑CTTCTAACA‑TCATCCACGCCACGCTTTGGCCAG‑‑CAGCTCTGGTGTCCCCCCTTGTGGTGCCAAGCCAACCCATTTCC‑AAACAACA‑GGCGGAAACAACTCCCC‑GACCGCCGCACATCACCg                                                                                                                                                              >  2:1626618/1‑238 (MQ=255)

                                                                                     tttt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTGTGATTCATCT‑CTTCTATCA‑TCACCGACTCGGAT‑TTGTGGCAG‑ATAGCGGTGATGTCCGCCATTTTGATGCCAAGCAAACCCATTTCG‑GACCAACA‑GGCGGCAAAAACGCCCA‑GACCGCCGCTGATCAcccccac                                                                                                                                                          >  1:1022928/1‑233 (MQ=255)

                                                                                     tttt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TTACCGACTCGATG‑TTGTGCCAG‑ATATCGGTGATGTCCCCCACTTTGATGCCAAACCACCCCACTTCG‑AAACAAAA‑GGCGGAAAAAACGCCCC‑GACCGCTGCTGGTCACCGCAAg                                                                                                                                                          >  1:336959/1‑238 (MQ=255)

                                                                                     tttt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCAACGACTCGATG‑TTGTGCCAG‑ATAGCGGGGATGTTCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACCCCCA‑GACCGCTGCTGATCACCCCAAg                                                                                                                                                          >  1:1380234/1‑238 (MQ=255)

                                                                                     tttt‑AAGGGTGTTGTTTAAGGCATCTCCCCGCCA‑AAAAAAGAAGGCCCCGTTTAGGGAGGTTTTTAGGTATCAGTTAAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGc                                                                                                                                                         <  2:1542677/239‑1 (MQ=255)

                                                                                      ttt‑AAGGTGGTTAATTAAAGCTTTTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTACCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGGCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGc                                                                                                                                                         <  1:711803/238‑1 (MQ=255)

                                                                                       tt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGTATCATCT‑CTTCTATCA‑TCACCCACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAACAAACCCACTGCG‑AACCAAAA‑AGCGGCAACAACGCCCC‑CAACCCTGCCCATCACCCCAAccc                                                                                                                                                        >  2:1915624/1‑238 (MQ=255)

                                                                                       tt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCCATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCCCAAAcc                                                                                                                                                        >  1:57097/1‑238 (MQ=255)

                                                                                       tt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGAATCGATG‑TTGTGCCAG‑ATAGCGGTGCTGTTCGCCATTGTGATGCCCAGCAACCCCATTTCG‑GACCAACA‑GGCGGAAACAACGCCCC‑GAACCCCGCTGATCACCACAAGccc                                                                                                                                                       >  1:1420895/1‑239 (MQ=255)

                                                                                       tt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGGCTTGTTCATCTCCTTCTTTC‑‑TCACCCCACCCATG‑TTGGGCCAGAAAAGCGGTG‑TATTCGGCACTGTGTAGCCCACGAACCCCCTTTCG‑AACCAACA‑GGGGGAAACAACGCCCA‑AACAGCCGCCGATCACCGccaacc                                                                                                                                                        >  2:1406904/1‑233 (MQ=255)

                                                                                       tt‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑GGGCTGATTCATCT‑CTTTTATCA‑TCACCCACCCGATT‑GTGGGGCAG‑ATAGCGGGGATGTCCGCCATTTTGATTCCCAGACACCCCACTGCGAAACCAAAA‑GGCGGAAA‑AACGCCCC‑GGCCCGCGCTCATCTCCGCaaagc                                                                                                                                                        >  2:915511/1‑235 (MQ=255)

                                                                                       tt‑AAGGAGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑GGTATGATTCATCT‑CTTCTTATA‑TCAACGACTCTATG‑GTGATCCAT‑ATAGCGGTGATTTCCGCCATTGTGCTGCCCAACAACCCCATTGCG‑AAACAACA‑AGCGGGAACAACGCCCA‑GACCGCTACTGATCACCGccaacc                                                                                                                                                        >  1:2043414/1‑233 (MQ=255)

                                                                                          aaGGTGTATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑GTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCTTTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCTCCGCAAGCCCGa                                                                                                                                                     <  1:915850/239‑1 (MQ=255)

                                                                                           aGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑GGTCTGATTTATTT‑CTTTTATCA‑TCACCCACTCGATG‑TTTTGCCAG‑ATAGCGGGGGTGTCCCCCCTTTTGGTGCCCAGCAACCCCCTTGCG‑AACCAAAAGGGCGGAAA‑AAAGCCCA‑AACCCCCGCCGAACCCCCCAAAccccc                                                                                                                                                     >  2:1843269/1‑236 (MQ=255)

                                                                                                gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTTTTCC‑CTTCTATCA‑TTAACCACTCGATC‑TTTTGCCAC‑ATTACTGGGCTGTCCGCCCTTGTGATGCCCAACCACCCCGTTGCG‑AACCAACA‑GACGGCAAAAAAACCCACAACCCCGCCGGATC‑CCGCAAGCCCCACAC‑AAt                                                                                                                                              >  2:285393/1‑239 (MQ=255)

                                                                                                gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATTT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCCTTGTGATGCCAAGCAACCCCCTTGGG‑AAACAAAA‑GGCGGAAAAAAAGCCCC‑AAACGCTGCTGATCACCCCAAGCCCCACAT‑aa                                                                                                                                               >  1:623091/1‑238 (MQ=255)

                                                                                                gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTTTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑aa                                                                                                                                               >  1:608389/1‑238 (MQ=255)

                                                                                                gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGAGGTGATGTCCGCCATTGTGATGCCAAACAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCC‑AACCGCTGCTGATCACCGCAAACCCCCCAT‑aa                                                                                                                                               >  1:1746092/1‑238 (MQ=255)

                                                                                                    atTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GAc                                                                                                                                                                             <  1:1477504/205‑1 (MQ=255)

                                                                                                     ttAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGa                                                                                                                                          <  1:1742560/238‑1 (MQ=255)

                                                                                                             tCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTTCGCCATTTTGATGCCAAGCAACCCCCTTTCG‑AACCAACA‑GGCGGAAACAACGCCCC‑GACCGCTGCTGATCACCCCCAACCCGACAA‑AATCAGACCGACGaaa                                                                                                                                 >  1:924455/1‑239 (MQ=255)

                                                                                                                   agccaTAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTAGCGG‑TGTCTTATTCAACT‑CTTCTATCA‑CCACCCACTCGATT‑TGTTGACAC‑ATAACAGCAGTGTACACCATTTAGATGCCCCGCCACCCCCTTCCG‑GACCAACC‑AGCAGGAACAACCCCAC‑AACCGCCACTGATCACCACCAGCCCgcacaaat                                                                                                                                               >  1:1585420/1‑213 (MQ=255)

                                                                                                                     ccaTAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTATGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑GGGTTGGTTCATTT‑CTTCTTTTA‑TCACCCACCCGAAG‑TTGTGGCCG‑ATAGTGGTGGAGTTCGCCACTTTGGTGCCAACCAAACCCATTGCG‑CAACAAAA‑GGCGGAAACAACCCACA‑GAACGCCGCTGCTCCACGCAAGCCCGACAC‑AATCAGAACGCCGCAAACgggagg                                                                                                                         >  2:564891/1‑235 (MQ=255)

                                                                                                                     ccaTAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑CGGATGATTTATCT‑TTTCTTACA‑TCACCGACAAGCTT‑TTGTGCCCG‑ATAACGGTGATGTTCGCCACTGAGCAGCCAAGAAACACCAATGCG‑AACCAACA‑GGAGGAAAAAACGCCCC‑GACCGCCGCTGATCACCGCAACACCCACAT‑ACTCAAACCCACGaaaaacgaagt                                                                                                                         >  2:651369/1‑232 (MQ=255)

                                                                                                                      caTAAAAAAGAAGCCTCCT‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCCCCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGAt                                                                                                                           <  1:935394/236‑1 (MQ=255)

                                                                                                                        tAAAAAAGAAGCCTCCT‑TTATGGAGCCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGa                                                                                                                      <  1:1792413/239‑1 (MQ=255)

                                                                                                                             aagaggCCTCCT‑TTAGGGAGTCTTTTACGTTTCAGGTAAAAACCATTGGCCATTGGGGCGG‑TGTCTGATTCATCT‑CTTCTTTCA‑TCCCCGACTAGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGtgtg                                                                                                                 <  2:1329886/239‑1 (MQ=255)

                                                                                                                               gaagCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTTTCA‑TCACCCACTCGATG‑TTGTGCCCG‑ATAGCGGTGATGTCCGCCATTGTTATGCCAAACAACCCCCTTGCG‑AACCCACA‑GGCCGAAACAACACCCA‑GACCCCCGCCGATCACCCCCAACCCGACATAAATCA‑ACCGACGAAAAAGGATGTTGGgtgt                                                                                                                  >  2:607869/1‑236 (MQ=255)

                                                                                                                               gaagCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGGCTGCTTCATTT‑CTTCTATCA‑TCACCGGCCCGAAG‑TTGTGCCAG‑AGAGCGGGGATGTCCGCCATTGTGAAGCCCAGCAACCCCATTGCG‑AACCAACA‑GGAGGAAACAACGACCA‑CACCGCTGCTCATCCCCGCAAACCCGACAT‑AATCAAACCCACGAAAAAGGAAGGTGAgtgt                                                                                                                  >  2:2042910/1‑236 (MQ=255)

                                                                                                                                  gCCTCCG‑TTGTGGGGGTTTCTACGTTTCAGTTCAAAACCTTGGGCCATTGTGGGGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGCCTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGc                                                                                                            <  2:2047583/239‑1 (MQ=255)

                                                                                                                                      cgg‑TTAGGGAGCGTTTTACGTTTCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGGGCCAG‑ATAGCGGGGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGa                                                                                                        <  2:1746092/237‑1 (MQ=255)

                                                                                                                                            ggggatgtttCTACGTTTCGGGTAAAAATCTTTGGCCATGGGGGGGG‑TGTCTGATTCACCT‑CTTCTATCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATGGCGGTGATGCCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAAcat                                                                                                   <  2:1380234/231‑1 (MQ=255)

                                                                                                                                                  gtttttACGTTTCAGGTCAAAATCATTGGCCTTTGTGGGGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTa                                                                                              <  2:618916/233‑1 (MQ=255)

                                                                                                                                                  gCTTCTCCGTATAAGGTCAAAATCATGGGCCTTTGTGGCGG‑TGTCTAATTCATCT‑CTTCTTCCA‑TCACCGACTCGAGG‑TTGTGCCAG‑ATAGCGGGGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTa                                                                                              <  2:1826054/238‑1 (MQ=255)

                                                                                                                                                  gCTCCTCCGTATCAGGTCAAAATCTTTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTa                                                                                              <  1:650803/238‑1 (MQ=255)

                                                                                                                                                   ttctaaCGTTTCAGTTCAAAATC‑TTGGGCCTTGTTGTGGTTGTCTGATTCATCT‑TTTCTATCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACg                                                                                            <  2:1046145/234‑1 (MQ=255)

                                                                                                                                                    ttCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGa                                                                                                                                          <  1:501931/192‑1 (MQ=255)

                                                                                                                                                       tACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGGGATGTCCGCCATTGTGATGCCAAGCAACCCCCTTGCG‑AACCAACA‑GGCGGAAAAAAAGCCCC‑GACCGCTGCTGATCACCCCAAGCCCGACAT‑AATCAGAACGACGaaaa                                                                                                                                >  1:4436/1‑199 (MQ=255)

                                                                                                                                                        tcGTATCAGGTCAAAATCATTGGCCATGTTGGCGG‑TGTCTGAGTCATCT‑CTTCTATCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACt                                                                                       <  1:836641/238‑1 (MQ=255)

                                                                                                                                                                ggTCAAAATCATTGGCCATTGTGGCGG‑TGTTTGATTTATCT‑CTTCTATCC‑TCCACCACTCGATT‑TTGTGCCAC‑ATAGCGGGGATGTCCGCCCTTGTTATGCCCAACAACCCCCATGCG‑AAACCACA‑GGCGGGAACAACGCCCA‑GACCGCCCCTGATCCACGCAAGCCCGACAT‑AATACAACCGACAAAAACGGAAGGTTAGCGGGCTGGCCAGCAACAAACTTACGGCCCtcaaaaagt                                                                               >  2:823770/1‑231 (MQ=255)

                                                                                                                                                                 gTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACgg                                                                                           <  1:1306499/226‑1 (MQ=255)

                                                                                                                                                                  acAAAATCATTGGCCTTGGGGGAGG‑TGTTGGATTCATCT‑CTTCTTTAA‑TCACCGACTCGATT‑GTGTGCCAG‑ATAGCGGTTAGGTCCGCCATTGGGAGGCCAAGCAGCCCCATAGCG‑AACCAACA‑GGCGTAAACAACGCCCA‑GCCCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGGCGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGc                                                                             <  2:336959/238‑1 (MQ=255)

                                                                                                                                                                     aaaTCATTGGCCTTTGTGGCGG‑TGTCTTATTCATCC‑CTTCTATCA‑CCGCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGAGGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAg                                                                          <  1:514032/239‑1 (MQ=255)

                                                                                                                                                                     aaaTCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCGTTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAg                                                                          <  1:285393/239‑1 (MQ=255)

                                                                                                                                                                          aTTGGCCATTGTGGCGG‑TGGCTGATTCATCT‑CTTCTATCA‑TCACCGCCTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCgaag                                                                     <  1:1582785/239‑1 (MQ=255)

                                                                                                                                                                                cATGGTGGCGG‑TGTCTGATCCACTT‑CTTCTACCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAAc                                                               <  1:1974474/239‑1 (MQ=255)

                                                                                                                                                                                  ttGGGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑CCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCa                                                             <  1:804573/239‑1 (MQ=255)

                                                                                                                                                                                   tgtgGCGG‑TGTCTGATTCATCT‑CTTCTATCT‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGc                                                                                                                                                         <  1:1239457/147‑1 (MQ=255)

                                                                                                                                                                                            tgtCTGATTCATCT‑CTTCTATCA‑TCCCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAgtg                                                     <  1:1843269/238‑1 (MQ=255)

                                                                                                                                                                                                   ttCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGt                                             <  1:1388531/239‑1 (MQ=255)

                                                                                                                                                                                                   ttCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCAAG‑ATAGCGGTGAGGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGt                                             <  1:1704651/239‑1 (MQ=255)

                                                                                                                                                                                                      aTCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACGACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCa                                                        <  1:267489/225‑1 (MQ=255)

                                                                                                                                                                                                      aTCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACGACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCa                                                        >  2:267489/1‑225 (MQ=255)

                                                                                                                                                                                                        ct‑ctTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGCGCCGGCCAGAAACATCATTACGGCAATAAGACGTGGCCAGAGAACAAGAACCACGCCAGTGGTGCGCGTAGCCa                                        >  2:1310569/1‑239 (MQ=255)

                                                                                                                                                                                                             tctATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAAAGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCaa                                    >  2:233184/1‑239 (MQ=255)

                                                                                                                                                                                                               tATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAAc                                  >  1:290948/1‑239 (MQ=255)

                                                                                                                                                                                                                 tca‑tcaCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACaa                                <  1:1647787/239‑1 (MQ=255)

                                                                                                                                                                                                                 tca‑tcaCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACaa                                >  1:773162/1‑239 (MQ=255)

                                                                                                                                                                                                                 tca‑tcaCCGACTCGAGG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACaa                                <  1:1722290/239‑1 (MQ=255)

                                                                                                                                                                                                                       cccGCCTCGATG‑TTGGCCCGG‑ATAGCGGGGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCct                           <  2:773162/238‑1 (MQ=255)

                                                                                                                                                                                                                       aCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTa                                                                                              <  1:1822719/172‑1 (MQ=255)

                                                                                                                                                                                                                       aCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTcc                            >  2:1578415/1‑238 (MQ=255)

                                                                                                                                                                                                                       aCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTcc                            >  1:1929584/1‑238 (MQ=255)

                                                                                                                                                                                                                       aCCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACACCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTcc                            >  2:84593/1‑238 (MQ=255)

                                                                                                                                                                                                                            cTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCCTTACGGCACTAAGAAGTTGCCAGCgaagaagaa                                                                >  1:150652/1‑197 (MQ=255)

                                                                                                                                                                                                                            cTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACg                                                                                            <  2:2074663/169‑1 (MQ=255)

                                                                                                                                                                                                                                 tg‑ttgtGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCg                                                                        <  2:169396/184‑1 (MQ=255)

                                                                                                                                                                                                                                     tgtgCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCg              <  1:4766/239‑1 (MQ=255)

                                                                                                                                                                                                                                     tgtgCCAA‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCg              <  1:1660686/239‑1 (MQ=255)

                                                                                                                                                                                                                                      gtgCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCg              >  2:1183802/1‑238 (MQ=255)

                                                                                                                                                                                                                                      gtgCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGa             >  2:73472/1‑239 (MQ=255)

                                                                                                                                                                                                                                       tgCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCa                                        <  1:2023235/211‑1 (MQ=255)

                                                                                                                                                                                                                                       tgCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGAt            >  1:1180390/1‑239 (MQ=255)

                                                                                                                                                                                                                                            g‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGc          >  2:146891/1‑236 (MQ=255)

                                                                                                                                                                                                                                                    gTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGa    >  1:325287/1‑235 (MQ=255)

CGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAA‑CAATAGGTTATATTGTTTT‑AAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCG‑TTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGG‑TGTCTGATTCATCT‑CTTCTATCA‑TCACCGACTCGATG‑TTGTGCCAG‑ATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCG‑AACCAACA‑GGCGGAAACAACGCCCA‑GACCGCTGCTGATCACCGCAAGCCCGACAT‑AATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACA  >  NC_007779/4278196‑4278668

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