BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,502,985	A→T	5.1%	S94C (AGC→TGC)	speB →	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,502,985	0	A	T	5.1%	361.1 / 5.3	116	S94C (AGC→TGC)	speB	agmatinase
Reads supporting (aligned to +/- strand): ref base A (57/53); new base T (2/4); total (59/57)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 118 total aligned reads displayed.

GGCCACCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGC  >  NC_002947/2502837‑2503133

ggCCACCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCA‑CTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAg                                                                                                                                                        >  3:308715/1‑149 (MQ=255)

   caccaGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTGATCGACTATGGCGACTGTCCCTTCGACTGCg                                                                                                                                                      >  3:310909/1‑149 (MQ=255)

   caccaGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCACGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTACGACAGCg                                                                                                                                                      >  2:168795/1‑149 (MQ=255)

       aGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCg                                                                                                                                                      <  4:123104/145‑1 (MQ=255)

         tAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgc                                                                                                                                                >  3:105007/1‑149 (MQ=255)

            ccGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                              >  1:286511/1‑148 (MQ=255)

            ccGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                              >  8:215527/1‑148 (MQ=255)

             cGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                              >  3:274186/1‑147 (MQ=255)

             cGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgc                                                                                                                                                <  1:195077/145‑1 (MQ=255)

             cGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgc                                                                                                                                                >  2:195077/1‑145 (MQ=255)

              gCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCAACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAg                                                                                                                                         >  2:321809/1‑149 (MQ=255)

               ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGt                                                                                                                                        <  2:63243/149‑1 (MQ=255)

               ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGt                                                                                                                                        >  7:309139/1‑149 (MQ=255)

                cccGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTc                                                                                                                                       <  6:211438/149‑1 (MQ=255)

                 ccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTc                                                                                                                                       >  1:48359/1‑149 (MQ=255)

                 ccGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCTGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCg                                                                                                                                      >  8:132455/1‑149 (MQ=255)

                   ggCGCGCGCTTCGGCCCGCTGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCGGCCACTGGCCGTGTGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCTACTGCGGCTTGCCG‑‑CTGTCgaa                                                                                                                                    <  8:44302/149‑3 (MQ=255)

                   ggCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCA‑CTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTa                                                                                                                                   <  4:308715/149‑1 (MQ=255)

                    gcgcgcgcTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGTCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTAATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCG‑‑CTGTCgaag                                                                                                                                   <  4:310909/149‑4 (MQ=255)

                    gcgcgcgcTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCG‑‑CTGTCgaag                                                                                                                                   <  1:337230/149‑4 (MQ=255)

                      gcgcgcTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCCTCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTAcc                                                                                                                                 >  7:288388/1‑149 (MQ=255)

                           cTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACa                                                                                                                            <  6:142209/149‑1 (MQ=255)

                                    gcgGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCTGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAgg                                                                                                                   <  7:132455/149‑1 (MQ=255)

                                    gcgGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAgg                                                                                                                   <  5:325870/149‑1 (MQ=255)

                                    gcgGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAgg                                                                                                                   <  2:394/149‑1 (MQ=255)

                                         cGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCAt                                                                                                             >  3:145254/1‑149 (MQ=255)

                                         cGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCAt                                                                                                             >  7:301675/1‑149 (MQ=255)

                                          gATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCG‑‑CTGTCgaa                                                                                                                                    <  3:118546/126‑3 (MQ=255)

                                          gATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCG‑‑CTGTCgaa                                                                                                                                    >  4:118546/1‑124 (MQ=255)

                                               gT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCCCCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAg                                                                                                       >  3:38081/1‑149 (MQ=255)

                                               gT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAg                                                                                                       <  2:286511/149‑1 (MQ=255)

                                               gT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAg                                                                                                       >  3:42646/1‑149 (MQ=255)

                                                t‑gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAgg                                                                                                                   >  4:60828/1‑137 (MQ=255)

                                                t‑gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAgg                                                                                                                   <  3:60828/137‑1 (MQ=255)

                                                     gcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACAt                                                                                                  >  7:151352/1‑149 (MQ=255)

                                                      ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACAtt                                                                                                 >  1:329047/1‑149 (MQ=255)

                                                         cGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTg                                                                                              <  4:274186/149‑1 (MQ=255)

                                                          ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCt                                                                                               <  8:309139/147‑1 (MQ=255)

                                                          ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTg                                                                                              >  1:176996/1‑148 (MQ=255)

                                                          ggTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGc                                                                                             <  5:341347/149‑1 (MQ=255)

                                                           gTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCa                                                                                            <  7:215527/149‑1 (MQ=255)

                                                              cagcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑c                                                                                                                 >  7:126277/1‑125 (MQ=255)

                                                              cagcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑c                                                                                                                 <  8:126277/125‑1 (MQ=255)

                                                              cagcagGCCTGGACCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGc                                                                                         >  4:154431/1‑149 (MQ=255)

                                                                      cTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTAcc                                                                                                                                 >  7:96477/1‑102 (MQ=255)

                                                                        gggCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTgcgc                                                                                >  4:71500/1‑148 (MQ=255)

                                                                           cccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATg                                                                            <  2:176996/149‑1 (MQ=255)

                                                                           cccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATg                                                                            <  2:87070/149‑1 (MQ=255)

                                                                           cccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATg                                                                            <  3:147130/149‑1 (MQ=255)

                                                                            ccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGc                                                                           >  7:260068/1‑149 (MQ=255)

                                                                            ccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGc                                                                           <  4:38081/149‑1 (MQ=255)

                                                                            ccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGc                                                                           >  6:160481/1‑149 (MQ=255)

                                                                            ccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGc                                                                           >  4:131870/1‑149 (MQ=255)

                                                                            ccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGAGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGc                                                                           >  5:331938/1‑149 (MQ=255)

                                                                              gCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCt                                                                          >  1:344117/1‑148 (MQ=255)

                                                                               ccACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGc                                                                                                           >  1:275122/1‑114 (MQ=255)

                                                                               ccACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGc                                                                                                           <  2:275122/114‑1 (MQ=255)

                                                                                cACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATg                                                                                                            <  8:150970/112‑1 (MQ=255)

                                                                                  cTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAccc                                                                     <  2:329047/149‑1 (MQ=255)

                                                                                     gCCGTGGGGGTTCTACCGGTTCTACCACCTGGCGTTAATGTATTATGGCTACTGTCCCTTGTACAGCGGAACGCCGATGCGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGAcc                                                                      <  1:107526/147‑1 (MQ=255)

                                                                                     gCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                  <  7:143263/149‑1 (MQ=255)

                                                                                      ccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGGCCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                  <  1:84220/149‑1 (MQ=255)

                                                                                      ccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtgg                                                                 >  6:267133/1‑149 (MQ=255)

                                                                                      ccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtgg                                                                 <  1:168795/149‑1 (MQ=255)

                                                                                        gTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggtg                                                               <  3:71500/149‑1 (MQ=255)

                                                                                         tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggtgg                                                              <  8:301675/149‑1 (MQ=255)

                                                                                         tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggtgg                                                              <  4:145254/149‑1 (MQ=255)

                                                                                           ggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCg                                                            <  5:215552/149‑1 (MQ=255)

                                                                                            gCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGa                                                           <  2:344117/149‑1 (MQ=255)

                                                                                              gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGAc                                                          >  3:351891/1‑148 (MQ=255)

                                                                                               ttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                        <  8:260068/149‑1 (MQ=255)

                                                                                                  gacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACtt                                                     >  5:225046/1‑149 (MQ=255)

                                                                                                  gacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACt                                                      >  4:6794/1‑148 (MQ=255)

                                                                                                   acccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCGTCGa                                                                                                                     <  5:92278/85‑1 (MQ=255)

                                                                                                   acccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCGTCGa                                                                                                                     >  6:92278/1‑85 (MQ=255)

                                                                                                       gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcatca                                                <  5:273691/149‑1 (MQ=255)

                                                                                                       gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcatca                                                <  8:290503/149‑1 (MQ=255)

                                                                                                           gaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTa                                            >  2:204071/1‑149 (MQ=255)

                                                                                                            accaccTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGc                                              >  5:109853/1‑146 (MQ=255)

                                                                                                              caccTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                  >  5:298058/1‑124 (MQ=255)

                                                                                                              caccTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                  <  6:298058/124‑1 (MQ=255)

                                                                                                                 cTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCgct                                      <  5:267133/149‑1 (MQ=255)

                                                                                                                   ggCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCgctgc                                    <  8:151352/149‑1 (MQ=255)

                                                                                                                      ggTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTc                                  >  7:222111/1‑148 (MQ=255)

                                                                                                                        tCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAg                               <  4:95995/149‑1 (MQ=255)

                                                                                                                        tCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAg                               <  3:131870/149‑1 (MQ=255)

                                                                                                                            cGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCt                            >  2:87863/1‑148 (MQ=255)

                                                                                                                              aCTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CGGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                         >  4:312172/1‑149 (MQ=255)

                                                                                                                              aCTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                         >  2:127326/1‑149 (MQ=255)

                                                                                                                                  tGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgccc                     <  4:327143/149‑1 (MQ=255)

                                                                                                                                      gACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCccgccg                 >  5:17162/1‑149 (MQ=255)

                                                                                                                                           tCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATg            >  7:277217/1‑149 (MQ=255)

                                                                                                                                             ccTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc          >  8:247017/1‑149 (MQ=255)

                                                                                                                                             ccTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc          >  4:87688/1‑149 (MQ=255)

                                                                                                                                             ccTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc          >  6:87725/1‑149 (MQ=255)

                                                                                                                                               ttCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                  >  6:246968/1‑91 (MQ=255)

                                                                                                                                               ttCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                  <  5:246968/91‑1 (MQ=255)

                                                                                                                                               ttCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGccc        >  1:353738/1‑149 (MQ=255)

                                                                                                                                                     aGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgc  >  8:273108/1‑149 (MQ=255)

                                                                                                                                                     aGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGgc  >  8:61575/1‑149 (MQ=255)

GGCCACCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGT‑GCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCG‑‑CAGTCGGTACCGGACAGCATCGAGG‑CCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGC  >  NC_002947/2502837‑2503133

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