BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,805,234	A→T	18.1%	intergenic (+161/‑588)	dsbG → / → CJ019_04513	Thiol:disulfide interchange protein DsbG/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,805,234	0	A	T	18.1%	138.8 / 5.8	100	intergenic (+161/‑588)	dsbG/CJ019_04513	Thiol:disulfide interchange protein DsbG/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (18/59); new base T (7/8); total (31/69)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

CGCCGTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGT  >  ACB143/4805099‑4805362

cgcCGTGAAGCAGGCGACGTGGGGCATGGCACCGGGTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTGTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑a                                                                                                                                          <  8:243843/136‑1 (MQ=255)

cgcCGTGAAGCAGGCGACGGGGGGCAGGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑a                                                                                                                                          <  7:324656/136‑1 (MQ=255)

 gcCGTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGTACCCCCTGCTACGGGGCTCGCCCAAACCGCTGGGGTTTTCGCCCCGCCCCGGCTAGCAAGGCCAGCGCGGGCCGCCG‑G‑A‑‑g                                                                                                                                       >  3:49061/1‑136 (MQ=255)

 gcCGTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑at                                                                                                                                         >  7:93905/1‑135 (MQ=255)

 gcCGTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑g                                                                                                                                       <  3:285279/136‑1 (MQ=255)

 ccgggcaaGCAGGCGCCGGGGGGCATGGCACCGGCTACGCCGGTGGTCGCGGGCGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCGGCGCGGTCGGCCG‑T‑A‑‑g                                                                                                                                       <  1:206470/130‑1 (MQ=255)

    gTGCAGCAGGCGACGTGGGGCAGGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAg                                                                                                                                    <  1:395366/136‑1 (MQ=255)

       aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGG‑C‑g‑g                                                                                                                              >  6:72325/1‑136 (MQ=255)

       aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑G‑‑GGGG‑C‑g‑g                                                                                                                              >  5:220371/1‑136 (MQ=255)

       aaGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑G‑‑GGAC‑C‑g‑g                                                                                                                              >  6:243655/1‑136 (MQ=255)

        aGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCC‑T‑A‑‑GGCG‑C‑c‑gg                                                                                                                             >  2:116256/1‑136 (MQ=255)

               gACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑AC‑GGGC‑C‑G‑GGTTTAcc                                                                                                                       >  4:346898/1‑136 (MQ=255)

                aCGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑TGGG‑C‑GCGGGTTAccg                                                                                                                      >  1:14677/1‑135 (MQ=38)

                aCGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCG‑CGG‑GGTTTAccc                                                                                                                      >  7:89914/1‑136 (MQ=255)

                aCGTGGGGCATGGCAAATTATACGCCGGTGTTCGCGGGTGACCCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGAGGTCGGCCG‑T‑AC‑GGGG‑G‑G‑GGTTTGccc                                                                                                                      >  3:303070/1‑136 (MQ=255)

                    ggggCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCC‑T‑A‑‑GGGG‑CGG‑GGTTTACCcccg                                                                                                                   >  1:384078/1‑135 (MQ=38)

                          tGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGG‑CGG‑GGGTTACCGCCGAAGGGGc                                                                                                            >  8:304955/1‑136 (MQ=37)

                              aCCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCGCC‑G‑GGTTTACCCCCGCATGGCCCgcc                                                                                                        >  6:162048/1‑134 (MQ=38)

                                   cTACGCCGGTGTTCGCTGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑AGTGGTG‑G‑G‑GGTTTACCCCGGACAGGGGCgccgagg                                                                                                    >  4:234537/1‑130 (MQ=255)

                                   cTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑G‑‑GGGG‑C‑G‑GGTTT‑CCCGCCCAAgggccgg‑gccaggg                                                                                                 >  6:189069/1‑132 (MQ=37)

                                   cTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCCGCGCGGTCGGCCG‑TGG‑‑GGGG‑C‑G‑GGGTTGCCCCCCAAGGGCCCggccagcc                                                                                                   >  6:118411/1‑130 (MQ=37)

                                      cGCCGGTGTTCGCGGGGTAACCCGCTCCTACGGGGCTCGCGCAAACGGCTGGGTTTTGCGCCCGGCCCGTGAGCGCACGGCCGGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGc                                                                                             <  5:293775/136‑1 (MQ=37)

                                              ttCGCGGGGGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTc                                                                                      <  7:379339/135‑1 (MQ=255)

                                                  cgGGTTAACCCGCTCCTACGGGGCGCGCGCTAACGGCGGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTccacca                                                                                 <  7:195932/136‑1 (MQ=255)

                                                     gttaccccGCTCCTACGGGGCGCGGGATACCGGGTGGGTTTTGCGCCCGGCCCGTTATGGCACGGCGTGCGGGGTGGGCCT‑T‑T‑‑GGTG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGcc                                                                              <  7:73580/131‑1 (MQ=16)

                                                     gttaaCCCGCTCCTACGGGGCGCGGGCTAAGGGGTGGGTTTTGCGCCTGGCCCGTTATCGCACGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGTG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGcc                                                                              <  6:354757/133‑1 (MQ=21)

                                                                 gtACGGGGCGCGCGCTAAGGGGGGGGTTTTGCGCCCGGCCCGTGATCGCAAGCCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgc                                                                  <  2:21900/135‑1 (MQ=21)

                                                                 cTAGGGGGCTCGCGCTACCGGCTGGGTTTGGCGCCAGGCCCGTGAGCGCACGGCCGGCGCGGGCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgc                                                                  <  8:261215/136‑1 (MQ=21)

                                                                    gggggCGCGCGCTACGGGCTGGGTTTTGCGCCTGCCCGGTTAGCGCACGGCGGGCGGGGTCGGCGG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgc                                                               <  3:355377/135‑1 (MQ=18)

                                                                      gggccggcgcTATCGGGTGGGGTTTGCGCCTGGCCCGTGATCGCATGGCCGGGGGGGCGGGCCT‑T‑A‑‑GTTG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                             <  7:230265/130‑1 (MQ=17)

                                                                      gggCTCGGGCAAACGGCGGGGTTTTGCGCCAGGCCGGTCATCGCACGGCCTGCGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                             <  6:181743/136‑1 (MQ=18)

                                                                      gggCTCGCGCTACCGGCGGGGTTTTGCGCCAGGCCGGTGATCGCAAGGCCTGCGCGGGCGGCCG‑T‑A‑‑GGCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                             <  7:259190/136‑1 (MQ=21)

                                                                      gggCTCGCGCTACCGGCGGGGTTTTGCGCCAGGCCGGTGATCCCACGGCCTGCCGGGTCGGCCG‑T‑A‑‑GGAG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                             <  6:438/136‑1 (MQ=21)

                                                                      gggCGCGCGCTAACGGGTGGGTTTTGCGCCAGGCCCGTGATCGCACGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                             <  7:371438/136‑1 (MQ=255)

                                                                      gggCGCGCCCTAAGGGCGGGGTTTTGCGCCAGGCCCGGTAGCGCACGGCGGGCGGGGTGGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccg                                                             <  8:43353/136‑1 (MQ=18)

                                                                       ggCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑G‑‑GGGGCC‑G‑GGGTTACCC‑CGCAAGGGGCGG‑GCCCGCGGGGCCAATTAACACCACCCGTTCCCCGGTgcccccc                                                             >  5:111645/1‑135 (MQ=255)

                                                                        gcgcgcgCTAACGCCTGGGGTTGGCGCCAGGCCCGTGATCGCAAGGCCGGCGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTgccgccgcc                                                           <  7:128954/133‑1 (MQ=21)

                                                                        gCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGGG‑C‑G‑GGGTTAACCGCGCAAGGGCCGG‑GCCGGC‑GGGGCCAATTCAACCAGCCGGTCACCCGTgccgccccc                                                           >  1:251729/1‑136 (MQ=255)

                                                                            gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCGTT‑G‑‑GGGG‑G‑G‑GGTTTTCCCCGGAAGGGGCCGC‑CCGGGG‑GGGCCAATTTGC‑CCAGCCGGGGACCGGTggcgccggccccg                                                       >  5:80154/1‑130 (MQ=255)

                                                                            gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGGCC‑G‑GGTTTACCGCCCAACGGGCGGG‑CCGGGG‑GGGCCAATTACC‑CCACCCGGGAACCCTTgccgcccccccag                                                       >  5:256827/1‑132 (MQ=255)

                                                                            gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGGCC‑G‑GGTTTACCC‑CGCAATGGGCGG‑CCCAGG‑TGGCCAATTTACCCCACCCGTTCAACGGTgccgcccccccag                                                       >  5:297717/1‑132 (MQ=255)

                                                                            gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCCCCCAAGGGGCGG‑ACCTGC‑TGGCCCATTTACCCCCGCCGTTGCCCGTTGCCCCCCCCCCTg                                                       >  8:273346/1‑136 (MQ=21)

                                                                            gcgcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGGTTACCCCCCAAATGGGCGG‑CCCAGC‑GGGCCAATTTCCCCCCGCCGGTCCCCCTTgccgccgccccag                                                       >  7:307852/1‑132 (MQ=21)

                                                                              gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑TGAG‑CGG‑GGTTTACCC‑CGCAAGGGGCGG‑CCCATC‑GGGGCAATTTTCACCCGGCGGTGACCGGTggcgccgccccagtt                                                     >  7:220033/1‑130 (MQ=255)

                                                                              gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGGCC‑G‑GGTTTACCCCGCAAGGGGCCGG‑CCGGGG‑GGGCCAATTTCCCCCGCCCGTGACCCGTT‑CCGCCGCCCCAGTa                                                     >  2:211403/1‑136 (MQ=255)

                                                                              gcTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCG‑CGG‑GGTTTACCCCGCAAGGGGCCGG‑CCCGCG‑GGGCCAATTGCC‑CCCGCCGGTGACCGGTgccgccgccccagtt                                                     >  4:304833/1‑130 (MQ=255)

                                                                              gcTAACGGCGGGGTTTTGCGCCAGGCCGGTTATCGCACGGCCGGCGCGGTCGGCGG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                     <  2:47952/136‑1 (MQ=18)

                                                                              gacagcGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                     <  2:322679/131‑1 (MQ=255)

                                                                               aaaaCGGCTGGGTTTTGCGCCAGGCCCGTGATCGCACGGCCGGCGGGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                    <  2:58528/134‑1 (MQ=21)

                                                                                         tgTTTTGCGCCAGGCCGTTTTGCGAAAGGCCTGCGCGGTGGGCCG‑T‑A‑‑GGTG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  6:85619/135‑1 (MQ=14)

                                                                                         ggTTTTGGGCCCGGCCCGTGATCCCAAGCCCTGGGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  7:40322/136‑1 (MQ=18)

                                                                                         ggTTTTGCGCCAGGCCGGTGAGCGCAAGGCCGGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  4:287217/136‑1 (MQ=255)

                                                                                         ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  8:89059/136‑1 (MQ=255)

                                                                                         ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  8:312986/136‑1 (MQ=255)

                                                                                         ggTTTTGCGCCAGGCCCGTGATCGCAAGGCCGGCGGGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  2:384078/136‑1 (MQ=255)

                                                                                         ggTTTTGCCCCCGGCAGTTGAGCCCAAGCCCTGCGCGGTGGGCCT‑T‑A‑‑GGAG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  5:219933/136‑1 (MQ=14)

                                                                                         ggTTTTCCGACCGGACCGTGATCGCACGGCCTGCGCGGTCGACCG‑T‑A‑‑GGGG‑C‑G‑GTTTTACCAGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCATCCGTTCACCGTTGACGCCGCCACTGTATTCGC‑GGTTaa                                         <  8:113411/136‑1 (MQ=17)

                                                                                         ggTTTGGCGCCAGGCCCGTGAGCGCAGGCCCTGCGGGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  6:93660/136‑1 (MQ=18)

                                                                                         ggTTTGGCCCCTGCCCCGTTATCCCACGGCCGGCGCGGTGGGCCT‑T‑A‑‑GGCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  2:307684/136‑1 (MQ=17)

                                                                                         ggTTTGCCGCCCGCCCGGTGACCGCACGGCCGGCGGGGTCGGCGG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaa                                         <  6:314782/136‑1 (MQ=17)

                                                                                          gTTTTGCGCCTGGCCCGTGATCGCATGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑G‑G‑GGTTTACCCGCGAAAGGGCGGG‑ACGGGC‑TGGCAAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTaaa                                        <  7:188788/136‑1 (MQ=18)

                                                                                          gTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑C‑GGTTTTCCCCCGAAAGGGCCCG‑GCCGGC‑GGGCCAATTTACACCCGCCCTTGACCCGTGCCGCCGCCACAGTATTAGC‑GGGGaa                                         >  8:312917/1‑135 (MQ=255)

                                                                                             ttGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGGG‑CGG‑GGTTTACCC‑CGCAAGGGGCGGCCCCGGG‑GGGCCAATTGCC‑CCCGCCGGTGACCGGTTGCGCCGCCCCTGTTTTAGCGGGGTAAAcc                                      >  2:241283/1‑136 (MQ=255)

                                                                                               gccccAGGCCGTTTATCCCACGGCCGGCGGGGTCGCCTT‑T‑T‑‑GCCG‑C‑G‑GTTTAACCGGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCg                                    <  7:209060/132‑1 (MQ=14)

                                                                                                 gcCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                     >  8:153932/1‑117 (MQ=255)

                                                                                                 gcCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTa                                                     <  7:153932/117‑1 (MQ=255)

                                                                                                   ccggCCGGTGATCTCACGGCCGGCGGGGCCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                               <  3:329533/134‑1 (MQ=17)

                                                                                                   ccggCCCGTGATCGCAAGGCCGGCGCGGTCGGCGG‑T‑C‑‑GGAG‑C‑G‑GTTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                               <  2:372510/134‑1 (MQ=17)

                                                                                                   cAGGCCCGTGATCGCCCGGCCTGCGCGGGCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                               <  4:306195/136‑1 (MQ=255)

                                                                                                   cAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑T‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                               <  2:292270/136‑1 (MQ=255)

                                                                                                   cAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                               <  8:106918/136‑1 (MQ=255)

                                                                                                   cAGGCCCGTGACCGCAAGCCCGCCGCGGGCGGCGT‑T‑C‑‑GCCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGATCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                               <  5:100578/136‑1 (MQ=17)

                                                                                                   acgccccGTGACCGCACGCCCTGCGGGGCCGGCCT‑T‑A‑‑GGAG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCt                               <  5:113202/132‑1 (MQ=17)

                                                                                                      gCCCGTGATCGCACGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACg                            <  8:195716/136‑1 (MQ=255)

                                                                                                        cggttaTCCCACGGCCTGCGCGTGCGGCGT‑T‑A‑‑GGAT‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCTTTCACCGTTGCCTCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCt                          <  8:157855/131‑1 (MQ=14)

                                                                                                        ccGTGATCGCACGGCCGGCGCGGTCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCt                          <  7:75202/136‑1 (MQ=255)

                                                                                                        ccGTCAGCGCACGGCCGGCGCGGTGGGCGG‑T‑A‑‑GCAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCt                          <  3:278435/136‑1 (MQ=17)

                                                                                                             agcGCACGGCCGGCGCGGGCGCCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCg                     <  4:294835/134‑1 (MQ=14)

                                                                                                             aTCGCACGGCCGGCGGGGGCGGCGG‑T‑A‑‑GCCG‑G‑G‑GTTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATc                      <  2:89742/135‑1 (MQ=17)

                                                                                                              tCGCATGGCGGGCGGGGGGGGCGG‑T‑T‑‑GGCG‑G‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCg                     <  5:196628/135‑1 (MQ=14)

                                                                                                              gcGCACGGCCGGCGGGGGCGGCCT‑T‑T‑‑GGAG‑G‑G‑GGTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCg                     <  6:122572/134‑1 (MQ=17)

                                                                                                               cGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCgg                    <  6:358066/135‑1 (MQ=255)

                                                                                                               cGCAAGGCCTGCGCGGGCGGCGG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCgg                    <  5:179880/135‑1 (MQ=255)

                                                                                                                 catggCCGGCGGGGTGGGCCG‑T‑T‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCggggg                 <  2:6294/133‑1 (MQ=17)

                                                                                                                 cacggCCGGCGGGGGCGGCGG‑A‑A‑‑GTTG‑G‑G‑GGTTACCCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCggggg                 <  4:41124/133‑1 (MQ=17)

                                                                                                                 cAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCggggg                 <  1:233371/136‑1 (MQ=255)

                                                                                                                   tggCCTGCGGGGTCGGCCG‑T‑A‑‑GGTG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTg               <  2:245645/135‑1 (MQ=255)

                                                                                                                   cgggcTGCGCGGTGGGCCG‑T‑A‑‑GGCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTg               <  5:127484/132‑1 (MQ=255)

                                                                                                                   aGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTg               <  1:248870/136‑1 (MQ=255)

                                                                                                                    ggCCTGCGCGGGCGGCCT‑T‑T‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTgg              <  6:161160/136‑1 (MQ=255)

                                                                                                                    gcccTGCGGGGTCGGCCG‑T‑A‑‑GCAG‑G‑G‑GTTTAACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTgg              <  1:214868/134‑1 (MQ=17)

                                                                                                                         gggcgGTCGGCCG‑T‑A‑‑GGCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGaca         <  3:246591/134‑1 (MQ=255)

                                                                                                                          cgcgGACGGCCG‑T‑A‑‑GTAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGaca         <  5:68428/135‑1 (MQ=255)

                                                                                                                            gggTGGGCGG‑T‑A‑‑GCCG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt      <  2:245430/135‑1 (MQ=255)

                                                                                                                            gggTCGGCGG‑A‑A‑‑GTCG‑G‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt      <  2:75684/135‑1 (MQ=17)

                                                                                                                            gggTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt      <  2:111005/135‑1 (MQ=255)

                                                                                                                            cgGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACAt      <  8:33648/136‑1 (MQ=255)

                                                                                                                             ggTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                    >  7:374611/1‑90 (MQ=255)

                                                                                                                             ggTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTAt                                                    <  8:374611/90‑1 (MQ=255)

                                                                                                                                gggCGG‑T‑A‑‑GGTG‑C‑G‑GTTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGTGTGGAGACACATTggg  <  2:100838/135‑2 (MQ=17)

CGCCGTGAAGCAGGCGACGTGGGGCATGGCACCGGCTACGCCGGTGTTCGCGGGTGAACCCGCTCCTACGGGGCTCGCGCTAACGGCTGGGTTTTGCGCCAGGCCCGTGATCGCAAGGCCTGCGCGGTCGGCCG‑T‑A‑‑GGAG‑C‑G‑GGTTTACCCGCGAAAGGGCCGG‑ACCTGC‑TGGCCAATTTCCACCAGCCGTTCACCGTTGCCGCCGCCACTGTATTCGC‑GGGTAAACCCGCTCCTACGCTGATCGGGGGTGGAGACACATTGGT  >  ACB143/4805099‑4805362

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