BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	160,311	A→G	5.1%	L26L (CTT→CTC)	PP_0149 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	160,311	0	A	G	5.1%	296.8 / 3.9	98	L26L (CTT→CTC)	PP_0149	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (48/45); new base G (2/3); total (50/48)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 102 total aligned reads displayed.

CGTCGGCAGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGG  >  NC_002947/160163‑160457

cGTCGGCAGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCa                                                                                                                                                     <  5:249408/149‑1 (MQ=255)

  tCGGCAGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAg                                                                                                                                                   <  3:199986/149‑1 (MQ=255)

        ggcaTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTTGAGCTGCAg                                                                                                                                             <  2:265208/149‑1 (MQ=255)

         gcaTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAgg                                                                                                                                            >  7:168724/1‑149 (MQ=255)

          caTGTTGTCGGTCCAGCTGAACCGCACGCGGCGGTGCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTTGAGCTGCAggc                                                                                                                                           <  3:122036/149‑2 (MQ=255)

          caTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTTGAGCTGCAggc                                                                                                                                           <  6:199092/149‑2 (MQ=255)

           aTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTc                                                                                                                                          <  7:312235/149‑1 (MQ=255)

           aTGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTc                                                                                                                                          >  6:344142/1‑149 (MQ=255)

               tgtCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑c                                                                                                                                      >  1:332908/1‑148 (MQ=255)

                gtCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCg                                                                                                                                         <  4:216583/145‑1 (MQ=255)

                gtCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCg                                                                                                                                         >  3:216583/1‑145 (MQ=255)

                      ccAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAggcgg                                                                                                                              <  6:24089/149‑1 (MQ=255)

                       cAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGGGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAggcgg                                                                                                                              <  8:233748/148‑1 (MQ=255)

                       cAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAggcggg                                                                                                                             >  8:226994/1‑149 (MQ=255)

                        aGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAggcgggc                                                                                                                            <  4:240654/149‑1 (MQ=255)

                             aCCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGcc                                                                                                                       >  6:91002/1‑149 (MQ=255)

                             aCCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGcc                                                                                                                       >  7:346216/1‑149 (MQ=255)

                               ccGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGcc                                                                                                                       >  7:11042/1‑147 (MQ=255)

                               ccGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGa                                                                                                                     <  5:91002/149‑1 (MQ=255)

                                cGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGAGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGag                                                                                                                    >  8:279739/1‑149 (MQ=255)

                                 gCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGagc                                                                                                                   <  3:186719/149‑1 (MQ=255)

                                 gCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGag                                                                                                                    >  7:10803/1‑148 (MQ=255)

                                        gcggTCCAGCTGTTGCTTGAAGCGCGGTGGAAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCa                                                                                                            <  2:120805/149‑1 (MQ=255)

                                         cggTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGGCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAg                                                                                                           >  6:242067/1‑149 (MQ=255)

                                          ggTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAgg                                                                                                          <  4:7570/149‑1 (MQ=255)

                                            tCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGgc                                                                                                        >  1:179310/1‑149 (MQ=255)

                                              cAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGgcgc                                                                                                      >  8:208190/1‑149 (MQ=255)

                                               aGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGgcg                                                                                                       >  7:24693/1‑147 (MQ=255)

                                               aGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGgcg                                                                                                       >  6:165568/1‑147 (MQ=255)

                                                gCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAg                                                                                                    <  2:332908/149‑1 (MQ=255)

                                                 cTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGGTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGc                                                                                                   >  3:201578/1‑149 (MQ=255)

                                                  tGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGcc                                                                                                  <  5:95833/149‑1 (MQ=255)

                                                    ttGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGGCCCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGcc                                                                                                  <  1:258286/148‑1 (MQ=255)

                                                    ttGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGcc                                                                                                  >  7:193622/1‑147 (MQ=255)

                                                    ttGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGcc                                                                                                  >  6:85788/1‑147 (MQ=255)

                                                    tggCTTGACGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGGTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCaa                                                                                                <  6:62085/147‑1 (MQ=255)

                                                      gCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAag                                                                                               >  4:232624/1‑148 (MQ=255)

                                                      gCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAag                                                                                               >  7:205294/1‑148 (MQ=255)

                                                       cTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAagc                                                                                              <  8:245734/148‑1 (MQ=255)

                                                        ttGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCa                                                                                                                                              >  6:343828/1‑100 (MQ=255)

                                                        ttGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCa                                                                                                                                              <  5:343828/100‑1 (MQ=255)

                                                          gAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAagcagc                                                                                           >  7:321146/1‑148 (MQ=255)

                                                           aaGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAagcagc                                                                                           >  3:360776/1‑147 (MQ=255)

                                                           aaGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAagcagc                                                                                           >  7:1975/1‑147 (MQ=255)

                                                            aGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAg                                                                                        <  7:338070/149‑1 (MQ=255)

                                                            aGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAg                                                                                        <  7:342939/149‑1 (MQ=255)

                                                             gcgcGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCGGCAGGTCGG‑CCAAGGCGGGCGTGCCGCGCAGCGTCATGGCGCATCCAAGCAGCCAgg                                                                                       >  7:309463/1‑149 (MQ=255)

                                                                cgGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGgcac                                                                                    <  5:324693/149‑1 (MQ=255)

                                                                 ggTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGgca                                                                                     <  7:279739/147‑1 (MQ=255)

                                                                    ggCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGgcacgcac                                                                                <  8:24693/149‑1 (MQ=255)

                                                                    ggCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGGCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGgcacgcac                                                                                <  3:41414/149‑1 (MQ=255)

                                                                        gCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACgcgc                                                                            <  8:168724/149‑1 (MQ=255)

                                                                            gCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCtt                                                                          <  3:134496/147‑1 (MQ=255)

                                                                             cccATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTca                                                                        >  2:369672/1‑148 (MQ=255)

                                                                              ccATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTcaca                                                                      <  7:208190/149‑1 (MQ=255)

                                                                              ccATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTca                                                                        >  1:15987/1‑147 (MQ=255)

                                                                               cATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAg                                                                     >  3:189348/1‑149 (MQ=255)

                                                                                  ggCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTgg                                                                  >  4:359114/1‑148 (MQ=255)

                                                                                     cTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGACAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAg                                                               >  6:191943/1‑149 (MQ=255)

                                                                                          cgAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATgg                                                          >  8:324605/1‑149 (MQ=255)

                                                                                           gAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGc                                                         <  7:226994/149‑1 (MQ=255)

                                                                                           gAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGGCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATgg                                                          <  2:134484/148‑1 (MQ=255)

                                                                                                  gCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGc                                                   >  4:239893/1‑148 (MQ=255)

                                                                                                   cTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTc                                                      <  8:175530/144‑1 (MQ=255)

                                                                                                   cTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTc                                                      >  7:175530/1‑144 (MQ=255)

                                                                                                        gAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATAcc                                            <  4:360776/149‑1 (MQ=255)

                                                                                                           gCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCaa                                                >  6:81412/1‑142 (MQ=255)

                                                                                                           gCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCaa                                                <  5:81412/142‑1 (MQ=255)

                                                                                                             tttAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTTCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGcacgcacg                                                                               >  7:173414/1‑109 (MQ=255)

                                                                                                             tttAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGcacgcacg                                                                               <  8:173414/109‑1 (MQ=255)

                                                                                                             tttAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTc                                       >  8:316068/1‑149 (MQ=255)

                                                                                                                 gTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCt                                    <  5:242067/148‑1 (MQ=255)

                                                                                                                  tttGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTgg                                  <  3:232624/149‑1 (MQ=255)

                                                                                                                      ttgtGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGcgcc                              >  5:29710/1‑149 (MQ=255)

                                                                                                                        gtgGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGcgccgc                            >  4:91136/1‑149 (MQ=255)

                                                                                                                          gggTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCGGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCcgcgc                          >  1:250065/1‑149 (MQ=255)

                                                                                                                           ggTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCcgcgc                          <  4:189112/148‑1 (MQ=255)

                                                                                                                            gTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCgcgcgc                        >  3:45078/1‑149 (MQ=255)

                                                                                                                                 ggCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCg                   <  8:1975/149‑1 (MQ=255)

                                                                                                                                  gCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCgg                  >  1:25743/1‑149 (MQ=255)

                                                                                                                                  gCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCgg                  <  7:324605/149‑1 (MQ=255)

                                                                                                                                   ccGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCTGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCgg                  >  3:222455/1‑148 (MQ=255)

                                                                                                                                    cGGCCGCCTGCAGCTCAAGCTGCAGGGCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTa                <  2:269555/149‑1 (MQ=255)

                                                                                                                                     ggACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTAc               >  5:363973/1‑149 (MQ=255)

                                                                                                                                       cgaccTGCAGCTTGAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCCAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActc             >  1:378245/4‑149 (MQ=255)

                                                                                                                                       aCGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActc             <  8:11042/149‑1 (MQ=255)

                                                                                                                                           cTGCAGCTTGAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc         >  2:228232/1‑149 (MQ=255)

                                                                                                                                           cTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCaa                                                                                                <  8:215626/62‑1 (MQ=255)

                                                                                                                                           cTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCaa                                                                                                >  7:215626/1‑62 (MQ=255)

                                                                                                                                           cTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc         <  4:201578/149‑1 (MQ=255)

                                                                                                                                           cTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc         >  5:185754/1‑149 (MQ=255)

                                                                                                                                           cTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctc         <  8:321146/149‑1 (MQ=255)

                                                                                                                                           cTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgc           >  6:284988/1‑147 (MQ=255)

                                                                                                                                           cTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgc           <  5:284988/147‑1 (MQ=255)

                                                                                                                                            tGCAGCTCAAGCTGCAGGCCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTActcgctcg        >  3:224505/1‑149 (MQ=255)

                                                                                                                                                 cTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAg   <  6:29710/149‑1 (MQ=255)

                                                                                                                                                 cTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAg   >  3:308582/1‑149 (MQ=255)

                                                                                                                                                    aaGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCTGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAgg  <  4:222455/147‑1 (MQ=255)

                                                                                                                                                    aaGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTg                                   <  2:53586/114‑1 (MQ=255)

                                                                                                                                                    aaGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTg                                   >  1:53586/1‑114 (MQ=255)

CGTCGGCAGGCATGTTGTCGCTCCAGCTGACCCGCACGCGGCGGTCCAGCTGTTGCTTGAAGCGCGGTGGCAGCTTGCCCATGGCCTCGTCGAGCAGGGCTTGGGAGGCTTTAGTTTGTTGTGGGTCGAGGCCGGACGCCTGCAGCTCAAGCTGCAGGTCGG‑CCAAGGCGGGCGTGCCGAGCAGCGTCAGGGCGCAGCCAAGCAGCCAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGG  >  NC_002947/160163‑160457

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: `‑`