BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,767,251	T→C	9.7%	pseudogene (32/945 nt)	rhsJ →	Rhs‑family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,767,251	0	T	C	9.7%	358.3 / 4.1	150	pseudogene (32/945 nt)	rhsJ	Rhs‑family protein
Reads supporting (aligned to +/- strand): ref base T (108/24); new base C (2/15); total (111/39)
Fisher's exact test for biased strand distribution p-value = 1.62e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 150 total aligned reads displayed.

TTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCA  >  NC_000913/3767015‑3767485

ttttATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGTTAAGCTGAAGCGTTCCTAGCGTGAGTTGCCTATGCCCAGTGGGGGA‑TT‑CCCG‑CCGGCACGG‑GGAGAGTGCCCAACCGGCGGGAAAGGGAAAAcc                                                                                                                                                                                                                                            >  2:653145/1‑239 (MQ=255)

ttttATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGGGCAGA‑TGCTCGACCGGCGGGAAAGTGAAATcc                                                                                                                                                                                                                                            >  2:680749/1‑239 (MQ=255)

          tGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCg                                                                                                                                                                                                                                  >  1:297070/1‑239 (MQ=255)

          tGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCg                                                                                                                                                                                                                                  >  1:807423/1‑239 (MQ=255)

          tGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑GGCTCGACCGGCTGGAACGGGAAACCCTGGCTGACCg                                                                                                                                                                                                                                  >  2:581362/1‑239 (MQ=255)

           ggtttAAAAAAAAATAGTTTTAAGTATGATTACATGGAAGCGGCGGGTGAATTAGGTGATAAACCCCTACTTCCTGTTGTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                                                                                                                                                                                                 <  2:1343430/236‑1 (MQ=255)

                      aaattgtttttatgtttttttttttGATGTTGGGGGTAATTTGGGGGAAAAACCGCTTTTTCTTTTTTTGAATATTTTTTTGTAAAGGTTGGAGGCAAAGGAAATAAAACCTTCCGTTATGGATAGTCGGAAGGGTCCATAGCGTGAGTTGCCTAGGCACAGGGGGGGA‑TC‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGa                                                                                                                                                                                                                      <  2:565092/216‑1 (MQ=12)

                      aaaTACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGTGAGTGAGGa                                                                                                                                                                                                                      >  1:1009136/1‑239 (MQ=255)

                              tttATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGGGCAGA‑TGCTCGACCGGCTGGGAAGTGAAATCCCGGCTGCCCGGGGGAGGGAGGAAAgccgcc                                                                                                                                                                                                              >  1:820729/1‑239 (MQ=255)

                              tttATGATGATATAATTGAAGCTGGGGGTGAATTAGGTGATAAAACGCTACTCCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagtg                                                                                                                                                                                                                          <  1:1195691/227‑1 (MQ=255)

                                 atgatgATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCCCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGGGAGGAAAGCcgccgcg                                                                                                                                                                                                           >  2:1157934/1‑238 (MQ=255)

                                        tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCCCAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCCCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGGGGGGAAAAccgccgggtctttaa                                                                                                                                                                                                    >  2:1195691/1‑230 (MQ=39)

                                        tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCCCGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTCACCGGGTGAGTGAGGAAAGCCGCCgctggctg                                                                                                                                                                                                      >  1:814857/1‑237 (MQ=255)

                                        tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGGGGAAAGCCGCCGCtggctggc                                                                                                                                                                                                    >  1:760158/1‑239 (MQ=255)

                                        tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCtggctggc                                                                                                                                                                                                    >  1:917067/1‑239 (MQ=255)

                                        tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGGGAAATCCTGGCGGACCGGGTGGGGAGGGAAAGCCGCCGCtggctggc                                                                                                                                                                                                    >  1:477274/1‑239 (MQ=255)

                                        tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGCGGGGGA‑TT‑CCCG‑CCGGCCCGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCGGGGTGACCGGGTGAGTGAGGAAAGCAGCCgccgggcgg                                                                                                                                                                                                     >  2:672506/1‑234 (MQ=255)

                                        tatCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGGTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGGGAAATCCGGGGTGACCGGGGGGGTGGGAAAAACCCCCGCcggcgggc                                                                                                                                                                                                    >  2:374121/1‑239 (MQ=255)

                                                      cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                                                                                                                                                                       >  1:883069/1‑238 (MQ=255)

                                                      cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACAGGGGGAGTGAGGAAATCCGCGGCTGGCTGGCATCGGTGGGCCTg                                                                                                                                                                                       >  1:587389/1‑238 (MQ=255)

                                                      cGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCCCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                                                                                                                                                                       >  1:411536/1‑238 (MQ=255)

                                                        ggTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGGAAGTGAAATCCTGGCTGACCGGGGGCGTGAGGAAAGACGCCGCCGGCTGGCATCGTGCGGCCGGAcc                                                                                                                                                                                    >  2:917423/1‑239 (MQ=255)

                                                        ggTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGCGGAAAGCCGCCGCtggctggccccggggtgcccgacc                                                                                                                                                                                    >  2:1038163/1‑223 (MQ=255)

                                                                     aTAAAACCCTCTTTCTTTTTTTAGATATATTTTTGTAAAGTTTTGATGAAAATGAAATTATACTTCCCGTTATGAAAAACCGGAGGGTTTCATGGGGGGATTTCCCTAGCCACGGTGGGGGA‑TC‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGa                                                                                                                                                                       <  2:1021219/239‑1 (MQ=9)

                                                                            gCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATGAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCGGGCCTCGTGCGGCCTGACCGTGGGGGAAATGaaaaaacc                                                                                                                                                                >  2:517790/1‑236 (MQ=255)

                                                                                                 tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtga                                                                                                                                                                                                                             >  2:136417/1‑157 (MQ=255)

                                                                                                 tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGCTGAAAAACCAGAT‑GGACCCGGTGTacctg                                                                                                                                            >  1:710892/1‑234 (MQ=255)

                                                                                                 tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGcc                                                                                                                                          >  1:566322/1‑239 (MQ=255)

                                                                                                 tATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGACAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGGAAGTGAAATCCTGGCTGACCGGGTGAGTTAGGAAAGCCGCGGCTGGCTGGCATCGTGCGGCCTGCCCGTGGGGCAGAGCAAAAAACAAAT‑ggcacaggggtaaaccca                                                                                                                                          >  2:758068/1‑223 (MQ=25)

                                                                                                   tgttgtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGCCTGCCCGGGGGAGGGAGGAAAGCCGCCGCGGGCTAGCATCGGGCGGCCGGACCGTGGGGCAGAGGAAAAACCaacaggaacagggggaacgcacgt                                                                                                                                         >  2:26929/1‑216 (MQ=25)

                                                                                                       ttaagttttttttaaataaaaaaTTATATCCTCCGTTTGTGAGAAGCTGAGTGTTTCAGCGGGTGAGGTCCATAGGCAGAGGGGGGGT‑TT‑CCCG‑CCGGCACGGTGAAGA‑TGCGCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGcgcg                                                                                                                                    <  2:775828/219‑1 (MQ=14)

                                                                                                       gtACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGGGTACAGCCCGGcgcg                                                                                                                                    >  1:906169/1‑239 (MQ=255)

                                                                                                        tACAAGTTTGATGACAATGAAATTGTAACCTACGCTATTGATAAGCTGAAGCGTTCATACCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGCCACGGTACAGA‑TGCTCCACGGGCTGGAAGGTGAAACCCTGGCTGACCGGGCGGGGGTGGCAATCCGCCCCCGGCGGGCATCGGCCGGCCTGACCCTGGGGCGGCAGAAAAACCAGATAGGACCCgagtacccgaccgcggcg                                                                                                                                    >  2:208449/1‑222 (MQ=255)

                                                                                                        tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGGGCCGGTGAACAACCAGAT‑GGACCCGGTGTACACGCCGGcgcgc                                                                                                                                   >  1:784579/1‑238 (MQ=255)

                                                                                                        tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGGGCAGATGAAAAACCAGAT‑GGACGCGGTGTACACGCCGGCGCGa                                                                                                                                   >  1:1087393/1‑239 (MQ=255)

                                                                                                        tACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGCTGAAAAACCAGAT‑GGACCCGGTGTaca                                                                                                                                              >  1:212209/1‑228 (MQ=255)

                                                                                                          cAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGCAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCAG‑CCGGCACGGTGCAGA‑TGCTCGAACGGCTGGAATTTGATATCCTGGCTGACCGGGCGAGGGAGGACAGCCGCCGCTGGCTGGTATCATGCGGCCTGACCGGGGAGAAGCTAAAAAACCAGAT‑GGCCCCGGGGGAAACACCGGCGGGaaa                                                                                                                                 >  2:1081665/1‑239 (MQ=255)

                                                                                                           aaGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATc                                                                                                                                                                                                                                             >  1:148662/1‑131 (MQ=255)

                                                                                                             ttttGATGAAAAGAAAATTAAAATTTCCGCTTTTGATAAGCTAAGGGTTTCATGGGGTAGGTTGCCTAGGCAAAGGGGGGGT‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAAt                                                                                                                              <  1:36085/238‑1 (MQ=17)

                                                                                                             gTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGaaaaa                                                                                                                               >  1:609179/1‑238 (MQ=255)

                                                                                                                  atgaCAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGGGTACACGCCGGCCCGCAAAATccc                                                                                                                           >  1:1157886/1‑236 (MQ=255)

                                                                                                                         tGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACACTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTAcc                                                                                                                   >  1:988911/1‑238 (MQ=255)

                                                                                                                                 taACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGAc                                                                                                           >  1:761741/1‑238 (MQ=255)

                                                                                                                                       ccGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGGGGACACGCCGGCGCGAAAAAACCACCTGAACAACTCGGAcccccggg                                                                                                    >  2:444845/1‑233 (MQ=18)

                                                                                                                                        cGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGCAAGCCGCCGCTGGCTGGCACCGTGCGGCCTGACCGGGGAGCAACTGGAAAACCAGAT‑GGACCCGGGGTCCACGCCGGCGCGAAAAACCAAACTGGACCACTggccccaaccgc                                                                                                    >  2:421144/1‑227 (MQ=18)

                                                                                                                                                    aggggagggttagtggggggggTTGCCGATGCACGGGGGGGGA‑TT‑CCCG‑CCGGCACGGGGCAGA‑TG‑‑CGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcg                                                                                      <  1:794398/217‑1 (MQ=14)

                                                                                                                                                      tGAAGGTTTCATGGGGTGAGTTGCTTATGCCCAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcg                                                                                      <  1:550945/238‑1 (MQ=255)

                                                                                                                                                       gAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCTGCCTGACCGTGGGGCAGATGAAAAACCAGAT‑GGACCCGGGGTGCCCGCCGGCCGGAAAAATCCCCCTGTACCACTGCGACCAACGCGCCCTCCCGCCGgggc                                                                                     >  2:657126/1‑238 (MQ=17)

                                                                                                                                                         aGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGGGTACACGCCGGCGCGAAAAATACACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcgccc                                                                                   >  1:1238507/1‑236 (MQ=255)

                                                                                                                                                             ttCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGGGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCCTCGCGGCCTGCCCCTGGCGCccctca                                                                               >  2:258914/1‑238 (MQ=255)

                                                                                                                                                               cATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACAACTGCGACCATCGCGGCCTGCCCAGGGCGCTCATCAGcc                                                                            >  1:403205/1‑238 (MQ=255)

                                                                                                                                                               cATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGGGTACACGCCGGCGCGAAAAATCCACCTGTACCCCTGCGCCCATCCCGGCCTGCCGCTGGGGCCCATCCGc                                                                             >  1:732364/1‑238 (MQ=255)

                                                                                                                                                                 tggggTGAGTTCCCTATGCACGGGGGGGCA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACg                                                                          <  1:444845/235‑1 (MQ=255)

                                                                                                                                                                   gaggtttttaaaagggacgggggggggA‑TT‑CCCC‑CCGGCCCGGTGCAGA‑TGCTCGACCGGCTGGAAGTGGAAATCCTGGCTGACCGGGTGGGTGAGAAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa                                                                        <  2:501186/217‑1 (MQ=11)

                                                                                                                                                                   aaggagccggcttaggcccggggggggA‑TC‑CCCG‑CCGGCACGGTGCAGA‑TGCCCGACCGGCTGGAAAGTGAAACCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa                                                                        <  1:275544/217‑1 (MQ=14)

                                                                                                                                                                    tggtgtttGCCTATGC‑GAGGGGGGGATTT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAACCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGGGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGa                                                                        <  1:258914/232‑1 (MQ=11)

                                                                                                                                                                     ggggtttccttcggcaaggggggggT‑TC‑CCCG‑GCGGACCGGGGCAGA‑TGCCCGGCCGGGTGGAAAGGGAAATCCGGGCTGGCCGGGGGAGTGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAg                                                                      <  2:815005/219‑1 (MQ=11)

                                                                                                                                                                     gTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGCCCGTGGAGCAGATGAAAAACCAGCT‑GGACCCGGTGTACCCGCCGGCGCGAAAAATCCACCTGTACCCCTGCGCCCATCGCGGCCTCCCCCTGGGGCccacaaacccagaaa                                                                      >  2:1112169/1‑227 (MQ=17)

                                                                                                                                                                      tGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGCT‑GGACCCGGGGTACCAGCCGGCCCGGAAAATCCACCGGAACCACTGCGACCATCGCGGCCGCCCCCTGGGGCCCACCAGCCCGGAAgg                                                                     >  2:917422/1‑239 (MQ=17)

                                                                                                                                                                        aGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGCCCATCGCGGCCTGCCGCTGGCGCCCATCAGCACGGAAggg                                                                    >  1:475206/1‑238 (MQ=255)

                                                                                                                                                                        aGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGc                                                                                                  >  1:1206723/1‑208 (MQ=255)

                                                                                                                                                                        aGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCCCATCAGCCCGGAAgggg                                                                   >  1:530519/1‑239 (MQ=255)

                                                                                                                                                                                gtgcccatggggggT‑TTCCCCG‑CCGGCGGGGGGGGGA‑TCCTCACCGGGGTGGAAAGTAAAATCCTGGCCGACCGGGTGGGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGc                                                           <  2:1279751/230‑1 (MQ=11)

                                                                                                                                                                                    acaGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCAACGCGACCATCGCGGCCTGCCGCTGGCGCTCCTCAGCACGGAAGGGG‑AAACAGccgggt                                                      >  1:607633/1‑234 (MQ=255)

                                                                                                                                                                                       gTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCCCGGAAAATCCACCTGTACCACTGCACCCATCCCGGCCTGCCCCTGGCGGCCCCCACCCCCGCAGGGG‑CAACAGCGGGGTgcg                                                   >  2:749050/1‑239 (MQ=14)

                                                                                                                                                                                        tGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGGGCAGATGAAAAACCAGAT‑GGACCCGGGGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCCCGGCCTCCCCCGGGCGCCCATCCGCACGGAAGGGG‑AAAAAGCggggggggc                                                  >  2:1020569/1‑234 (MQ=14)

                                                                                                                                                                                        tGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGCCCCGGTGTCCACCCCGGCGCGCAACATCCACCTGTACCACTGCGACCATCCCGGCCTGCCCCCGGCGCCCACAACCACGGCAGGGG‑CAACacgggggggcgc                                                  >  2:309907/1‑228 (MQ=14)

                                                                                                                                                                                        tGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTgcgc                                                  >  1:130713/1‑239 (MQ=40)

                                                                                                                                                                                        tGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGAACCACTGCGACCATCGCGGCCTGCCGCTGGGGCCCCTCAGCACGGAAAGGG‑CAAAAACGCGGGgcgc                                                  >  2:288444/1‑239 (MQ=14)

                                                                                                                                                                                        tGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGGACACGCCGGCGCGAAAAATCCACCCGTACAACCGCGACCACCGCGGCCCGCCCCTGGCGCTCACACGCCCGGAAGGGG‑AAACAACggggggggc                                                  >  2:854733/1‑234 (MQ=14)

                                                                                                                                                                                        tGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGGGCACACGCCGGGGGGAAAAATCCACCTGTACCACTGCGACCACCCCGGCCGGCCGCTGCGGCCCACACCCACGGAAGGGG‑Caaaacaggggggcgc                                                  >  2:496841/1‑226 (MQ=14)

                                                                                                                                                                                        tGGGGGA‑TC‑CCCG‑CCGGCACGGTGCAAA‑TGCTCGACCGGTTGGAAAGTGAAATCCTGGCTGACGGGGTGAGTGAGGAAAGCCGGCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTgcg                                                   <  1:1157934/238‑1 (MQ=255)

                                                                                                                                                                                        ggggggA‑TT‑CCCG‑CCGGACCGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTgcg                                                   <  1:979898/237‑1 (MQ=37)

                                                                                                                                                                                         gggggg‑TT‑CCCGCCCGGCAGGGGGGAGG‑GGCTCGACGGGCGGGAAAGGAAAACCCTGGGGGCCGGGGGGGGGGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACCCGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTgcg                                                   <  2:761741/238‑1 (MQ=11)

                                                                                                                                                                                          gggggctgccccc‑ccGCCACGGTGCAGT‑TGCCGGCCCGGCTGAAAAGGGAATCCCGGGCTGCCGGGGAGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTgcgc                                                  <  2:1086504/230‑1 (MQ=11)

                                                                                                                                                                                           gggagtcccccG‑CCGGCGGGGGGAAGG‑TGCTGGGCCGGCTGGAAAGTGAAATCCGGGGGGACCGGGGGGTGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGGGGAGAAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTgcgc                                                  <  2:297070/230‑1 (MQ=11)

                                                                                                                                                                                           gggA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCTCGGCCTGCCGCTGGCGCGCATCAGCACGGAAGGGG‑CACCAGCGTGTGGCGCAGa                                               >  1:1119556/1‑239 (MQ=255)

                                                                                                                                                                                             ggttt‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAACCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGa                                               <  1:668946/235‑1 (MQ=37)

                                                                                                                                                                                             gA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCCCATCAGCACGACAGGGG‑CACAACCGTGGGGCGCAGaaa                                             >  1:1229826/1‑238 (MQ=255)

                                                                                                                                                                                                 t‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAGGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGGGCGAAAAATCCACCTGGACCACTGCGACCATCGCGGCCTGCCGCCGGGCCCCAACACCCCGGAAGGGG‑Caaaaagggggggcgcaaaaaaaa                                          >  2:532748/1‑218 (MQ=11)

                                                                                                                                                                                                 t‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCAACCATCGCGGCCTGCCGCTGCCGCCCCTCAGCCCGGCAGGGG‑Cacaaccggggggggccgaaaagg                                          >  2:1071214/1‑217 (MQ=11)

                                                                                                                                                                                                 t‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGAGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTCCACGCCGACGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGACTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGAgcgccagaaaag                                          >  2:315438/1‑231 (MQ=255)

                                                                                                                                                                                                    ccg‑ccgGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACa                                                             >  1:501186/1‑218 (MQ=37)

                                                                                                                                                                                                    ccg‑ccgGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGAGCGCAGAATaagat                                        >  1:879009/1‑239 (MQ=25)

                                                                                                                                                                                                             cgggTGAAGA‑TGCTCGACCGGCGGGAAAGTAAAACCCTGGGTGGCCGGGTGAGTGAGGAAAGCCGCCGGTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATggg                                 <  1:854733/236‑1 (MQ=255)

                                                                                                                                                                                                                ggggaagtggcccggcggggtGGAAGGGGAACCCCTGCCTGCCCGGTGGGGTGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAGCCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCa                              <  2:607633/219‑1 (MQ=11)

                                                                                                                                                                                                                  gCAGA‑TGCCCGACCGGCTGGAAGGTGAAATCCTGGCTGCACGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAAc                            <  1:1328848/238‑1 (MQ=9)

                                                                                                                                                                                                                        agggcgccgggggggAAAGGGAAACCCTGGCTGGCCGGGTGAGGGGAGAAAGCCGGCGGCGGCTGGGATCGTGCGGCCTGGCCGGGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACCCGCCGGCGCGAAAAATCCACCTGTACCACTGGGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACctgct                       <  1:517790/225‑1 (MQ=11)

                                                                                                                                                                                                                            cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGGGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑GAACAGCGGGGGGCGCAGAAAAAGATAAAATGGGAAACctgctggggg                  >  2:1220114/1‑235 (MQ=11)

                                                                                                                                                                                                                            cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGCTGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCACCGCGGCCTGCCCCTGGCGCTACACAGCACGGAAGGGG‑CAAAACCGTGGGGCGCAGACTATGGAGAAAGGGGGCACctggggaggg                  >  2:215693/1‑232 (MQ=11)

                                                                                                                                                                                                                            cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGAACCACTGCGACCATCCCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑ACACAGCGGGGGGCGCAGAAAATGGTGAATGGGGCAACATGCTGAggg                  >  2:1327813/1‑237 (MQ=11)

                                                                                                                                                                                                                            cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGCAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGCGCAGATGACAAACCAGAT‑GGACCCGGGGTACACACCCGCGCGAAAAATTCACCTGCACCACTGCGACCCTCGCAGCACTCCGCTGGCGACCCTCCGCACGGAAGGCG‑CAACACCGAGGGGCCCCGAATAAGAAAAATTGGAGAACCTCCCGAgag                  >  2:605937/1‑237 (MQ=2)

                                                                                                                                                                                                                               ccGGGTGGAAGGTAAAACCCGGGCGGCCCGGGTGAGGGGGAAAGGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTGCACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGtgat                <  1:657126/238‑1 (MQ=11)

                                                                                                                                                                                                                               ccGGCTGGAAAGGAAAATCCGGGGTGCCCGGGTGGGGGAGGAAAGCCGCCGCTGGCTGGCATCGGGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGtgatg               <  2:1281646/239‑1 (MQ=11)

                                                                                                                                                                                                                                 gggtGGAAAGGGAACTCCTGGCTGACCGGGGGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATgag             <  1:391960/236‑1 (MQ=255)

                                                                                                                                                                                                                                    tGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCCCTGGCGCTCATCAGCCCGGAAGGGGCCAAGAG‑GGGGGGGGCAGAAATGGATGAAAGGGGAAACCGTCTGAGggagggaaa           >  2:391960/1‑232 (MQ=11)

                                                                                                                                                                                                                                        aaGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGc      <  1:236931/239‑1 (MQ=33)

                                                                                                                                                                                                                                         aGTGAAACCCTGGCGGACCGGGGGAGTGAGGAAGGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGca     <  1:1071214/239‑1 (MQ=16)

                                                                                                                                                                                                                                             aaaTCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATatgat                                        >  1:185431/1‑200 (MQ=31)

                                                                                                                                                                                                                                             aaaTCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCACCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CACCAGCGGGGTGCGCAGAGAATGATGAATGGGGCAACctgcgacgggagaagaaaccgcacca  >  2:176892/1‑216 (MQ=11)

TTTTATGCGATGAGTTAAAAAAAAATACTGTTTATGATGATATCATTGAAGCTGCGGGTGAATTAGGTGATAAAACGCTACTTCCTGTTTTAGATACTATGTTGTACAAGTTTGATGACAATGAAATTATAACTTCCGCTATTGATAAGCTGAAGCGTTCATAGCGTGAGTTGCCTATGCACAGTGGGGGA‑TT‑CCCG‑CCGGCACGGTGCAGA‑TGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGAAAAACCAGAT‑GGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTCATCAGCACGGAAGGGG‑CAACAGCGTGGTGCGCAGAATATGATGAATGGGGCAACCTGCTGAGTGATGAGAACCCGCATCA  >  NC_000913/3767015‑3767485

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: `‑`