BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,770	T→C	F140F (TTT→TTC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,770	0	T	C	96.4%	258.2 / ‑0.3	85	F140F (TTT→TTC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (37/44); total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-02

ACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTG  >  NC_000913/71523‑72010

aCGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATtcatt                                                                                                                                                                                                                                                 >  1:799652/1‑251 (MQ=255)

  gCCTGGCTGGCAGGGCTTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATTCTGGATCTGCCGGGGGCAAAATTCATTAc                                                                                                                                                                                                                                               >  1:735934/1‑251 (MQ=255)

    ctggctggCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGATGCTGGATCTGCCGGGGGCAAAATTCATTACGc                                                                                                                                                                                                                                             >  2:663762/1‑251 (MQ=255)

         tggcagggattgtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCACAATTCATTACGCCGAat                                                                                                                                                                                                                                        >  1:830481/25‑251 (MQ=255)

          ggcagggattgtcggttgttagctGGGCGACTGGATCTCCTTATGGCTGGGATGGCGGTTTAAAAAGCCGTTGCATCGCTGGTCATGTCTGAAGAAAAACAAATCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGGCGGGCGCTTTGTGGGCCCGACGCGTCCGCTGGTGCCAATGGGGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAata                                                                                                                                                                                                                                       <  2:216042/228‑1 (MQ=255)

          ggcagggattgtcgggtgtctgctGGGCGACTGGATCTCCTTCGGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAATCATTACGCGATAAAACTGAACAGGCGTTGCATCAACACAGCATGTGCGCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAata                                                                                                                                                                                                                                       <  2:47003/228‑1 (MQ=255)

            cagggattgtcggttgtctgctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTTGATCTGCCGGTGGCAAAATTCATTACGCCGAAtatt                                                                                                                                                                                                                                     <  1:542979/230‑1 (MQ=255)

                    gtcggttgtctgctGGGCGACTGGATCTCCTTCTGGATGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGct                                                                                                                                                                                                                              <  1:276296/237‑1 (MQ=255)

                      cggttgtctgctggtcGACTGGATCTCCTTCGGGCGTGGCTGGCGTTTTAAAA‑GCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGGTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGctgcct                                                                                                                                                                                                                          <  2:479180/236‑1 (MQ=255)

                               gctGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCtgt                                                                                                                                                                                                                     <  1:898311/246‑1 (MQ=255)

                                 ggggcgatgggatctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATGTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTGGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcc                                                                                                                                                                                                                 <  2:315832/233‑1 (MQ=255)

                                  gggcgaccggtgttccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGTGGGTCATTTCGGAAGAAAACCAAAGCATTACTCGATAAAACTGAACATTCGTTGCAGCAACACAGCATGTTCACCTTTCTGGTCGGTCTCTTTGTTGGCCCGACGCGTCCGCTGGTTCAAATGGTGGCGGGAATGCTGGTACTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcc                                                                                                                                                                                                                 <  2:883500/233‑1 (MQ=255)

                                     cGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATCTTATCGGCTGCCTGCTGTGGCCACCGt                                                                                                                                                                                                            >  1:592820/1‑251 (MQ=255)

                                      gattGGATCTCCTTCTGGCTGGGCTGGTTGTTTAAAATGCCGTTGCATCGCTGGGCATTTCTTAAGAAAAACAAAGCATTCCCCGATAAAACTGAACATGCTTTGCATCAACAAAGCATGTTCACCATTCTGGGCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAAGGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgcaaccgtt                                                                                                                                                                                                           <  2:657920/248‑8 (MQ=255)

                                          ggatctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTAc                                                                                                                                                                                                       >  1:285485/9‑251 (MQ=255)

                                            atctccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACtt                                                                                                                                                                                                     >  1:937636/7‑251 (MQ=255)

                                                ccttctgtgtggGGTGGCTTTTTAAAAAGCCGTTGCAACGCTGGTAATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTGGGCCCGACGGGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTGTACTTCCTg                                                                                                                                                                                                 <  2:552719/242‑1 (MQ=255)

                                                ccttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTg                                                                                                                                                                                                 >  1:677797/3‑251 (MQ=255)

                                                  ttctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTAGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGGGGCGGGGATGCTGGGTCTGCCGGGGGCAAAATTCCTTACGCCGATTATTTTCGGCTGCCTTATGGGGCCCCCGTTTTACTTTCTGcc                                                                                                                                                                                               >  2:73009/1‑251 (MQ=255)

                                                   tctGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGc                                                                                                                                                                                                <  1:210389/249‑1 (MQ=255)

                                                           ggctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAATCAATACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGTCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCt                                                                                                                                                                                      <  1:398158/248‑1 (MQ=255)

                                                           ggctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCt                                                                                                                                                                                      <  1:971861/248‑1 (MQ=255)

                                                            gTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTg                                                                                                                                                                                     <  1:308489/251‑1 (MQ=255)

                                                             gcaaccatttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTGCCGGTGGCACCATTCATTACACCGAATATTATCGGCTGCCCGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTgg                                                                                                                                                                                    >  2:505415/8‑251 (MQ=255)

                                                                   ttttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcggg‑cg                                                                                                                                                                             >  1:336791/1‑251 (MQ=255)

                                                                    tttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcggg‑cg                                                                                                                                                                             <  1:691650/250‑1 (MQ=255)

                                                                    tttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGg‑cgc                                                                                                                                                                            >  1:715484/1‑251 (MQ=255)

                                                                         aaaGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGAGTCTGGATCTGCCGGTGGCAAAATTCATTACCCCGAATATTATCGGCCGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGagtcgtgcgggggcgggcg                                                                                                                                                                        >  2:356411/1‑233 (MQ=255)

                                                                          aaGCCGTTGCATCGCTGGTCATTTCTGTAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGcgat                                                                                                                                                                      <  2:677797/251‑1 (MQ=255)

                                                                            ggcGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCTGG‑CGCGGcgatcg                                                                                                                                                                    <  1:585811/249‑1 (MQ=255)

                                                                                    atcccTGGTCATTTTTTAAGAAAAACAAAGCATTACCCGAGAAAACTGAACATGCGTTGCACTACCACAGCATGTTCACCATCCGGGTCGGTCGCTGTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCTCCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑Tcca                                                                                                                                                           <  2:592820/247‑2 (MQ=255)

                                                                                      cGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGGACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGACCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGCTGGCAAAATTCCTTACGCCGAATATTATCGGCTGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGGGGG‑CGCGGGGATCGATAt‑ttccgc                                                                                                                                                         >  2:132311/1‑246 (MQ=255)

                                                                                       acTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGACGCTTTGTTGGCCCGACGCGTCCGCTGGTGCAAATGGTGGCGGGAATGCTGGATATGACGGTAGCAAAATTCATTACGCCGCATATTATAGGCTGCCTTCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcggggctgggggtcgatcattcagcc                                                                                                                                                          >  2:508318/2‑229 (MQ=255)

                                                                                          gggcATTTCGGAAGAAAAAAAAAGCATTACTCGATAAAACGGAACAGGCGTTGCATCAACACAGCATGGTCACCATTCGGGGCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGGGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCg                                                                                                                                                       <  2:285485/246‑1 (MQ=255)

                                                                                          ggTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGGTAT‑TCCCGCCggg                                                                                                                                                     >  1:114978/1‑250 (MQ=255)

                                                                                             cATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGGTCGATAT‑TCCAGCCGGAATg                                                                                                                                                  >  1:232588/1‑251 (MQ=255)

                                                                                              aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGTATGTTGGGTCTGCCGGTGGCCACATTCGTTACGCCGGATATTATCGGCGTCCTGCTGTGGCCACCGTTTTAACTCCTTCCTGGGCATCTGGCGGG‑CGGGGcgagcgggttttcagccgggaagg                                                                                                                                                  >  2:16574/1‑230 (MQ=255)

                                                                                              aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGGGCGATAT‑TCCAACCGGAATgg                                                                                                                                                 >  1:648349/1‑250 (MQ=255)

                                                                                                ttCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTTCCTGGGATTCTGGCGGG‑CGCGGCGATCGGTAT‑TCCAGCCGGGATGc                                                                                                                                                 >  1:635148/1‑249 (MQ=255)

                                                                                                 cccgccgcccttccgatctttactCGATAAAACTGAACATGCGTTGCATCAACCCAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCaga                                                                                                                                              <  2:503087/231‑1 (MQ=255)

                                                                                                     aagaaAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCgg                                                                                                                                          <  1:572563/251‑1 (MQ=255)

                                                                                                       gaaAAACAAAGCATTACGCGATAAAACTTAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGGCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCggcg                                                                                                                                        <  1:132311/251‑3 (MQ=255)

                                                                                                        aaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGGTTCTGGCGGG‑CGCGGCGCTCGCTAT‑TTCAGCCGGAATGCAGAGCggcgg                                                                                                                                       >  1:674943/1‑248 (MQ=255)

                                                                                                        aaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGGTAT‑TCCAGCCGGAATGCAGAGCggggg                                                                                                                                       >  1:891673/1‑248 (MQ=255)

                                                                                                         aaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGCCGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGAGTCTGGCGGG‑CGCGGCGGTCGAATT‑TTCAGCCGGAAATCAGAGCggcggg                                                                                                                                      >  1:97749/1‑247 (MQ=255)

                                                                                                             cAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCggcgcgttt                                                                                                                                   >  1:983021/1‑243 (MQ=255)

                                                                                                                  cattactCGATAAAACTGAACATGCGATTCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTACGCTGGTGCCAATGGTGGCGGGAATGCCGGAACTGCCGGTGGCAAAATTCATTACGCCGAATAATATCGGATGCCCTCTATGGCCACCGTTTTAATTCCTACCTGGGATTCTGGCGGG‑CGCGGAGATCGA‑ATATCCAGCCGGAATGAAGAGCgggaggattaaatgg                                                                                                                             >  1:917214/4‑238 (MQ=255)

                                                                                                                   attactCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGCATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑GGGGGGGGTCGATAT‑TCCCGCCCGAATgacgaggggcgaggtttaatggt                                                                                                                            >  2:307788/3‑229 (MQ=255)

                                                                                                                    ttactCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGCaga                                                                                                                                              >  1:503087/2‑232 (MQ=255)

                                                                                                                     tacccgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGtt                                                                                                                           <  1:920128/245‑1 (MQ=255)

                                                                                                                      actcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACACCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTtgc                                                                                                                         <  1:505415/247‑1 (MQ=255)

                                                                                                                        tcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGCCGAGCGGCGGGTTTAAATgggttctg                                                                                                                       >  1:677428/3‑245 (MQ=255)

                                                                                                                           aTAAACCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                    <  2:336791/251‑1 (MQ=255)

                                                                                                                           aTAAAACTGAACATGCGTTGCATCAACACAGCATTTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                    <  1:191263/251‑1 (MQ=255)

                                                                                                                            aaaaactaaaccggcgctGCATCCACACTGCAGGTGCACCATTCGGGGCGGGCGCTTTGTTGGCCCGACGCTGCCGCGGGTGCCAATGTGGGCGGGAAGGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTgg                                                                                                                   <  2:388327/234‑1 (MQ=255)

                                                                                                                               aaCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGCAGAGCGGGGGGTTTAAATGGTTGCTGCTGGcgg                                                                                                                >  1:708062/1‑249 (MQ=255)

                                                                                                                                   gAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATg                                                                                                                                                  <  1:630070/213‑1 (MQ=255)

                                                                                                                                    aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGGGGGCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGAATGAAGAGCGGCGGGGTTAAATGGTTGGTGCTGGcggaggcg                                                                                                           >  1:142919/1‑244 (MQ=255)

                                                                                                                                       aTGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTAGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGATGGGTATGCCGGTGGCAAAATTAATTAAGCCGAATATTATCGGCTGACTGGTGTGGACACCGTTTTAATTCCCGCCTGGGGTTCTGGGGGGGCGCGG‑GGTCGATAT‑TCCAGGCGGCATCCAGAGCGGGGGGTTTTAATGGTtgctggcggggacggaggc                                                                                                         >  1:415293/1‑236 (MQ=255)

                                                                                                                                         gctttGCATCAACCCAGCAGTTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtat                                                                                                       <  2:114978/247‑3 (MQ=255)

                                                                                                                                         gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGAGAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtatt                                                                                                      <  1:637119/251‑4 (MQ=255)

                                                                                                                                             tGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGCATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATTT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTTAATGGTTGCTGCTggggacggcggtatttgtg                                                                                                  >  1:81655/1‑234 (MQ=255)

                                                                                                                                                        aGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGGATGCAAAGCGGCGAGTTTAAATGGTTTCTTACGGCGACGGCGGTATTTGTGTGGGTttgtgg                                                                                       >  1:611649/1‑251 (MQ=255)

                                                                                                                                                              ttCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGGTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTTGTggctggct                                                                                 >  1:269377/1‑251 (MQ=255)

                                                                                                                                                                      tCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑Tccag                                                                                                                                                          <  1:651795/170‑3 (MQ=255)

                                                                                                                                                                          gacggtcgCTTTGTGGGCCCGCCGCGTCCGCGGGTGCCAATGGTGGCGGGAATGGTGGATCTGCCGGTGGCAAAATTCATTGCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCGGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGtt                                                                     <  2:203310/249‑1 (MQ=255)

                                                                                                                                                                            cggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGGTCGATAT‑TCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGGCGGCGGGATTTGTGTGGGGTGGTGGCTGGCTGTGCTGGCgggttt                                                                   >  1:603635/1‑247 (MQ=255)

                                                                                                                                                                             ggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTGGGGGGTGGCTGTGCTGGCgggtgtt                                                                  >  1:230587/1‑246 (MQ=255)

                                                                                                                                                                                   tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCCGCCGGGATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTTGGTGGGGTGGGGGCTGGCTGtggtgggggggtttgggcgg                                                            >  1:117252/1‑237 (MQ=255)

                                                                                                                                                                                   tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGCGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTTGTGGCTGGCTGTGCTGGCgggtttggggggg                                                            >  1:92557/1‑240 (MQ=255)

                                                                                                                                                                                     tggtgtcccggcgggtcctcgggTGCAAATGGTTGCGGGATTGTTGGATCTGCCGGTGGCAAAATGCATTACGCCGAATATTATCGTCTGCCGGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCTGCGATTGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGTTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                           <  2:64203/229‑1 (MQ=255)

                                                                                                                                                                                     ggtttgcccaccgcgTCCGCTGGTGCCAATGGTGGCGGGAAGGCTGGCTCTGCCGGTGTCCAAATTCACTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGGTTGACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGGTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                           <  2:394733/239‑1 (MQ=255)

                                                                                                                                                                                          gCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TTCCGCCGGGATGCAGAGCGGCGAGGTTAAATGGGTGCTGGTGGCGGCGGCGGTTTTTTTGTGGGGTGGGGGCTGGCGGTGCTTGCgggttttgtggtccgggg                                                       >  2:307855/1‑233 (MQ=255)

                                                                                                                                                                                                 gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                              <  1:473692/251‑6 (MQ=255)

                                                                                                                                                                                                 gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                              <  1:149110/251‑6 (MQ=255)

                                                                                                                                                                                                    cccGCTGGTGCCACTGGTGGCGGGAATGCTGGAGCTGCCGGTGGCAAAAGTCATTTCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccgatc                                           <  2:427747/250‑9 (MQ=255)

                                                                                                                                                                                                                   ggtggCGGGAATGCTGGATCTGCCGGTGGCAAAAGTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTATTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCattat                             <  1:993768/250‑1 (MQ=255)

                                                                                                                                                                                                                         gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                          <  2:118778/167‑1 (MQ=255)

                                                                                                                                                                                                                         gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                          >  1:118778/1‑167 (MQ=255)

                                                                                                                                                                                                                            aaTGCTGGAGCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcgcg                   <  1:865995/251‑1 (MQ=255)

                                                                                                                                                                                                                             aTGCTGGATCTGCAGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCgg                  <  1:783426/251‑1 (MQ=255)

                                                                                                                                                                                                                                        gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTTGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGTCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttg           <  1:242311/247‑1 (MQ=255)

                                                                                                                                                                                                                                        gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttgt          <  1:137733/248‑1 (MQ=255)

                                                                                                                                                                                                                                         ccGGTGGCAAAATTCATTACGCCGAATATTTTCGGCTGCCTGCTGTGGCCACCGTTTTTCTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTtgtg      <  2:648349/251‑1 (MQ=255)

                                                                                                                                                                                                                                          cGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTgg     <  1:4234/251‑1 (MQ=255)

                                                                                                                                                                                                                                              ggcaaaaTTCATTACGACGAAGAGTATCGGCTGGCTGCTGTGGCCATCGTTTTACTGCCGGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTg  <  2:183427/246‑1 (MQ=255)

ACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTG  >  NC_000913/71523‑72010

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