BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,990,454	C→G	P71P (CCG→CCC)	ycaC‑III ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,990,454	0	C	G	97.2%	3061.6 / 66.0	1584	P71P (CCG→CCC)	ycaC‑III	putative hydrolase
Reads supporting (aligned to +/- strand): ref base C (21/23); new base G (779/756); total (801/783)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 1638 total aligned reads displayed.

CACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGCCA  >  NC_002947/5990307‑5990598

 aCGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCggg                                                                                                                                                >  1:75229/1‑149 (MQ=18)

  cGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCggg                                                                                                                                                >  6:244882/1‑148 (MQ=18)

  cGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCggg                                                                                                                                                >  1:378525/1‑148 (MQ=18)

  cGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGc                                                                                                                                               <  4:141497/149‑1 (MQ=18)

   gCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCa                                                                                                                                              >  3:44275/1‑149 (MQ=18)

   gCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑CCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCACGGCGGGCAATGTACGCCGCATCCGGGAACATCACTTACAGCTCGGGCa                                                                                                                                              >  3:113632/1‑149 (MQ=18)

   gCAAACCTCGGTCACAACACCTGCAATGAACAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCa                                                                                                                                              >  1:230258/1‑149 (MQ=18)

    cAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCAc                                                                                                                                             <  8:273097/149‑1 (MQ=18)

      aaCCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCa                                                                                                                                           <  4:341674/149‑1 (MQ=18)

       aCCGCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAg                                                                                                                                          <  8:4697/149‑1 (MQ=17)

        ccTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGc                                                                                                                                         >  7:337507/1‑149 (MQ=18)

         cTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCg                                                                                                                                        >  5:189023/1‑149 (MQ=17)

          tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                                       >  3:315036/1‑149 (MQ=17)

           cGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCggg                                                                                                                                      >  2:344173/1‑149 (MQ=17)

            ggTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCCGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                                       <  6:216600/147‑1 (MQ=17)

            ggTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                                       <  2:65795/147‑1 (MQ=17)

            ggTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                                       <  2:36337/147‑1 (MQ=17)

             gTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGcc                                                                                                                                    <  5:141272/149‑1 (MQ=17)

              tCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGcca                                                                                                                                   <  3:372144/149‑2 (MQ=17)

              tCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGACTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGcca                                                                                                                                   >  1:91458/1‑148 (MQ=17)

               cacaACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGccat                                                                                                                                  >  1:187583/1‑147 (MQ=17)

                acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCATTCACCGCATTGACGAAATACTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATACGGGAACATCTCTTTAAGCTCGGGCACCAGCGGGCCAtt                                                                                                                                 >  8:207643/1‑149 (MQ=17)

                 caacaCGTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAttc                                                                                                                                <  5:397998/149‑2 (MQ=17)

                  aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAttcg                                                                                                                               >  2:116712/1‑147 (MQ=17)

                  aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAttcg                                                                                                                               >  4:307471/1‑147 (MQ=17)

                  aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAttcg                                                                                                                               <  5:215073/149‑3 (MQ=17)

                   acacCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTcgg                                                                                                                              >  2:150467/1‑149 (MQ=17)

                   acacCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTcgg                                                                                                                              <  5:262542/149‑1 (MQ=17)

                   acacCTGCAATGATCAGCTGCTTCTTGCCGGGC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTcgg                                                                                                                              <  4:214750/149‑1 (MQ=17)

                    cacCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTcggg                                                                                                                             <  2:202434/149‑1 (MQ=17)

                    cacCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTcggg                                                                                                                             <  6:68818/149‑1 (MQ=17)

                     acCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGc                                                                                                                            <  6:371820/149‑1 (MQ=17)

                      ccTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGc                                                                                                                                         <  6:109957/135‑1 (MQ=17)

                      ccTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGcc                                                                                                                           <  4:102039/149‑1 (MQ=17)

                        tGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCtt                                                                                                                         >  7:39029/1‑149 (MQ=14)

                         gCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAGTGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCAttc                                                                                                                                <  1:169969/141‑2 (MQ=17)

                         gCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTGGGGCACCAGCGGGCCATTCGGGc                                                                                                                            <  6:35341/145‑1 (MQ=14)

                         gCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGc                                                                                                                                     <  8:246962/136‑1 (MQ=17)

                         gCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTg                                                                                                                        >  5:102213/1‑149 (MQ=14)

                          cAATGATCAGCTGCTTCTTGCCGGTCGGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTg                                                                                                                        >  1:5141/1‑149 (MQ=14)

                          cAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCtt                                                                                                                         <  8:360243/147‑1 (MQ=14)

                          cAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTg                                                                                                                        <  5:165646/148‑1 (MQ=14)

                          cAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGc                                                                                                                       >  2:48229/1‑149 (MQ=14)

                          cAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGc                                                                                                                       <  8:322128/149‑1 (MQ=14)

                          cAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGc                                                                                                                       >  7:159208/1‑149 (MQ=14)

                          cAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGc                                                                                                                       <  4:379717/149‑1 (MQ=14)

                           aaTGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCt                                                                                                                      >  3:141808/1‑149 (MQ=14)

                           aaTGATCAGCTGCTGCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCt                                                                                                                      <  5:326630/149‑1 (MQ=14)

                            aTGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGTTTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                                                                                                     <  1:102248/149‑1 (MQ=14)

                               aTCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCg                                                                                                                    >  1:330213/1‑147 (MQ=14)

                               aTCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                                  <  8:176920/149‑1 (MQ=14)

                               aTCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                                  <  4:49588/149‑1 (MQ=14)

                                tCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCTGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAg                                                                                                                 <  6:189023/149‑1 (MQ=14)

                                tCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                                       >  1:30927/1‑127 (MQ=17)

                                tCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                                  <  4:166161/148‑1 (MQ=14)

                                tCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                                  <  2:92496/148‑1 (MQ=14)

                                tCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAg                                                                                                                 >  3:247209/1‑149 (MQ=14)

                                tCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAg                                                                                                                 >  6:41875/1‑149 (MQ=14)

                                  aGCTGCTTCTTGCCGGTC‑GCCTTCGCCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                               >  5:295668/1‑149 (MQ=14)

                                  aGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                               >  7:232609/1‑149 (MQ=14)

                                   gctgctTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTTGCCAGGGCGGGCAATGTACGGCGCTTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtg                                                                                                              <  2:245323/149‑1 (MQ=14)

                                   gctgctTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtg                                                                                                              <  1:274036/149‑1 (MQ=14)

                                    ctgctTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGc                                                                                                                            >  1:356230/1‑134 (MQ=14)

                                    ctgctTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtgg                                                                                                             <  7:70431/149‑1 (MQ=14)

                                     tgctTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggt                                                                                                            >  7:31209/1‑149 (MQ=14)

                                      gctTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                                                                                           >  8:68577/1‑149 (MQ=14)

                                      gaaatatGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                                                                                           <  4:248233/143‑1 (MQ=11)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGATGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtgg                                                                                                          <  5:348182/149‑1 (MQ=14)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                                       <  7:101740/120‑1 (MQ=14)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGc                                                                                                                       >  6:31770/1‑136 (MQ=14)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCt                                                                                                                      >  5:37129/1‑137 (MQ=14)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCt                                                                                                                      >  3:204857/1‑137 (MQ=14)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtgg                                                                                                          >  7:179144/1‑149 (MQ=14)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                                                                                           <  2:316048/148‑1 (MQ=14)

                                       cttcttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggt                                                                                                            <  5:53512/147‑1 (MQ=14)

                                          cttGCCGGTC‑GCCTTCATCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGcc                                                                                                                           >  4:847/1‑129 (MQ=11)

                                          cttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTcgg                                                                                                                              >  8:148687/1‑126 (MQ=11)

                                           ttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCGCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCggg                                                                                                                                      >  8:302120/1‑117 (MQ=14)

                                           ttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                                       >  3:79905/1‑116 (MQ=14)

                                           ttGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCt                                                                                                                      >  1:141186/1‑133 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAATCTGGTGGTGAgg                                                                                                     <  2:273520/148‑1 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                               >  3:340702/1‑138 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                       <  3:281696/146‑1 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                       <  6:135890/146‑1 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                                       <  7:208868/146‑1 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAgg                                                                                                     <  2:179943/148‑1 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                                    >  7:132579/1‑149 (MQ=11)

                                             gCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                                    >  6:166740/1‑149 (MQ=11)

                                              ccGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                                    <  1:208707/148‑1 (MQ=11)

                                              ccGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                                                                   <  7:260713/149‑1 (MQ=11)

                                              ccGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                                                                   >  4:272961/1‑149 (MQ=11)

                                              ccGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                                                                   <  5:243287/149‑1 (MQ=11)

                                              ccGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGAt                                                                                                   >  6:282865/1‑149 (MQ=11)

                                              ccGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGGGAGGAt                                                                                                   >  4:124627/1‑149 (MQ=11)

                                               cGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATg                                                                                                  <  8:117542/149‑1 (MQ=11)

                                               cGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATg                                                                                                  >  3:169508/1‑149 (MQ=11)

                                                 gTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                                <  5:76522/149‑1 (MQ=11)

                                                      ccTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGcc                                                                                                                                    >  2:296734/1‑109 (MQ=11)

                                                        ttCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGccat                                                                                                                                  <  7:293522/109‑3 (MQ=11)

                                                                                                                           cATCCGTGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATcg                       <  7:97320/149‑1 (MQ=255)

CACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTC‑GCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCAGGTCAGCCAGCGCCAGCACGTTGTTCTTGAACTTGTCTGGCTCGATATCGCGCACCAGCGACAGCAGGCCA  >  NC_002947/5990307‑5990598

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