BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	944,837	T→C	5.1%	E263G (GAG→GGG)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	944,837	0	T	C	5.1%	766.7 / 17.6	234	E263G (GAG→GGG)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand): ref base T (106/116); new base C (6/6); total (112/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 236 total aligned reads displayed.

TGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCGCTGTTG  >  NC_002947/944688‑944986

tGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCt                                                                                                                                                       >  2:3775099/1‑150 (MQ=255)

  ccAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCg                                                                                                                                                     <  1:504728/150‑1 (MQ=255)

  ccAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCg                                                                                                                                                     <  1:2284934/150‑1 (MQ=255)

    aGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGTGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGTGCCGAGCCCGGc                                                                                                                                                   <  1:4613259/150‑1 (MQ=255)

    aGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGc                                                                                                                                                   <  1:2591879/150‑1 (MQ=255)

     gCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCCCCGCAGTGCCCCGCAATACGCCACCTGCGCCTCGCACGGAACCCGTCGCGCCGGGCTCGGcc                                                                                                                                                  >  1:5390025/1‑149 (MQ=255)

     gCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCa                                                                                                                                                  <  2:2692109/150‑1 (MQ=255)

      ccGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCAc                                                                                                                                                 <  1:992796/150‑1 (MQ=255)

        gcgTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGa                                                                                                                                               <  2:2952085/150‑1 (MQ=255)

         cgTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCACGCACGGAAACCGTCGTGCCGAGCCCGGCTCGAc                                                                                                                                              <  2:5555696/150‑1 (MQ=255)

             tGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGt                                                                                                                                          >  1:3277121/1‑150 (MQ=255)

               cggcgCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTaa                                                                                                                                        <  1:1718356/150‑1 (MQ=255)

               cggcgCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTaa                                                                                                                                        <  1:4956961/150‑1 (MQ=255)

                ggcgCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCCCCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAg                                                                                                                                       >  2:2190443/1‑150 (MQ=255)

                ggcgCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAg                                                                                                                                       <  1:1963338/150‑1 (MQ=255)

                 gcgcTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAgg                                                                                                                                      <  2:3455265/150‑1 (MQ=255)

                   gcTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGgt                                                                                                                                    >  1:2518965/1‑150 (MQ=255)

                    cTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGgtg                                                                                                                                   >  2:3268118/1‑150 (MQ=255)

                     tGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGgtgt                                                                                                                                  >  1:5444121/1‑150 (MQ=255)

                     tGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGgtgt                                                                                                                                  >  1:4720805/1‑150 (MQ=255)

                      gATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGtt                                                                                                                                 <  1:344714/150‑1 (MQ=255)

                      gATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGtt                                                                                                                                 <  1:1563361/150‑1 (MQ=255)

                       aTCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGGGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTa                                                                                                                                <  2:4920565/150‑1 (MQ=255)

                         ctcCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTAcc                                                                                                                              <  2:4997484/150‑1 (MQ=255)

                         ctcCAGGGCCGGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTAcc                                                                                                                              <  2:4925147/150‑1 (MQ=255)

                            cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTg                                                                                                                           >  1:3439455/1‑150 (MQ=255)

                            cAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTg                                                                                                                           <  2:2378307/150‑1 (MQ=255)

                               ggCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTc                                                                                                                        <  1:3513859/150‑1 (MQ=255)

                                gCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCa                                                                                                                       >  2:2008919/1‑150 (MQ=255)

                                gCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCa                                                                                                                       >  1:3920605/1‑150 (MQ=255)

                                gCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGACTCGGCACGACCGGTAAGGGTGTTACCCTGCTCa                                                                                                                       >  1:3201535/1‑150 (MQ=255)

                                 ccTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCAt                                                                                                                      >  1:553403/1‑150 (MQ=255)

                                   tGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACAGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTg                                                                                                                    >  1:2594548/1‑150 (MQ=255)

                                     gCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGTCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGAt                                                                                                                  <  2:114997/150‑1 (MQ=255)

                                      ccGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATa                                                                                                                 <  2:4285721/150‑1 (MQ=255)

                                      ccGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATa                                                                                                                 <  1:3625161/150‑1 (MQ=255)

                                      ccGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATa                                                                                                                 <  2:742977/150‑1 (MQ=255)

                                        gTCGCTTTGTGGCGGTTGCTGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGc                                                                                                               <  1:2693928/150‑1 (MQ=255)

                                        gTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGc                                                                                                               <  1:5624231/150‑1 (MQ=255)

                                         tCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGcc                                                                                                              >  1:779481/1‑150 (MQ=255)

                                          cGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCa                                                                                                             <  1:3397387/150‑1 (MQ=255)

                                           gCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAg                                                                                                            >  2:2988064/1‑150 (MQ=255)

                                           gCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAg                                                                                                            >  2:2060107/1‑150 (MQ=255)

                                                gtgGGGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAggtcgg                                                                                                       >  1:4034611/1‑150 (MQ=255)

                                                  ggCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGATACCCTGCTCATTGATAGCCAggtcggtc                                                                                                     >  2:4756163/1‑150 (MQ=255)

                                                       ttGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGTCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGgcgc                                                                                                >  2:4583232/1‑150 (MQ=255)

                                                         gCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCa                                                                                              <  2:1144996/150‑1 (MQ=255)

                                                          cAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAg                                                                                             >  2:2406774/1‑150 (MQ=255)

                                                          cAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAg                                                                                             <  2:864564/150‑1 (MQ=255)

                                                          cAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAg                                                                                             <  2:864508/150‑1 (MQ=255)

                                                            gggTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtg                                                                                           >  1:2433991/1‑150 (MQ=255)

                                                             ggTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtgg                                                                                          <  1:1007338/150‑1 (MQ=255)

                                                             ggTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCCGGCTCgac                                                                                                                                                   >  1:4043449/1‑91 (MQ=255)

                                                              gaaaTGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAgtggt                                                                                         <  2:3054697/148‑1 (MQ=255)

                                                              gTAATGCTGAACTGGCCGGCCGCACCGGCCACCGCAGTGCCCAGCACGACGCCACCTGCGCCGCGCACGGAACCCGTCGAGCCGGGCTCGGCACGACCGGTACGGGTGTTACCCTGCTCATGGAGAGCCAGGTCGGTCGGCGCCAgtggt                                                                                         <  2:3624998/150‑1 (MQ=255)

                                                                    cTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGtt                                                                                   <  1:301631/150‑1 (MQ=255)

                                                                     tGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTg                                                                                  >  2:3121174/1‑150 (MQ=255)

                                                                          tGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTcc                                                                             >  2:4734402/1‑150 (MQ=255)

                                                                           ggCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCg                                                                            >  2:1734544/1‑150 (MQ=255)

                                                                           ggCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCg                                                                            >  2:2234895/1‑150 (MQ=255)

                                                                             ccGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGt                                                                          >  2:4249279/1‑150 (MQ=255)

                                                                               gTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGAGTCCGGTGc                                                                        >  1:3101147/1‑150 (MQ=255)

                                                                                  cGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGt                                                                     >  1:5127850/1‑150 (MQ=255)

                                                                                       cGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAg                                                                >  2:4658596/1‑150 (MQ=255)

                                                                                        ggCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCACCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAgc                                                               <  2:1783504/150‑1 (MQ=255)

                                                                                         gCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAgcg                                                              >  2:2828489/1‑150 (MQ=255)

                                                                                            aCCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTa                                                                                                                                <  2:1313632/81‑1 (MQ=255)

                                                                                                cAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGgcg                                                       >  2:1689316/1‑150 (MQ=255)

                                                                                                    gCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCga                                                   <  2:4034322/150‑1 (MQ=255)

                                                                                                      ccAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCgaga                                                 >  1:1231854/1‑150 (MQ=255)

                                                                                                      ccAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCgaga                                                 >  1:2420012/1‑150 (MQ=255)

                                                                                                       cAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCAGTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGAc                                                >  1:2768696/1‑150 (MQ=255)

                                                                                                           aaTACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTa                                            <  1:3034990/150‑1 (MQ=255)

                                                                                                             tACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTAcc                                          <  1:3735201/150‑1 (MQ=255)

                                                                                                                gCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACcggc                                       >  2:4605469/1‑150 (MQ=255)

                                                                                                                        agcCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGACAGGCGCCGAGACATTACCGGCGGCgtcag                               >  2:2476910/2‑150 (MQ=255)

                                                                                                                             cGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAgg                          <  1:2139525/150‑1 (MQ=255)

                                                                                                                              gCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAggg                         <  2:4200020/150‑1 (MQ=255)

                                                                                                                                aCGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTg                       >  1:5260562/1‑150 (MQ=255)

                                                                                                                                  ggAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGtt                                                                                   <  2:3141531/88‑1 (MQ=255)

                                                                                                                                  ggAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGtt                                                                                   >  1:3141531/1‑88 (MQ=255)

                                                                                                                                       cggTCGTGCCGAGCCCGGCTCGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCaa                >  2:1528277/3‑150 (MQ=255)

                                                                                                                                       ccGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCaa                <  1:3268118/150‑1 (MQ=255)

                                                                                                                                       ccGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCaa                <  1:3658744/150‑1 (MQ=255)

                                                                                                                                        tgTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCCGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAg               >  2:3837415/2‑150 (MQ=255)

                                                                                                                                        ggTCGTGCCGAGCCCGGCTCGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAg               >  2:1519244/2‑150 (MQ=255)

                                                                                                                                        cGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAg               >  2:3343503/1‑150 (MQ=255)

                                                                                                                                         gTCGTGCCGAGCCCGGCTCGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAgg              >  2:3693595/1‑150 (MQ=255)

                                                                                                                                          tCGTGCCGAGCCCGGCTCGACGGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCaa                >  2:5333896/1‑147 (MQ=255)

                                                                                                                                          tCGTGCCGAGCCCGGCTCGACGGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCaa                <  1:5333896/147‑1 (MQ=255)

                                                                                                                                          tCGTGCCGAGCCCGGCTCGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGt             >  2:3346011/1‑150 (MQ=255)

                                                                                                                                          tCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGt             >  1:2964611/1‑150 (MQ=255)

                                                                                                                                          tCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGt             >  2:689985/1‑150 (MQ=255)

                                                                                                                                             aGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGtttt          <  1:959667/150‑1 (MQ=255)

                                                                                                                                              gCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTc         >  2:2630130/1‑150 (MQ=255)

                                                                                                                                                cGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCGc       <  2:4720805/150‑1 (MQ=255)

                                                                                                                                                 gggCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCGCt      <  2:5246765/150‑1 (MQ=255)

                                                                                                                                                   gCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCGCTGt    <  2:1916598/150‑1 (MQ=255)

                                                                                                                                                    cTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCGCTGtt   <  1:4541659/150‑1 (MQ=255)

                                                                                                                                                     tCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCGCTGTTg  <  1:4434362/150‑1 (MQ=255)

TGCCAGCCGCGTCTGCGGCGCTGATCTCCAGGGCCTGGCCGTCGCTTTGTGGCGGTTGCAGGGTAATGCTGAACTGGCCGTCCGCACCGGCCACCGCAGTGCCCAGCAATACGCCACCTGCGCCTCGCACGGAAACCGTCGAGCCGGGCTCGGCACGACCGGTAAGGGTGTTACCCTGCTCATTGATAGCCAGGTCGGTCGGCGCCAGTGGTGCAGTTGCGTCCGGTGCAGTGACAGCGCCAGGCGCCGAGACATTACCGGCGGCGTCAGTCAGGGTGACATCAAGGTTTTCGCTGTTG  >  NC_002947/944688‑944986

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