BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,806	G→A	P152P (CCG→CCA)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,806	0	G	A	93.4%	219.8 / NA	82	P152P (CCG→CCA)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand): ref base G (1/2); major base A (31/40); minor base C (4/1); total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTG  >  NC_000913/71560‑72052

cGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATCTTATCGGCTGCCTGCTGTGGCCACCGt                                                                                                                                                                                                                                                      >  1:592820/1‑251 (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                                                                                                                                                    >  1:118778/1‑167 (MQ=255)

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

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

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

                                                                                                                                                                                              gttctgtgtgtgtcaaaagtcattaCGGCGAATATTATCGGCTGCCGGCTGTGGTATTCGTTTTTCTTCATGACTGCGATTCTGGCGGG‑CGCGTCGATCGATGG‑TCCAGCCGGAATTCAGAGTGGCGAGTTGAAATGTTTGCTGATGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGAGGTAATGCAACCGATCGTTTAAGTCATGATTTGTCGCGCGGCCGttt                                                      <  2:756836/230‑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)

                                                                                                                                                                                                                          cGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGcc                          <  1:553288/251‑1 (MQ=255)

                                                                                                                                                                                                                              gattagCGGCTGCCTGCTGGGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAG‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACTGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTAGGGCGGAGTGGTAAAGCAACCGATCGTTTAAGTCATTATTTATCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCg                      <  1:616216/250‑1 (MQ=255)

                                                                                                                                                                                                                                ttaccGGCTGCCTTCTGTGGCCACCTTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATTGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCTATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCgg                     <  2:830481/246‑1 (MQ=255)

                                                                                                                                                                                                                                 tatCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGc                    <  1:73009/250‑1 (MQ=255)

                                                                                                                                                                                                                                         ggctgtggtgGCCTCGGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgcg           <  2:81655/244‑1 (MQ=255)

                                                                                                                                                                                                                                         gcctgcTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgc            <  2:344918/250‑1 (MQ=255)

                                                                                                                                                                                                                                         gcctgcTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCTGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGtggtgg             <  1:924854/249‑1 (MQ=255)

                                                                                                                                                                                                                                                  ggCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGgtggtg  <  1:841082/251‑1 (MQ=255)

CGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGG‑CGCGGCGATCGATAT‑TCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTG  >  NC_000913/71560‑72052

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