BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,797,890	A→T	6.4%	A49A (GCT→GCA)	vdh ←	vanillin dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,797,890	0	A	T	6.4%	299.7 / 8.0	96	A49A (GCT→GCA)	vdh	vanillin dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (46/42); new base T (4/2); total (51/45)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 106 total aligned reads displayed.

CCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATGCT  >  NC_002947/3797744‑3798039

ccATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGgcagca                                                                                                                                                     <  2:1546032/150‑1 (MQ=255)

 cATACCAGTTGGCCATGCCTGTTTTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGcagcag                                                                                                                                                    <  1:126020/150‑1 (MQ=255)

  aTACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCagcagc                                                                                                                                                   <  1:1048672/150‑1 (MQ=255)

   tACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGcc                                                                                                                                                  >  1:788953/1‑150 (MQ=255)

       aGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACcg                                                                                                                                              >  1:1323572/1‑150 (MQ=255)

        gTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACcgc                                                                                                                                             >  2:1548830/1‑150 (MQ=255)

            gCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggca                                                                                                                                         >  1:1508150/1‑150 (MQ=255)

             ccATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcat                                                                                                                                        >  2:2171523/1‑150 (MQ=255)

             ccATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCAGCCGcggcat                                                                                                                                        >  1:1910277/1‑150 (MQ=255)

              cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACcg                                                                                                                                              >  1:2364547/1‑143 (MQ=255)

              cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                       <  2:610130/150‑1 (MQ=255)

              cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                       >  1:1561677/1‑150 (MQ=255)

                 ggCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatcggc                                                                                                                                    <  2:788953/150‑1 (MQ=255)

                   ccccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatcggcat                                                                                                                                  >  1:1038003/1‑150 (MQ=255)

                    cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                       >  1:1841354/1‑144 (MQ=255)

                    cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                       <  2:1841354/144‑1 (MQ=255)

                    cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCggcatcggcatc                                                                                                                                 >  2:1280875/1‑150 (MQ=255)

                    cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCggcatcggcatc                                                                                                                                 >  2:1280257/1‑150 (MQ=255)

                         tttCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCa                                                                                                                            <  2:2204220/150‑1 (MQ=255)

                            cACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAgtg                                                                                                                         <  2:757229/150‑1 (MQ=255)

                            cACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAgtg                                                                                                                         >  1:1672190/1‑150 (MQ=255)

                               cggcggcggCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCagcagc                                                                                                                                                   <  2:1157221/121‑1 (MQ=255)

                               cggcggcggCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCagcagc                                                                                                                                                   >  1:1157221/1‑121 (MQ=255)

                                   ggcggcgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCCAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCGACCGCGGCATCGGCATCGGCCAGTGTggcgg‑c                                                                                                                 <  2:847737/150‑1 (MQ=255)

                                        cgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCNGCCAGCNAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCCGCCGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAg                                                                                                            >  1:1588701/1‑150 (MQ=255)

                                          aTGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGgcagca                                                                                                                                                     <  2:2301166/108‑1 (MQ=255)

                                            gAACTCGGCGGCCGTCGCTTGCAACGGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTg                                                                                                        <  2:1588701/150‑1 (MQ=255)

                                            gAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTg                                                                                                        <  1:689184/150‑1 (MQ=255)

                                             aaCTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGc                                                                                                       <  1:1440466/150‑1 (MQ=255)

                                                tCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGac                                                                                                    >  2:322746/1‑150 (MQ=255)

                                                tCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGac                                                                                                    >  2:1503232/1‑150 (MQ=255)

                                                     ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGGTGGCTGCTGCCAGCGcggc                                                                                                                                          <  2:784166/108‑1 (MQ=255)

                                                     ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGGTGGCTGCTGCCAGCGcggc                                                                                                                                          >  1:784166/1‑108 (MQ=255)

                                                     ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGAcaccac                                                                                               >  2:609003/1‑150 (MQ=255)

                                                     ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGAcaccac                                                                                               <  1:2171523/150‑1 (MQ=255)

                                                        ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGACAAAAGCCGGGCT‑‑CAGCAGCCACCGCGGCCTCGGAATCGGCCTGCATGTCCG‑CGGCAGCCTGCGACACCACCTc                                                                                            >  1:656329/1‑150 (MQ=255)

                                                        ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACTCCACCTc                                                                                            <  1:1118287/150‑1 (MQ=255)

                                                        ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTc                                                                                            >  2:673379/1‑150 (MQ=255)

                                                         cTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTTCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCg                                                                                           >  2:155174/1‑150 (MQ=255)

                                                          tCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCGGCCTGCGACACCACCTCGc                                                                                          >  2:2382611/1‑150 (MQ=255)

                                                            gCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCa                                                                                        <  2:1038003/150‑1 (MQ=255)

                                                               tGCAACAGATGAGCGCCGGCCAGCAAGCGGCTGCGCCGCTCGCCCGGGGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTc                                                                                     <  2:1508150/150‑1 (MQ=255)

                                                               tGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTc                                                                                     >  1:1055398/1‑150 (MQ=255)

                                                                 cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGtgtg                                                                                                                       >  2:2355664/1‑115 (MQ=255)

                                                                 cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGtgtg                                                                                                                       <  1:2355664/115‑1 (MQ=255)

                                                                 cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTggcgg‑cggc                                                                                                              >  2:1777312/1‑123 (MQ=255)

                                                                 cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTggcgg‑cggc                                                                                                              <  1:1777312/123‑1 (MQ=255)

                                                                     aGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCggg                                                                               <  2:1323572/150‑1 (MQ=255)

                                                                     aGATCAGCGCATGCCAGCAACCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCggg                                                                               <  2:552452/150‑1 (MQ=255)

                                                                        tCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTa                                                                            >  1:235662/1‑150 (MQ=255)

                                                                          aGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACa                                                                          <  1:997014/150‑1 (MQ=255)

                                                                          aGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACa                                                                          <  2:1561677/150‑1 (MQ=255)

                                                                           gcgcCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAg                                                                         >  1:1910876/1‑150 (MQ=255)

                                                                            cgcCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAgc                                                                        <  1:1280875/150‑1 (MQ=255)

                                                                               cTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCt                                                                     >  1:1367742/1‑150 (MQ=255)

                                                                                   cAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGaa                                                                 <  1:343575/150‑1 (MQ=255)

                                                                                        aGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTc                                                            <  2:1910876/150‑1 (MQ=255)

                                                                                         gCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCg                                                           >  2:2333520/1‑150 (MQ=255)

                                                                                           ggCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACGTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTc                                                         <  1:1952093/150‑1 (MQ=255)

                                                                                           ggCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTc                                                         >  1:1995867/1‑150 (MQ=255)

                                                                                              tGCGCCGCTCGCCCGGTGCCAGGGGGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCAt                                                      <  1:2320528/150‑1 (MQ=255)

                                                                                              tGCGCCGCTCGCCCGGCGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGa                                                                                                     <  1:1210595/103‑1 (MQ=255)

                                                                                              tGCGCCGCTCGCCCGGCGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGa                                                                                                     >  2:1210595/1‑103 (MQ=255)

                                                                                                 gccgcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGc                                                                                          <  1:2075723/111‑1 (MQ=255)

                                                                                                 gccgcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGc                                                                                          >  2:2075723/1‑111 (MQ=255)

                                                                                                  ccgcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCt                                                  <  2:2255141/150‑1 (MQ=255)

                                                                                                    gcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTgg                                                >  2:1357491/1‑150 (MQ=255)

                                                                                                       cGCCCGGTGCCAGGGCGGCCTAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGGATGGTCCATTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCgg                                             <  2:76912/150‑1 (MQ=255)

                                                                                                       cGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAg                                                                >  1:113521/1‑131 (MQ=255)

                                                                                                       cGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAg                                                                <  2:113521/131‑1 (MQ=255)

                                                                                                          ccGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGgc                                          >  1:2364800/1‑150 (MQ=255)

                                                                                                          ccGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGgc                                          >  1:1637662/1‑150 (MQ=255)

                                                                                                            ggTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACGGCGCTCGAAGGTTCGTCCATCGCTGGCGGGgcgc                                        <  1:1280257/150‑1 (MQ=255)

                                                                                                              tGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGa                                      <  1:121453/150‑1 (MQ=255)

                                                                                                               gCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGAc                                     <  1:2333520/150‑1 (MQ=255)

                                                                                                                 cAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTg                                   >  1:2142557/1‑150 (MQ=255)

                                                                                                                  aGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGc                                  >  2:776265/1‑150 (MQ=255)

                                                                                                                         cccAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCa                                                                                        <  1:1624947/89‑1 (MQ=255)

                                                                                                                         cccAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCa                                                                                        >  2:1624947/1‑89 (MQ=255)

                                                                                                                          ccAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCGATCGCTGTCGGGGCGCGACTGCCCGCCAAt                          >  1:1924584/1‑150 (MQ=255)

                                                                                                                           cAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATc                         <  1:155174/150‑1 (MQ=255)

                                                                                                                            aGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGacac                                                                                                  <  2:1498946/76‑1 (MQ=255)

                                                                                                                            aGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGacac                                                                                                  >  1:1498946/1‑76 (MQ=255)

                                                                                                                                 gg‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCaaa                   <  2:1637662/150‑1 (MQ=255)

                                                                                                                                    atctagccgaGCTGGCGGCAGCCACCGCGGCATCGGCATCGGCCAGTGTAGAGGACGACGTACGTCCGGTTGAAGT‑GCCAGTCACCGTGTTACGGCGCTCGAAGATACGACAAACGCTGGTTGGGCGCGACTGCCCGCCAATCAGCAAAg                  <  1:2351547/145‑1 (MQ=255)

                                                                                                                                     aaaCGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAgg                 >  2:1328247/1‑150 (MQ=255)

                                                                                                                                         gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCaaa                   <  2:1156266/144‑1 (MQ=255)

                                                                                                                                         gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCaaa                   >  1:1156266/1‑144 (MQ=255)

                                                                                                                                         gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAggagat             >  2:2351547/1‑146 (MQ=255)

                                                                                                                                         gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCAcc             >  2:2216243/1‑150 (MQ=255)

                                                                                                                                         gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGCCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCAcc             >  2:439858/1‑150 (MQ=255)

                                                                                                                                          ccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGGGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCt            >  2:825429/1‑150 (MQ=255)

                                                                                                                                          ccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCt            <  2:1367742/150‑1 (MQ=255)

                                                                                                                                           cgcgGTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCgg                                             >  2:1689592/1‑116 (MQ=255)

                                                                                                                                           cgcgGTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCgg                                             <  1:1689592/116‑1 (MQ=255)

                                                                                                                                            gcggTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGacac                                                                                                  >  2:2308188/1‑62 (MQ=25)

                                                                                                                                            gcggTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGc          >  1:236412/1‑150 (MQ=255)

                                                                                                                                            gcgcTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGc          <  2:2364800/150‑1 (MQ=255)

                                                                                                                                               cTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAAc       <  2:656329/150‑1 (MQ=255)

                                                                                                                                                  gcagcaGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACAACTTCGCCAGTCAAGGAGGTACAGCGCGCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATg    >  1:900368/1‑150 (MQ=255)

CCATACCAGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATGCT  >  NC_002947/3797744‑3798039

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