BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,925,395	C→T	6.9%	L250L (CTG→CTA) ‡	gcvT‑II ←	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,925,395	0	C	T	6.9%	250.9 / 7.1	88	L250L (CTG→CTA) ‡	gcvT‑II	aminomethyltransferase
Reads supporting (aligned to +/- strand): ref base C (41/40); new base T (2/4); total (44/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CACGCACGACCTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGC  >  NC_002947/5925247‑5925538

cacgcacgACCTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGAGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGccc                                                                                                                                                  <  3:55396/149‑1 (MQ=255)

     acgACCTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCCATGCTCCAGCACAGGtt                                                                                                                                             >  3:101322/1‑149 (MQ=255)

         ccTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAg                                                                                                                                        <  1:332995/149‑1 (MQ=255)

           tGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGt                                                                                                                                      >  7:4793/1‑149 (MQ=255)

            gATGGGCGCGCAGTACAACGCGCTCTTCCAGCACCAGACACACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCCACACTGCCGAAGAAGTTGCGCTCGGCCGGAT‑‑CCCAGGCGATGCTGCTGCCCAGGTT‑TGAGGt                                                                                                                                      >  4:238401/1‑148 (MQ=255)

             aTGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCT‑GGAgcatc                                                                                                                                    <  2:54198/149‑5 (MQ=255)

             aTGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCT‑GGAgcatc                                                                                                                                    <  3:221020/149‑5 (MQ=255)

                ggCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCT‑GGAgcatcgc                                                                                                                                  <  1:2463/148‑7 (MQ=255)

                 gcgcgcAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCgg                                                                                                                                <  1:181446/149‑1 (MQ=255)

                      aaGTACACCGAGCTCTTCCAGCACGAGACCCACTAGCTTTTCCTGCACGCCGTGATCGATTTCAGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTGGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGgtg                                                                                                                            <  4:101322/147‑1 (MQ=255)

                       aGTACACCGCGCTCTTCCAGTACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGtgtg                                                                                                                          <  2:159241/149‑1 (MQ=255)

                       aGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGtgtg                                                                                                                          <  5:246946/149‑1 (MQ=255)

                         tACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGt                                                                                                                        >  6:326643/1‑149 (MQ=255)

                            acCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGtttt                                                                                                                     >  4:333845/1‑149 (MQ=255)

                            acCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGtttt                                                                                                                     <  5:326643/149‑1 (MQ=255)

                             ccGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTtc                                                                                                                    >  7:263292/1‑149 (MQ=255)

                             ccGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTtc                                                                                                                    >  1:96895/1‑149 (MQ=255)

                              cgcgCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTtcg                                                                                                                   >  3:172391/1‑149 (MQ=255)

                                 gcTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTtcgtcg                                                                                                                <  3:18082/149‑1 (MQ=255)

                                   tctTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              >  6:299348/1‑149 (MQ=255)

                                     ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              <  3:313805/147‑1 (MQ=255)

                                     ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              <  2:273367/147‑1 (MQ=255)

                                     ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGAt                                                                                                              >  1:273367/1‑147 (MQ=255)

                                     ttCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGt                                                                                                            >  5:263229/1‑149 (MQ=255)

                                         aGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTg                                                                                                        <  4:172391/149‑1 (MQ=255)

                                         aGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCT‑GGAgcatcgc                                                                                                                                  <  4:318915/123‑7 (MQ=255)

                                         aGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAACCTGGGCT‑GGAgcatcgc                                                                                                                                  >  3:318915/1‑117 (MQ=255)

                                          gCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTgg                                                                                                       <  2:96895/149‑1 (MQ=255)

                                            aCCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCACGCAGGCCGACACGGCCGATGAAGTTGCGCTCGCCCGGTT‑‑CCCAGGCGATGCTCCAGCCGAGGTT‑TGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGcc                                                                                                     >  2:316302/1‑149 (MQ=255)

                                              cAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGATGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGATCGGGGTGTGGTTTTCGTCGATGTCCTGGcc                                                                                                     <  6:200357/147‑1 (MQ=255)

                                              cAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGt                                                                                                   <  4:238932/149‑1 (MQ=255)

                                                 aCCCACCAGCTTTTCCTGCACGCCGTGTTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGGGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTa                                                                                                  >  8:151705/1‑147 (MQ=255)

                                                   ccaccaGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAgg                                                                                              >  3:115607/1‑149 (MQ=255)

                                                    caccaGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGt                                                                                             >  6:104704/1‑149 (MQ=255)

                                                         gCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGTTGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATg                                                                                        <  3:11683/149‑1 (MQ=255)

                                                               ccTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCt                                                                                  >  2:118604/1‑149 (MQ=255)

                                                                cTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCtt                                                                                 <  8:263292/149‑1 (MQ=255)

                                                                cTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCt                                                                                  <  1:146470/148‑1 (MQ=255)

                                                                  gCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTc                                                                                >  6:357720/1‑148 (MQ=255)

                                                                    aCGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAg                                                                             <  8:78329/149‑1 (MQ=255)

                                                                     cGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAgg                                                                            <  5:357720/149‑1 (MQ=255)

                                                                      gCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGgc                                                                           >  6:157085/1‑149 (MQ=255)

                                                                          tGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAg                                                                       <  7:350851/149‑1 (MQ=255)

                                                                          tGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAg                                                                       <  3:273915/149‑1 (MQ=255)

                                                                          tGCTCGATTTCCGAAAGCAGGCAGAAACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CACAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGCGCGGGGTGTGGTTTTCGACGATGTCATGGCAGTACAGGTTCATGCCGGCTTCCAGGCGCAg                                                                       >  6:95624/1‑149 (MQ=255)

                                                                             tCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCTCGCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGt                                                                    <  3:200184/149‑1 (MQ=255)

                                                                              cGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACGGGTTCATGCCGGCTTCCAGGCGCAGGGTa                                                                   >  5:333701/1‑149 (MQ=255)

                                                                              cGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGGTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTa                                                                   >  7:343938/1‑149 (MQ=255)

                                                                               gATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTAt                                                                  >  3:32785/1‑149 (MQ=255)

                                                                                     cgaaatcAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCTTGCCGGCTTCCAGGCGCAGGGTATCGCggg                                                            <  3:238401/143‑1 (MQ=255)

                                                                                     cGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCgg                                                             <  1:118604/148‑1 (MQ=255)

                                                                                      gCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgc                                                           >  7:150395/1‑149 (MQ=255)

                                                                                      gCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgc                                                           >  4:144924/1‑149 (MQ=255)

                                                                                      gCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgc                                                           >  4:304204/1‑149 (MQ=255)

                                                                                        aaGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGgcgc                                                         >  5:76734/1‑149 (MQ=255)

                                                                                              gCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTTAGGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAgg                                                    >  1:293465/1‑149 (MQ=255)

                                                                                              gCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGc                                                   <  8:4793/149‑1 (MQ=255)

                                                                                                   cacGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGgcgc                                                                         >  3:312001/1‑122 (MQ=255)

                                                                                                   cacGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGgcgc                                                                         <  4:312001/122‑1 (MQ=255)

                                                                                                    acGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCTGCTTg                                             <  5:299348/149‑1 (MQ=255)

                                                                                                    acGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTg                                             >  6:200107/1‑149 (MQ=255)

                                                                                                     cGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTgg                                            >  4:140353/1‑149 (MQ=255)

                                                                                                      ggCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTggg                                           >  8:133050/1‑149 (MQ=255)

                                                                                                      ggCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGGT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTggg                                           <  3:52753/149‑1 (MQ=255)

                                                                                                       gCCGATGAAGGAGCGCTCGGCCGGTT‑‑CCCAGGCGCTGCTCCAGCCCAGGGT‑GGAGGTGAGCGGGGTGGGGGTTTCGTCGATGTCCTGGCCGTACAGGTTCCTGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCGTgggg                                          <  4:297446/149‑1 (MQ=255)

                                                                                                          gATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTgcggcg                                        >  8:50653/1‑148 (MQ=255)

                                                                                                              atgTTGCGCTGGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGTGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATAc                                    <  5:157085/146‑1 (MQ=255)

                                                                                                               aGTTGCGCTCGGCCGGTA‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCg                                  >  7:171731/1‑149 (MQ=255)

                                                                                                                 ttGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCTGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGgc                                <  1:230581/149‑1 (MQ=255)

                                                                                                                    cgcTCGGCCGGCTCCCCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTATCGTCGATCTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGgc                                <  1:316302/148‑1 (MQ=255)

                                                                                                                       tCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGCTACCGGCGCCGAc                          >  5:287870/1‑149 (MQ=255)

                                                                                                                           ccGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAgg                      <  5:200107/149‑1 (MQ=255)

                                                                                                                           ccGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAgg                      >  2:126982/1‑149 (MQ=255)

                                                                                                                            cGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTgg                                            <  3:233764/126‑1 (MQ=255)

                                                                                                                            cGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGt                     >  4:326145/1‑149 (MQ=255)

                                                                                                                              gTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCg                   >  5:163987/1‑149 (MQ=255)

                                                                                                                                   cccAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGtt                 <  6:163987/148‑1 (MQ=255)

                                                                                                                                   cccAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGtt                 <  6:287870/148‑1 (MQ=255)

                                                                                                                                     cAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTg                >  7:189660/1‑147 (MQ=255)

                                                                                                                                      aGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTgaag             <  3:144924/149‑1 (MQ=255)

                                                                                                                                        gctcaccTCCAACCTGGGCT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTgaag             >  3:102617/7‑147 (MQ=255)

                                                                                                                                        gCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGTTGTCCTGGGCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTTTCGCGGGCGCCGAGGCCGGTTGGGGCGATACCGGCGCCGACCAGGTCGTTgaagaa           <  7:50653/149‑1 (MQ=255)

                                                                                                                                         cGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGaagaag          <  2:293465/149‑1 (MQ=255)

                                                                                                                                           aTGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGc        <  7:133050/149‑1 (MQ=255)

                                                                                                                                             gCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCa      >  5:213992/1‑149 (MQ=255)

                                                                                                                                               tCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGCCCAGGTCGTTGAAGAAGGCCAcc    >  4:238018/1‑149 (MQ=255)

                                                                                                                                               tCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCAcc    >  6:341889/1‑149 (MQ=255)

                                                                                                                                                 cAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  >  3:283833/1‑149 (MQ=255)

                                                                                                                                                 cAGCCCAGGTG‑GGCGGGGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  <  6:263229/149‑1 (MQ=255)

                                                                                                                                                      cAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  >  7:50411/1‑144 (MQ=255)

                                                                                                                                                      cAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGc  <  8:50411/144‑1 (MQ=255)

CACGCACGACCTGATGGGCGCGCAGTACACCGCGCTCTTCCAGCACCAGACCCACCAGCTTTTCCTGCACGCCGTGCTCGATTTCCGCAAGCAGGCCGACACGGCCGATGAAGTTGCGCTCGGCCGGTT‑‑CCCAGGCGATGCTCCAGCCCAGGTT‑GGAGGTGAGCGGGGTGTGGTTTTCGTCGATGTCCTGGCCGTACAGGTTCATGCCGGCTTCCAGGCGCAGGGTATCGCGGGCGCCCAGGCCGCTTGGGGCGATACCGGCGCCGACCAGGTCGTTGAAGAAGGCCACCGC  >  NC_002947/5925247‑5925538

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