BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,596	C→G	5.4%	Q66H (CAG→CAC)	PP_1399 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,598,596	0	C	G	5.4%	693.6 / 3.0	229	Q66H (CAG→CAC)	PP_1399	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (115/93); new base G (5/11); total (122/107)
Fisher's exact test for biased strand distribution p-value = 7.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

GGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGAAGGCGC  >  NC_002947/1598449‑1598744

  tGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGgc                                                                                                                                                  <  1:1163438/150‑1 (MQ=34)

     ggTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGgccgc                                                                                                                                               >  2:1279759/1‑150 (MQ=35)

      gTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGTTGGCCACCGGTGGGGTCGTAGCCGTCTTCGGCGATGGGGTGgcctgg                                                                                                                                              >  2:3903098/1‑147 (MQ=11)

       tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGccgccg                                                                                                                                             <  2:6987854/150‑1 (MQ=37)

       tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCACTGGccgccg                                                                                                                                             <  1:2420807/150‑1 (MQ=37)

        cGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCgg                                                                                                                                            <  2:501144/150‑1 (MQ=37)

         gAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCGGGCCGCCGGt                                                                                                                                           >  2:1196748/1‑150 (MQ=37)

              cGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGt                                                                                                                                      <  1:3059019/150‑1 (MQ=37)

                   tCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAg                                                                                                                                 >  2:29797/1‑150 (MQ=37)

                     aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCCCCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGcc                                                                                                                               <  2:2554127/150‑1 (MQ=37)

                     aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGcc                                                                                                                               <  2:3287084/150‑1 (MQ=37)

                      gCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGc                                                                                                                                <  1:7036271/148‑1 (MQ=37)

                              ttGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGc                                                                                                                      >  2:2769778/1‑150 (MQ=37)

                                ggCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGa                                                                                                                    >  2:5983970/1‑150 (MQ=37)

                                 gCCACCGGTGCTGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGaa                                                                                                                   <  1:893506/150‑1 (MQ=37)

                                  ccACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGACGGTGCGGTCGAAGCCGTCTTCGGCGAAg                                                                                                                  >  1:3381133/1‑150 (MQ=37)

                                      cGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGtt                                                                                                              <  1:6462771/150‑1 (MQ=38)

                                      cGGGGCGGCCGACGCCCTCTGCAGCGAAGCGTTGCCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGtt                                                                                                              <  1:5974017/150‑1 (MQ=11)

                                       ggTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTg                                                                                                             <  2:6310459/150‑1 (MQ=39)

                                       ggTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTg                                                                                                             <  2:113396/150‑1 (MQ=39)

                                       ggTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTg                                                                                                             <  2:5041081/150‑1 (MQ=39)

                                       ggTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTg                                                                                                             <  2:472217/150‑1 (MQ=39)

                                         tGCGGTCTAACCCGTATTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGATGGCCTCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGc                                                                                                           <  2:2483666/150‑1 (MQ=11)

                                            ggTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGTAGCCGTCTTCGGCGAAGCGGTGGcccg                                                                                                        >  2:1261291/1‑148 (MQ=255)

                                            ggTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCAc                                                                                                        <  2:5954184/150‑1 (MQ=39)

                                            ggTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCAc                                                                                                        >  2:4876834/1‑150 (MQ=38)

                                             gTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGATGGCCAcc                                                                                                       >  1:3803425/1‑150 (MQ=255)

                                                 aaGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTg                                                                                                   >  1:3960910/1‑150 (MQ=39)

                                                  aGCCGTCTGCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGCCCACCGGTGCGGTCGAAGCCGACTGCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGACATCGGCCAAGCGGTGGCCACCGGTGc                                                                                                  >  2:791270/1‑150 (MQ=14)

                                                   gCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCg                                                                                                 <  2:1734199/150‑1 (MQ=37)

                                                   gCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGGTGGCCCCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCg                                                                                                 >  2:2695193/1‑150 (MQ=255)

                                                    ccGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTg                                                                                                             >  1:1245840/1‑137 (MQ=40)

                                                         ttCAGCGAATCTTTGGCCTACGGTGCTGTCTAAGCCGTCTTCAGCGAAGCGTTGGCCACAGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGaa                                                                                           <  2:4151269/150‑1 (MQ=11)

                                                         ttCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGGTGGCCGACGGTGCGGTCGAAGCCGTCTTCGGCTAAGCGTTTGCCACCGGTGAGGTctta                                                                                           >  1:178669/1‑147 (MQ=11)

                                                          tCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAg                                                                                          >  2:1970206/1‑150 (MQ=39)

                                                            aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGcc                                                                                        >  1:1828785/1‑150 (MQ=39)

                                                            aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGcc                                                                                        >  1:614016/1‑150 (MQ=39)

                                                              cGACGCGTTGGCCCCCGGTGCGGTCGAAGCCGTATTCAGCGAAGCGTTGGCGACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGt                                                                                      <  2:6423856/150‑1 (MQ=255)

                                                              cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGt                                                                                      >  1:2825058/1‑150 (MQ=39)

                                                              cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGt                                                                                      >  1:7242143/1‑150 (MQ=39)

                                                              cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGt                                                                                      >  2:7099939/1‑150 (MQ=39)

                                                               gAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTc                                                                                     >  2:6756512/1‑150 (MQ=39)

                                                                aaGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCt                                                                                    >  1:7183664/1‑150 (MQ=39)

                                                                   cGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCa                                                                                 <  2:5901558/150‑1 (MQ=40)

                                                                        gCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTa                                                                            >  1:1027098/1‑150 (MQ=255)

                                                                         ccACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAg                                                                           >  1:5727898/1‑150 (MQ=255)

                                                                           aCCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCg                                                                         <  2:1599400/150‑1 (MQ=255)

                                                                             cGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCggg                                                                       >  2:2202503/1‑150 (MQ=255)

                                                                               gTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGc                                                                     >  2:3694351/1‑150 (MQ=255)

                                                                                   ggTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTg                                                                                                                                                    <  2:5743992/67‑1 (MQ=40)

                                                                                    gTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTg                                                                >  1:4209316/1‑150 (MQ=255)

                                                                                     tCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCAGCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTgg                                                               >  2:4008204/1‑150 (MQ=255)

                                                                                       gAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGt                                                             >  1:6427329/1‑150 (MQ=255)

                                                                                        aaGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGtt                                                            <  2:5651583/150‑1 (MQ=255)

                                                                                         aGCCTTCTGCCGCGAAGCGTTGGCGCCGGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGttt                                                           <  2:1755309/150‑1 (MQ=17)

                                                                                          gCCGTCTTCAGCGATGCGTTGGCGACGGGTGCGGTCTAAGCCGTCTTAGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTg                                                          <  2:2892105/150‑1 (MQ=17)

                                                                                          cgctccTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCCGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTg                                                          <  1:1461710/145‑1 (MQ=255)

                                                                                             gTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTg                                                                                                                                          <  2:2266431/67‑1 (MQ=255)

                                                                                             gTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGATGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGCCAAACGGTTAGGACTAAGCCGTATTATTCGTAGAGGGGCCACTGGATTTCAAc                                                       >  1:5240314/1‑150 (MQ=14)

                                                                                               cctCAGCGAAGCGGTGGGCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTg                                                     <  2:4507996/148‑1 (MQ=255)

                                                                                                ttCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGc                                                    >  2:6792106/1‑150 (MQ=255)

                                                                                                ttCAGCGAAGCGGTGGCGACCGGTGCTGTCGAAGACGTATAAGGCGAAGCGGTGGGCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGc                                                    <  2:897011/150‑1 (MQ=17)

                                                                                                 tCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGAGGGGCGACTGGGTGAGAACTGCg                                                   >  2:3288776/1‑150 (MQ=255)

                                                                                                  cAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCgg                                                  >  1:676892/1‑150 (MQ=255)

                                                                                                   aGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGc                                                 <  1:3211166/150‑1 (MQ=255)

                                                                                                    gCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCa                                                >  1:248833/1‑150 (MQ=255)

                                                                                                    gCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCa                                                >  1:6088937/1‑150 (MQ=255)

                                                                                                     cGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCaa                                               >  2:5582632/1‑150 (MQ=255)

                                                                                                      gAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGACGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCaaa                                              <  1:3125188/150‑1 (MQ=255)

                                                                                                           gTTGGCCTCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGgtgt                                         <  1:7149056/150‑1 (MQ=255)

                                                                                                            ttGGCCACCGGTGCGGTCGAAGCCGTCTTCTGCGAAGCGCTGGCCGCCGGTGCTGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGtgtg                                        <  1:4893036/150‑1 (MQ=255)

                                                                                                              ggCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAg                                             <  1:3564904/143‑1 (MQ=255)

                                                                                                              ggCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAg                                             >  2:3564904/1‑143 (MQ=255)

                                                                                                                ccACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCg                                    >  2:1324718/1‑150 (MQ=255)

                                                                                                                ccACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCg                                    >  1:5016482/1‑150 (MQ=255)

                                                                                                                ccACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCg                                    >  2:6451662/1‑150 (MQ=255)

                                                                                                                 cgttggTGCGGTCGAAGCCGTCTTCAGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGt                                   <  1:5761140/146‑1 (MQ=255)

                                                                                                                 cACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGt                                   >  1:2537815/1‑150 (MQ=255)

                                                                                                                   ccGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcg                                 <  1:1594831/150‑1 (MQ=255)

                                                                                                                   ccGGCGCGGGCGCAGCCGTCTTCGGCGCAGCGCTGGCCGCCGGAGCGGGCGAAGCCGTCTTCCGCGAAGCGTGGGCCTCCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcg                                 <  2:2343245/150‑1 (MQ=21)

                                                                                                                    cGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCa                                                                                 >  2:5240417/1‑101 (MQ=255)

                                                                                                                       tGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcggtcc                             >  1:2587399/1‑150 (MQ=255)

                                                                                                                         cGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcggtccgg                           <  2:7206556/150‑1 (MQ=255)

                                                                                                                           gTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcggtccggtc                         >  1:6767486/1‑150 (MQ=255)

                                                                                                                           gTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcggtccggtc                         >  1:4774057/1‑150 (MQ=255)

                                                                                                                           gTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTcggtccggtc                         >  1:233832/1‑150 (MQ=255)

                                                                                                                            tCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCg                        >  2:263302/1‑150 (MQ=255)

                                                                                                                            tCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCg                        >  1:1517694/1‑150 (MQ=255)

                                                                                                                              gAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGaa                      >  1:7100869/1‑150 (MQ=255)

                                                                                                                              gAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGaa                      <  1:722121/150‑1 (MQ=255)

                                                                                                                              gAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGaa                      <  2:7242143/150‑1 (MQ=255)

                                                                                                                                 gCCGTCGTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGcc                   <  2:6019654/150‑1 (MQ=255)

                                                                                                                                    gTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATc                >  1:6316399/1‑150 (MQ=255)

                                                                                                                                    gTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTGGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGGCGGTCCGGTCGAAGCCATc                >  1:588347/1‑150 (MQ=255)

                                                                                                                                     tcttctGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCt               <  1:1911324/150‑1 (MQ=255)

                                                                                                                                       ttCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTc             >  1:6529550/1‑150 (MQ=255)

                                                                                                                                       ttCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGCGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTc             >  1:4349331/1‑150 (MQ=255)

                                                                                                                                        tCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCg            <  1:4599030/150‑1 (MQ=255)

                                                                                                                                         cGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCgg           >  2:6199622/1‑150 (MQ=255)

                                                                                                                                         cGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCgg           <  1:1196613/150‑1 (MQ=255)

GGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGCTGGCCGCCGGTGCGGTCGAAGCCGTCTTCGGCGAAGCGTTGGCCACCGGTGCGGTCGAAGCCGTCTTCAGCGTAGCGGGGCGACTGGGTTTGAACTGCGGCAAAGGTGTGAGCGTCGGTCCGGTCGAAGCCATCTTCGGCGAAGGCGC  >  NC_002947/1598449‑1598744

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