BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	789,959	T→G	13.0%	T7T (ACA→ACC)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,959	0	T	G	13.0%	254.6 / 2.3	105	T7T (ACA→ACC)	galK	galactokinase
Reads supporting (aligned to +/- strand): ref base T (36/52); new base G (14/3); total (50/55)
Fisher's exact test for biased strand distribution p-value = 2.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 108 total aligned reads displayed.

CAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCA  >  NC_000913/789727‑790197

cAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCCGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGt‑ttt                                                                                                                                                                                                                                              <  1:418814/238‑1 (MQ=255)

  aTGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCCCCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGGGTGAGTGGCAGGGTACCCAAATGGGTTGGCAAAAAGAGATTTGGt‑tttttt                                                                                                                                                                                                                                           >  2:175488/1‑238 (MQ=255)

       cgcATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTAATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGGGTTCACCAATAAAATTCACCCGGCCAGGCCC‑CCACATGGGGTGAGTGGCAGGGGGACCAACTGCCTTGGCCAACAGAAGTTGGGt‑tttttttttta                                                                                                                                                                                                                                      >  2:397235/1‑237 (MQ=255)

          aTCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAGTCAAATTCACGCGGCCAGGCGC‑CTGAATGGGTGGGGTGGCGGGGAGGCACAATCGGTTGGCACAAAGAGATTGTGT‑TTTTTTTTTTAgc                                                                                                                                                                                                                                    >  2:789645/1‑237 (MQ=255)

            cGAGGGAAAACTCGTCGGGCTGATTTACATAATCGGCTGCCGTCAGGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAGCCGTCGTTGTAGTCGGTGGGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTc                                                                                                                                                                                                                                 <  1:75083/239‑1 (MQ=255)

            cGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTc                                                                                                                                                                                                                                 <  1:889518/239‑1 (MQ=255)

              aGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAAAAGAGATTGTGT‑TTTTTCTTTCAGACTCAt                                                                                                                                                                                                                               >  1:821699/1‑239 (MQ=255)

               gggAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATAAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGAGTGTTCACCAATCAAATTAACGCGGCCAGGCGC‑CCGAATGGTGGGAGTAGCAGGGGAGCCCAATGAGCTGGCGAGCCAAGATTGGGT‑TTTTTTTTGCAGACACAtt                                                                                                                                                                                                                              >  2:963979/1‑239 (MQ=255)

                                 gATTTTCATAATCGGGTGCCATCAGGGGAACTTTACGGTCATCGCGTGGGGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTcg                                                                                                                                                                                                            <  1:771093/239‑1 (MQ=255)

                                   ttttCATAAGCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTcgcg                                                                                                                                                                                                          <  1:849147/239‑1 (MQ=255)

                                            tCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGat                                                                                                                                                                                                 >  1:377053/1‑239 (MQ=255)

                                            tCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGat                                                                                                                                                                                                 <  1:285777/239‑1 (MQ=255)

                                                          gcgAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTCAGACTCATTTTTTACCCTCCGGATTCGCGAAAACGGATATCGCTGACTgcg                                                                                                                                                                                    >  1:151871/1‑238 (MQ=255)

                                                            gatcttttaggTCATCGGGTGGTTGACAACTTGCCACGTTTTGATAATGAATCGCGGAGGGCAGAACGAAAGCGTCGTTGTAGTCGGTGTTTTCACTCATCAAATTCACGCGGGCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGACAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTgcgcgc                                                                                                                                                                                 <  2:676109/230‑1 (MQ=255)

                                                            gAACTTTACGGTCATGGCGTTGGGCAAAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTgcgcgc                                                                                                                                                                                 <  2:1017339/239‑1 (MQ=255)

                                                              aCTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGCTTGTGT‑TTTTTCTTTCAGACTCATTTCTTGCACTCCGGATTCGCGAAAATGGGAAACGCTGACTGCCGGCaa                                                                                                                                                                               >  1:663411/1‑239 (MQ=255)

                                                                 ttACGGTCATCGCGTGTTCGAACACTGATCACGGTTTGATAATCAAGCGCGCCGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAACCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACg                                                                                                                                                                            <  1:1003366/239‑1 (MQ=255)

                                                                  tACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACg                                                                                                                                                                            >  1:523574/1‑238 (MQ=255)

                                                                   aCGGTCATCGGGTGGTGCACAATTGAACACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGct                                                                                                                                                                          <  1:1008505/239‑1 (MQ=255)

                                                                        cATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACCCTctgc                                                                                                                                                                      >  1:53579/1‑238 (MQ=255)

                                                                                  ttcacaACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTc                                                                                                                                                                                                                                 <  1:898956/167‑1 (MQ=255)

                                                                                      caACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGGCGTTGTAGTCGGTGTGTGCGCCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGTCTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGt                                                                                                                                                      <  1:101781/239‑1 (MQ=255)

                                                                                      aaattgttCATGGTTTTATAATCACTCGCACAGGGCCGAACGAAACCGTGGTTGTATTCGGGTGGTTCCCCAATCAAATTCCCGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTCCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGt                                                                                                                                                      <  2:432444/232‑1 (MQ=255)

                                                                                         cTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCTCGACAATGGATATCGCTGACTTCGCGTAAAACATCTCC‑TGCCCGTTCTGCGGTCA‑g                                                                                                                                                  >  1:251018/1‑239 (MQ=255)

                                                                                           gATCCCGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGt                                                                                                                                                <  1:680962/239‑1 (MQ=255)

                                                                                            aTCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCCAACGCTCTGC‑TGCCTGTTCTCCGGTCA‑gg                                                                                                                                                 >  1:191487/1‑237 (MQ=255)

                                                                                            aTCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGAAATCGCTGACTTCGCGCAAACGCCCTGC‑TGCCTGTTCTGCGGTCA‑GGt                                                                                                                                                >  1:369926/1‑238 (MQ=255)

                                                                                            aTCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGCTTCGCGAAAATGGATATCGCGGACTTCTCGCAACCGCTCCGC‑TGCCTGTTCTGCGGTCA‑GGt                                                                                                                                                >  1:745751/1‑238 (MQ=255)

                                                                                            aTCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATCCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACCCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCCCTGC‑TGCCGGtttgtcggtggggg                                                                                                                                                 >  2:1014206/1‑226 (MQ=255)

                                                                                            aTCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAACCAAATGCGTAGGCCAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTtca                                                                                                                                                                                                                       >  2:898956/1‑169 (MQ=255)

                                                                                            aTCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGAAGCCCAATGCGTTGGCAAACCGAGATTGTGG‑TTTTTCTTTCAGACTCATTTCTTTCACTCCGGATTCGCGAAAATTGATATTCCTGCCTGCGCGCAAACGC‑CTGCGTGCCGGTTcgcggggcaggt                                                                                                                                                 >  2:894228/1‑227 (MQ=255)

                                                                                            aTCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGAGGCAGGGTAGCCAAATGAGTTGGCACACAGAGATTGTGT‑TTTTTCTTTACGACTCATTTCTTACACCCCGGCTTCGCGAAAATGGATATCGCTGCCTGCGCGCAAACGCCCTGC‑TGCCTGTTCTGCGGTCA‑GGt                                                                                                                                                >  1:543982/1‑238 (MQ=255)

                                                                                              cACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTTTTTCAGACTCATTTCTTACCCTCCGGATTCGGGAAAATGGATATAGCTGCATCCGCGCAAACCCTCTGC‑GCCCTGTTTTCCGGGAA‑GGgctc                                                                                                                                             >  2:676529/1‑239 (MQ=255)

                                                                                                 ggttttcgattcaattCCGCAGGGCGGGACGAATCCGTCGTTGGAGTTGGTTGGTCCACCAATCAAATTTCCGCGTCCAGTCACATTGAGTGGTGAGGGT‑GCAGGGTAGCCAAATGGCTTGGCAAACAGAGATTGTGT‑GTTTTCTTTCAGACTCATTTGTTACTCTCCGGAGTCGCGAAAAAGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCt                                                                                                                                          <  2:543982/229‑1 (MQ=255)

                                                                                                 ggTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCt                                                                                                                                          <  1:409612/239‑1 (MQ=255)

                                                                                                      gATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGGGTTGGCAAACAGAGGTTGGGT‑TTTTTCTTTTAGACTCATTTCTTACACTCCGGATTCGTGAAAATGGATATCGCTGACTGCGCGCAACAGCTCCTC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTggg                                                                                                                                       >  1:1063010/1‑237 (MQ=255)

                                                                                                           tcaatcGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTTTg                                                                                                                                 >  1:1596/1‑238 (MQ=255)

                                                                                                              atcGCGCAGGGCAGAACGAAACCGTTGTTGTAGTCGGTGTGTTTACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAg                                                                                                                             <  2:290379/239‑1 (MQ=255)

                                                                                                              atcGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTTTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGATTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAg                                                                                                                             <  1:338745/239‑1 (MQ=255)

                                                                                                               tcGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGCGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTTCCCTCCGGATTCGCGAAAATGGGTATCCGTGACTTCGCGCCAACCCTCTGC‑TGCCTTTTCTGCGGTGA‑GGCCTACCTCGGGCTCCGcccg                                                                                                                             >  2:361815/1‑236 (MQ=255)

                                                                                                               tcGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGCAAATGGATATCGCTGACTGCGCGCAAACGCCCTGC‑TGCCTGTTCTTCGGTCA‑GGGCTCGCTGGGTCTCTGCCAGc                                                                                                                            >  1:400506/1‑239 (MQ=255)

                                                                                                                       ggCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCtgc‑tgc                                                                                                                                                                  <  1:281819/194‑1 (MQ=255)

                                                                                                                           gAACGAAACCGTCGTTGTAGGCGGTGTGTGCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATaacc                                                                                                                 <  1:1014206/238‑1 (MQ=255)

                                                                                                                             aCGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATaacc                                                                                                                 <  1:997959/236‑1 (MQ=255)

                                                                                                                                 aaCCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGCGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAACCGCTCTGC‑TGCCTGTTCTGCGGTCC‑GGTCTCGCTGGGTCTCTGCCAGCCTTTATTaaccaaccctc                                                                                                          >  1:118694/1‑236 (MQ=255)

                                                                                                                                 aaCCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGa                                                                                                                                                                                                                                    >  2:220151/1‑119 (MQ=255)

                                                                                                                                 aaCCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGa                                                                                                                                                                                                                                    <  1:220151/119‑1 (MQ=255)

                                                                                                                                 aaCCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCCCt                                                                                                                                          >  1:563864/1‑207 (MQ=255)

                                                                                                                                 aaCCGTCGTTGTAGTCGGTGTGTTCACCAAGCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCt                                                                                                                                          <  2:563864/207‑1 (MQ=255)

                                                                                                                                      tcgtTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTCGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACCCTCCGGATTCGCGAAAATGGATATCCCTGACTGCGCGCAAACGCTCTGC‑CGCCTGGT‑TTCGGGCACGGTCCCGCTGGGTCTCCGCCAGCCTTTACAAACAAACCAAAAAatt                                                                                                     >  2:428916/1‑239 (MQ=255)

                                                                                                                                      tcgtTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGG‑TTTTTCTTTCAAGCTCATTTCTTACACTCCGGATTCGCGAAAATGGAAATCGCTGACTGCCCGCAAACGCTCCGC‑TGCCTGTTCTGCGGTCA‑GGGTCCCCTGGGTCCCTCCCACCATTTCATAAACAACCAAAAAatt                                                                                                     >  2:925330/1‑239 (MQ=255)

                                                                                                                                      tcgtTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTCAGACTCATTTCTTACCCTCCGGCTTTGCGCAAATGGATATCGCTTACTGCCCGCAAACGCCCTTC‑TGCCCGTTTTGCGGTCA‑GGTCGCGCTGGGGCTTTCCCAGCCTTTCAAAACAAACAAAAAAtat                                                                                                     >  2:837905/1‑237 (MQ=255)

                                                                                                                                      tcgtTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGCCTTCGCGCAAACCCCCTGC‑TGCCTGTTCTGCGGTCA‑GGTCCCCTTGGGTCTCCGCCAGAATTTCATAACCAACCAAAAAatt                                                                                                     >  2:877088/1‑239 (MQ=255)

                                                                                                                                         ttgttggcgttttgtTCACCAATCAAATTCCCGGGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacg                                                                                                  <  2:523574/227‑1 (MQ=255)

                                                                                                                                           gttgtgggTTTGTTCCCCAATCAAATTCCCGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGG‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgta                                                                                                <  2:821699/233‑1 (MQ=255)

                                                                                                                                           gTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGTTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgta                                                                                                <  1:175488/239‑1 (MQ=255)

                                                                                                                                            ttgtcgggtgtttCACCAATCAAATTCAAGGGGTCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTTTGT‑TTTTTCTTTCAGACTCATTTCTTACCCTCCGGTTTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgtac                                                                                               <  2:266081/228‑1 (MQ=255)

                                                                                                                                            tggtgggTTTGTTCCCCAAAAAAATTCCAGCGGGCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTGCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATAGCGCGGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgtac                                                                                               <  2:64788/234‑1 (MQ=255)

                                                                                                                                            tAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGGAGCCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTCAGACTCATTTCTTTCCCTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGGCCTGC‑TGCCTGGTCTGCGGCCG‑GGTCTCGCTGGGTCTCTGCCAGAATTTATTAACCACCCAAAAAttacagaaa                                                                                               >  2:605466/1‑232 (MQ=255)

                                                                                                                                                 ggTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGGAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTAGGGTgg                                                                                          >  1:275848/1‑239 (MQ=255)

                                                                                                                                                      gttacccAATCAAATTTACGGGGGCTGGGGG‑CTGAGTGGGGTGAGTGGCAGGGGAGCCAAATTCGTGTGCAAATGGGGAGTGTTG‑TTTTTCTTTCAGACGCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCGGC‑TGCCTTTTCTGCTGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCATCCATAAATTTACGTACGGTGG‑CGGAg                                                                                    <  2:1053176/234‑1 (MQ=255)

                                                                                                                                                       ttCACCAATCAAATTCACGTGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCCCTGC‑TGCCTGTTCTTCGGTCA‑GGTCTCGCTGGGTCCCTCCCACCATTTCAAAAACAACCATAAATTTACGTAACGGGG‑Cggggg                                                                                   >  2:973261/1‑236 (MQ=255)

                                                                                                                                                          aCCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGtc                                                                                                                                     <  1:833797/187‑1 (MQ=255)

                                                                                                                                                             aaTCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAAAGGATATCGCTGGCTGCGCGCAAACGCCCCGC‑TTCCTGTTCTGCGGGCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCAAAAATTATCGACCGGGGG‑GGGGGCgcagca                                                                             >  2:1045457/1‑239 (MQ=255)

                                                                                                                                                              aTCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTTGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGcagcag                                                                            <  1:421882/239‑1 (MQ=255)

                                                                                                                                                              aTCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGcagcag                                                                            <  1:1045457/239‑1 (MQ=255)

                                                                                                                                                               tCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACCGGGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGCCTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCC‑AGGCCCGTCGGGTCCCCGCACGCATTTACTAACCACCCATAACTTTACGAAGGATGGCCGGAGcacaaagg                                                                            >  2:370523/1‑235 (MQ=255)

                                                                                                                                                               tCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAg                                                                                    >  1:525169/1‑230 (MQ=255)

                                                                                                                                                               tCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAg                                                                                    <  2:525169/230‑1 (MQ=255)

                                                                                                                                                               tCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGGGAAGCaga                                                                           >  1:499658/1‑239 (MQ=255)

                                                                                                                                                               tCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGGGG‑CGGGGCGcagcagg                                                                           >  2:1090610/1‑238 (MQ=255)

                                                                                                                                                               tCAAATTCACGCGGCCAGGCGC‑CTGAATGGCGTGAGCTGCAGGGCAGTCAAATGCGTTGGCAAACAGAGATTGGGT‑TTTTTCTTTACGATTCATTTCTTACACTCCGGATTCGCGAAAATGGAACACCCTGACTGCGCGCAAACGCTCCTT‑CGCCCGCTGTGCGGCCA‑GGTCTTCATGGGTCTATCACAGGATTTTATAACAAACCGCAATTTTACCGACGGGGG‑Cgaggacccccaag                                                                           >  2:756660/1‑227 (MQ=255)

                                                                                                                                                                taaaTTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCagag                                                                          <  1:1090610/238‑1 (MQ=255)

                                                                                                                                                                cAAATTCCCCCGGCCAGGCGC‑CTGAATGGTGGTGGGGGCAGGGTGGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACCCTCCGGATTCGCGAAAATGGATAGCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGGCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCagag                                                                          <  2:607146/239‑1 (MQ=255)

                                                                                                                                                                cAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCATCCATAAATTTACGTACGGTGG‑CGGAGCGCAGCagag                                                                          <  1:605466/239‑1 (MQ=255)

                                                                                                                                                                  attccccccgcccAGGCTC‑CTGAATTGTGTGTGTGGCAGGGTAGCCAAATCCGTTGGCAAACAGAGATGTTGT‑TTTTTCTTTCGGACCCCGTTCTTACCCTCCGGATTCCCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTATGCCGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTTG‑CGGATCGCAGCAGAGGc                                                                        <  2:941065/228‑1 (MQ=255)

                                                                                                                                                                         gcgGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCCAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCCTAAATTTACGTACGGTGG‑CCGTGCGCAGCAGAGGCGGTTaaa                                                                 >  1:218624/1‑239 (MQ=255)

                                                                                                                                                                                      cTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCagag                                                                          <  1:651928/218‑1 (MQ=255)

                                                                                                                                                                                               tttggttgCAGGTTTGCCAATTGCGTTGGGAAACAGAGATT‑TGTGTTTTTCTTTCAGCCGCATTTCTTACCCCCCGGATTCGCGAAAATGGATATCGCGGGCTGCGCGCAAACGCTCTGC‑TGCCTTTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAg                                             <  2:576687/231‑1 (MQ=255)

                                                                                                                                                                                                tgAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTTCCAGCATTTCATAACCCACCATAAATTTACGTCCGGTGG‑CGGGGCGCAGCACAGGGGGGAaaaaacagtgctttaaaaaaacccg                                           >  2:651928/1‑218 (MQ=255)

                                                                                                                                                                                                            tcgCCTAATGCGTTGGCAAACGGAGATTGTGT‑TTTTTCTTCCAGACTCATTTCTACCCCTCCGGATTCCCGAAAACGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCtt                                <  2:151871/236‑1 (MQ=255)

                                                                                                                                                                                                            tAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgta                                                                                                >  1:319581/1‑175 (MQ=255)

                                                                                                                                                                                                            tAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgta                                                                                                <  2:319581/175‑1 (MQ=255)

                                                                                                                                                                                                            tAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCtt                                <  1:373242/238‑1 (MQ=255)

                                                                                                                                                                                                            tAGCCAAATGCGTTGGCAAACAGAGATTGTGG‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGACCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCCGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGGGGGCAGCAGCGGGGGACAAAAAGGCGAGGGCAGATGcccggggttgttcactt                                >  2:957735/1‑223 (MQ=255)

                                                                                                                                                                                                             ggccaaagtggTTGGCAAATAGAGATTGTGT‑TTTTTCTGTCGGTCTCATTGCTTAATCGCCGGATTGGCGAAAATGGATATCGCTGACTGCTCGCAAACGCTCTGC‑TGCCTGTGCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCg                              <  2:839533/229‑1 (MQ=255)

                                                                                                                                                                                                                    tGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGCGG‑CGGAGCGCAGCAGGGGCGGATAAAAGGGCGCGTGCAGCTGCCCGTGTTGATTCTCTTCGCCCTTaa                       >  1:749055/1‑239 (MQ=255)

                                                                                                                                                                                                                       ttttcgaaataggtttgttgt‑TTTTTCTTCTCGACTCTTTTTTTCCCCTCGGGATTCCAGAAAAGGGATATCGCTGGCGGCGCGCAAACGCTGTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGGTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATg                    <  2:377053/221‑1 (MQ=255)

                                                                                                                                                                                                                         gcgacaaCAGTGATTGTGT‑TTTCTGTTTCGGACTGATTTCTTACACCCTCGATTCAGGAAAATGGATGTCGCTGACTGCGACCTAACGCTCTGC‑TGCCTTTTCGGCGGTGA‑GAGCTCGCTGGGTCTCTGCAAGCATTGCAAAACAAACCATCAAGTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGGGCGCGTGCGGCGGCCAGTGTTGATTCTCTTCGCCATTAAATgg                   <  2:908778/233‑1 (MQ=255)

                                                                                                                                                                                                                             aaaaagggtttGTGT‑TTTTTCTTTCGGACTCATTTCTTACACTCCGGATTCGAGAAAATGGATAGCGCGTACTGCGCGCAAACGCTCTGC‑TGCCTGTTCCGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGgcgcg               <  2:812123/229‑1 (MQ=255)

                                                                                                                                                                                                                                 agagATTGTGT‑TTTGTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTg          <  1:361815/239‑1 (MQ=255)

                                                                                                                                                                                                                                    gaTTGTGT‑TTTTTCTTTCAGCCTCATTTCTTACCCCCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGt           <  1:782082/235‑1 (MQ=255)

                                                                                                                                                                                                                                     aTTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAACTTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGAGTCTTTTCCCCCTTAAATGGCGCGCCGTGcac       >  1:841181/1‑236 (MQ=255)

                                                                                                                                                                                                                                     aTTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCCTTAAATGGCGCGCCGTGCCa       >  1:876462/1‑238 (MQ=255)

                                                                                                                                                                                                                                     aTTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGAGCGCGTGCAGCGGCCAGTGTTGATTCCCTTGCCCATTAAAAGGAGCGCCGGGCCa       >  2:49988/1‑238 (MQ=255)

                                                                                                                                                                                                                                     aTTGGGT‑TTTGTGTTTCAGACCCATTTTTTACACTCCAGATTCGGGAAAAGGGAACGCGCTTATCGCGGGCAATCGCTCTGC‑TGCCGGTTTGGCGGTCA‑GGTCGCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGGCGTGCCa       <  2:745751/238‑1 (MQ=255)

                                                                                                                                                                                                                                      ttGTGT‑TTTTTCTTTCAGACTCATTTCTTACCCTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCa       <  1:789645/237‑1 (MQ=255)

                                                                                                                                                                                                                                      ttGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCa       <  1:735534/237‑1 (MQ=255)

                                                                                                                                                                                                                                         tgt‑TTTTTCTTTCCGACTCTTTTCTTACCCTCCGGATTCGCGAAAAGGGATATCGCTGACGTCTCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCACGCTGGTTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCa  <  2:663411/239‑1 (MQ=255)

                                                                                                                                                                                                                                         tgt‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGGGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGcccc  >  1:1042912/1‑238 (MQ=255)

CAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGTTTGATAATCAATCGCGCAGGGCAGAACGAAACCGTCGTTGTAGTCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGC‑CTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGT‑TTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGC‑TGCCTGTTCTGCGGTCA‑GGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTACGGTGG‑CGGAGCGCAGCAGAGGCGGATAAAAGTGCGCGTGCAGCTGCCAGTGTTGATTCTCTTCGCCATTAAATGGCGCGCCGTGCCAGCCCA  >  NC_000913/789727‑790197

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