BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,186,152	A→T	6.2%	V263E (GTG→GAG)	PP_4560 ←	ribonuclease BN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,186,152	0	A	T	6.2%	938.6 / 34.8	274	V263E (GTG→GAG)	PP_4560	ribonuclease BN
Reads supporting (aligned to +/- strand): ref base A (126/131); new base T (7/10); total (133/141)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 293 total aligned reads displayed.

TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCGGG  >  NC_002947/5186003‑5186301

     ggCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCa                                                                                                                                                  <  2:4641840/150‑1 (MQ=255)

      gCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAg                                                                                                                                                 >  1:2844188/1‑150 (MQ=255)

       cTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGc                                                                                                                                                <  1:2897236/150‑1 (MQ=255)

        ttGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGcc                                                                                                                                               >  1:279175/1‑150 (MQ=255)

            tCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                                           <  2:2060882/150‑1 (MQ=255)

            tCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                                           >  2:3636365/1‑150 (MQ=255)

            tCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGaaa                                                                                                                                           >  1:5886118/1‑150 (MQ=255)

               ggTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGAAatga                                                                                                                                        <  1:4878134/150‑3 (MQ=38)

               ggTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAatat                                                                                                                                        <  1:4839680/150‑1 (MQ=255)

                gTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGCTGCCGTGCTGCTGCTCGGC‑GCCGAAatgaa                                                                                                                                       <  2:6089476/150‑4 (MQ=38)

                gTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAAtata                                                                                                                                       <  2:6030632/150‑1 (MQ=255)

                   tGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGaa                                                                                                                                    <  2:7419899/150‑1 (MQ=255)

                       tGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATac                                                                                                                                <  1:4520856/150‑1 (MQ=255)

                       tGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATac                                                                                                                                >  2:3499430/1‑150 (MQ=255)

                       tGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATac                                                                                                                                >  2:3720171/1‑150 (MQ=255)

                        gCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATaca                                                                                                                               <  1:2311870/150‑1 (MQ=255)

                          tCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagc                                                                                                                             >  1:3277037/1‑150 (MQ=255)

                            gTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaa                                                                                                                           >  1:499475/1‑150 (MQ=255)

                            gTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaa                                                                                                                           >  1:5316630/1‑150 (MQ=255)

                             tGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaac                                                                                                                          >  2:4354361/1‑150 (MQ=255)

                                cccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTACGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacagc                                                                                                                       >  1:1067276/1‑150 (MQ=255)

                                 ccAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATacagcaacagca                                                                                                                      >  2:1587255/1‑150 (MQ=255)

                                    gTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACga                                                                                                                   <  1:4624155/150‑1 (MQ=255)

                                     tACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACgat                                                                                                                  >  2:6344502/1‑150 (MQ=255)

                                       cAGTCAAGGTGTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGGTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACgatga                                                                                                                <  2:196911/150‑1 (MQ=255)

                                        aGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACgatgat                                                                                                               >  1:2972715/1‑150 (MQ=255)

                                          tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTg                                                                                                             >  2:2856939/1‑150 (MQ=255)

                                          tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTg                                                                                                             >  2:3915662/1‑150 (MQ=255)

                                          tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTg                                                                                                             >  2:5728955/1‑150 (MQ=255)

                                          tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTg                                                                                                             >  1:6531409/1‑150 (MQ=255)

                                          tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTg                                                                                                             <  1:1608273/150‑1 (MQ=255)

                                          tCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTg                                                                                                             >  2:3146743/1‑150 (MQ=255)

                                            aaGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCa                                                                                                           >  2:2295255/1‑150 (MQ=255)

                                            aaGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCa                                                                                                           >  2:2372284/1‑150 (MQ=255)

                                             aGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                                          <  1:6182347/150‑1 (MQ=255)

                                             aGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAc                                                                                                          <  2:4389553/150‑1 (MQ=255)

                                              ggTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCAcc                                                                                                         <  1:103688/150‑1 (MQ=255)

                                                tttCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCaa                                                                                                       <  1:1569281/150‑1 (MQ=255)

                                                 ttCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAt                                                                                                      <  1:3489450/150‑1 (MQ=255)

                                                 ttCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAt                                                                                                      >  1:4035988/1‑150 (MQ=255)

                                                  tCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAtg                                                                                                     <  1:1464413/150‑1 (MQ=255)

                                                   ccTTGGGCTGGTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAAggc                                                                                                    <  2:4062794/150‑1 (MQ=255)

                                                       gggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGcc                                                                                                <  2:1405080/150‑1 (MQ=255)

                                                        ggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCg                                                                                               >  2:5862324/1‑150 (MQ=255)

                                                        ggCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCg                                                                                               >  1:6627700/1‑150 (MQ=255)

                                                          cTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTa                                                                                             <  2:6515570/150‑1 (MQ=255)

                                                           ttttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTAc                                                                                            <  2:2101548/150‑1 (MQ=255)

                                                           ttttCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTAc                                                                                            >  2:844588/1‑150 (MQ=255)

                                                              tCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATg                                                                                         <  1:5693439/150‑1 (MQ=255)

                                                               ccccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATgg                                                                                        >  1:4406316/1‑150 (MQ=255)

                                                               ccccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATgg                                                                                        >  2:1793992/1‑150 (MQ=255)

                                                                 ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCg                                                                                      <  1:1326287/150‑1 (MQ=255)

                                                                 ccGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCg                                                                                      <  2:4755430/150‑1 (MQ=255)

                                                                  cGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGt                                                                                     <  1:5948526/150‑1 (MQ=255)

                                                                   gTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGtt                                                                                    <  1:2114614/150‑1 (MQ=255)

                                                                   gTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGtt                                                                                    <  2:2581236/150‑1 (MQ=255)

                                                                      gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTa                                                                                 <  1:3915662/150‑1 (MQ=255)

                                                                      gAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTa                                                                                 >  2:3850820/1‑150 (MQ=255)

                                                                          tCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCg                                                                             <  1:3607192/150‑1 (MQ=255)

                                                                            tttCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGc                                                                           >  1:3832103/1‑150 (MQ=255)

                                                                            tttCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGc                                                                           <  1:4238646/150‑1 (MQ=255)

                                                                                tcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                                       <  2:4951493/150‑1 (MQ=255)

                                                                                tcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                                       >  2:3269696/1‑150 (MQ=255)

                                                                                tcGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAg                                                                       >  1:5115414/1‑150 (MQ=255)

                                                                                  gCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGt                                                                     >  2:6233182/1‑150 (MQ=255)

                                                                                   cccGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTc                                                                    <  2:5175579/150‑1 (MQ=255)

                                                                                     cGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCtt                                                                  <  1:572766/150‑1 (MQ=255)

                                                                                           cTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgta                                                            <  1:5434821/150‑1 (MQ=255)

                                                                                           cTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgta                                                            <  1:4225361/150‑1 (MQ=255)

                                                                                            ttGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtag                                                           >  2:6785676/1‑150 (MQ=255)

                                                                                              gCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagta                                                         >  2:7080949/1‑150 (MQ=255)

                                                                                               ccTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGtagtag                                                        >  2:789710/1‑150 (MQ=255)

                                                                                                cTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAgg                                                       >  1:2326826/1‑150 (MQ=255)

                                                                                                 ttCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGc                                                      >  1:8229/1‑150 (MQ=255)

                                                                                                       tGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGcc                                                >  2:3372033/1‑150 (MQ=255)

                                                                                                       tGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGcc                                                >  1:1720286/1‑150 (MQ=255)

                                                                                                        gACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGccc                                               <  1:4133801/150‑1 (MQ=255)

                                                                                                          cATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAg                                             <  2:5359331/150‑1 (MQ=255)

                                                                                                          cATGTGCTCGATCACCGCATTCATTCCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAg                                             <  2:1067276/150‑1 (MQ=255)

                                                                                                           aTGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAgg                                            <  1:6785676/150‑1 (MQ=255)

                                                                                                             gtgCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGa                                          <  2:7361599/150‑1 (MQ=255)

                                                                                                              tgCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAg                                         >  2:717297/1‑150 (MQ=255)

                                                                                                              tgCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAg                                         <  1:6016724/150‑1 (MQ=255)

                                                                                                                cTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGTACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGc                                       >  2:6878497/1‑150 (MQ=255)

                                                                                                                 tCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagta                                                         >  1:2328718/1‑131 (MQ=255)

                                                                                                                 tCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGcc                                      >  2:5460303/1‑150 (MQ=255)

                                                                                                                        ccGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATcc                               >  1:2456682/1‑150 (MQ=255)

                                                                                                                        ccGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATcc                               >  1:761001/1‑150 (MQ=255)

                                                                                                                         cGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagta                                                         >  1:2620470/1‑123 (MQ=255)

                                                                                                                         cGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACgtagta                                                         <  2:2620470/123‑1 (MQ=255)

                                                                                                                           cattcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccaca                            <  1:6237389/150‑1 (MQ=255)

                                                                                                                           cattcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccaca                            <  2:4406316/150‑1 (MQ=255)

                                                                                                                            attcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacac                           <  1:5599838/150‑1 (MQ=255)

                                                                                                                            attcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacac                           <  1:1897900/150‑1 (MQ=255)

                                                                                                                             ttcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacacc                          >  1:727559/1‑150 (MQ=255)

                                                                                                                              tcattTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATccacacca                         >  1:7002758/1‑150 (MQ=255)

                                                                                                                                      gctgcCGTGCTGCTGCTCGGC‑GCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGcc                 >  2:4276721/4‑150 (MQ=255)

                                                                                                                                              gcagcagcaCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCg          <  2:5886118/150‑1 (MQ=255)

                                                                                                                                                     caCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCgg   >  2:5614801/1‑150 (MQ=255)

                                                                                                                                                     caCGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCgg   <  1:6435150/150‑1 (MQ=255)

TATGAGGCTTGCTCGGGTTTGGGTGCTCGTGACCCAGTACAGTCAAGGTTTCCTTGGGCTTTTCCCCGTCGAAGTCTTTCTCGCCCGGGTTCTTGCCTTCGGCTGACATGTGCTCGATCACCGCATTCATTTCGGC‑GCCGAGCAGCAGCACGGCAGCCGAAATATAGAAATACAGCAACAGCACGATGATTGCACCAATGCTGCCGTACATGGCGTTGTAGTCGGCAAAGGTCTTCACGTAGTAGGCAAAGCCCAGGGAGGCCACGATCCACACCACAACCGCCAGCACCGAGCCCGGG  >  NC_002947/5186003‑5186301

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