BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,299,502	A→T	5.7%	L862Q (CTG→CAG)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,299,502	0	A	T	5.7%	613.0 / 17.5	182	L862Q (CTG→CAG)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand): ref base A (77/94); new base T (3/7); total (81/101)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 182 total aligned reads displayed.

TGGCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTC  >  NC_002947/2299356‑2299650

  gCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCg                                                                                                                                                 <  1:506058/150‑1 (MQ=255)

  gCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCg                                                                                                                                                 >  2:4141364/1‑150 (MQ=255)

  gCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCg                                                                                                                                                 >  1:3292576/1‑150 (MQ=255)

   ccTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCgg                                                                                                                                                <  1:5206006/150‑1 (MQ=255)

   ccTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCgg                                                                                                                                                <  1:5205262/150‑1 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                               <  2:1224329/150‑1 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                               <  2:2607970/150‑1 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                               >  1:965689/1‑150 (MQ=255)

    cTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGc                                                                                                                                               <  1:1001075/150‑1 (MQ=255)

      tGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGa                                                                                                                                             <  1:411598/150‑1 (MQ=255)

      tGAGCTCGCTGGCCAGTTGCACGCCTCGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGa                                                                                                                                             <  1:906636/150‑1 (MQ=255)

       gAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCGGCGaa                                                                                                                                            <  2:3212145/150‑2 (MQ=255)

        aGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGt                                                                                                                                           <  2:5370134/150‑1 (MQ=255)

         gCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGGGTCAAGGGTGCGGGCGTGCTCCAGCGGGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCAGGTGCGCCGAGCGGTTCGGCGCGTg                                                                                                                                          >  2:5022344/1‑150 (MQ=255)

         gCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTg                                                                                                                                          >  2:5022380/1‑150 (MQ=255)

           tCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAACTGCTCGGCGaacatg                                                                                                                                        <  1:1322508/150‑6 (MQ=255)

               tGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCa                                                                                                                                    >  2:482997/1‑150 (MQ=255)

                ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                   <  1:4088807/150‑1 (MQ=255)

                ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                   <  1:1660314/150‑1 (MQ=255)

                ggCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAg                                                                                                                                   <  1:1460122/150‑1 (MQ=255)

                 gCCAGTTGCACGCCTTGGGTGTGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgg                                                                                                                                  <  2:5424452/150‑1 (MQ=255)

                 gCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgg                                                                                                                                  >  2:2191313/1‑150 (MQ=255)

                 gCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgg                                                                                                                                  <  2:4408398/150‑1 (MQ=255)

                   cAGTTGCACGCGTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAgggg                                                                                                                                >  2:155140/1‑149 (MQ=255)

                    aGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCt                                                                                                                               <  2:4515129/150‑1 (MQ=255)

                       tGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGc                                                                                                                            <  1:12579/150‑1 (MQ=255)

                            gCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCTCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCGTGCTGGCg                                                                                                                       >  2:2944295/1‑150 (MQ=255)

                             ccTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGt                                                                                                                      <  1:4048978/150‑1 (MQ=255)

                                tGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGc                                                                                                                   <  2:191188/150‑1 (MQ=255)

                                   gTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCtcg                                                                                                                >  2:1648706/1‑150 (MQ=255)

                                         aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAgg                                                                                                          <  1:1026282/150‑1 (MQ=255)

                                         aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAgg                                                                                                          <  1:2685790/150‑1 (MQ=255)

                                         aGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTCGCTCGTCGAgg                                                                                                          <  2:1953590/150‑1 (MQ=255)

                                            tCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGggcg                                                                                                       >  1:1325286/1‑150 (MQ=255)

                                            tCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGggcg                                                                                                       <  2:5281137/150‑1 (MQ=255)

                                                  aGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGt                                                                                                 >  2:1736270/1‑150 (MQ=255)

                                                    gcgcTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGttt                                                                                               <  2:2832923/150‑1 (MQ=255)

                                                         gggCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctgct                                                                                          <  1:5253519/150‑1 (MQ=255)

                                                          ggCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctgctg                                                                                         <  2:3685844/150‑1 (MQ=255)

                                                            cGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTa                                                                                       >  2:4950174/1‑150 (MQ=255)

                                                                   ggCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCa                                                                                                   <  2:3165777/131‑1 (MQ=255)

                                                                   ggCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCa                                                                                                   >  1:3165777/1‑131 (MQ=255)

                                                                       tCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctg                                                                                            <  1:4695392/134‑1 (MQ=255)

                                                                       tCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTctg                                                                                            >  2:4695392/1‑134 (MQ=255)

                                                                          cgggTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCCTGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGGTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCt                                                                         >  2:3564396/2‑150 (MQ=255)

                                                                             gTGCGGGCGTGCTCCAGCGCGGGGTCCATGTGCTGCTGCCAGGCTTCGGAGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGc                                                                      >  1:4453675/1‑150 (MQ=255)

                                                                                  ggCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAggg                                                                                                         <  1:2405601/110‑1 (MQ=255)

                                                                                   gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCtt                                                                        <  1:5353415/142‑1 (MQ=255)

                                                                                   gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCtt                                                                        >  2:5353415/1‑142 (MQ=255)

                                                                                   gCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCAcgcgc                                                                >  1:2390717/1‑150 (MQ=255)

                                                                                     gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGt                                                                                                                                           <  2:14610/73‑1 (MQ=255)

                                                                                     gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACgcgcgc                                                              >  2:2029055/1‑150 (MQ=255)

                                                                                     gTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACgcgcgc                                                              <  1:5022344/150‑1 (MQ=255)

                                                                                         tCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgc                                                          >  2:3828169/1‑150 (MQ=255)

                                                                                           cAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCctgctg                                                        <  2:766073/150‑1 (MQ=255)

                                                                                            aGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGt                                                       >  1:22960/1‑150 (MQ=255)

                                                                                             gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGTTGGCCTTGCACGCGCGCctgcaaga                                                      <  2:1325286/150‑5 (MQ=255)

                                                                                             gcgcGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGtt                                                      <  1:4228702/150‑1 (MQ=255)

                                                                                              cgcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCgg                                                                                                                                                >  2:2274192/1‑59 (MQ=255)

                                                                                              cgcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTg                                                     <  1:1073516/150‑1 (MQ=255)

                                                                                               gcgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTgg                                                    <  2:1116453/150‑1 (MQ=255)

                                                                                                cgGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGt                                                   <  1:1972554/150‑1 (MQ=255)

                                                                                                     tCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGc                                              >  2:2439964/1‑150 (MQ=255)

                                                                                                       cATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAg                                            >  2:1038568/1‑150 (MQ=255)

                                                                                                         tgtgCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCt                                          >  1:4338352/1‑150 (MQ=255)

                                                                                                             ctgctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTc                                      >  2:5391951/1‑150 (MQ=255)

                                                                                                             ctgctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTc                                      <  1:1736270/150‑1 (MQ=255)

                                                                                                              tgctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTcc                                     >  2:2084615/1‑150 (MQ=255)

                                                                                                               gctgcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCa                                    >  1:884913/1‑150 (MQ=255)

                                                                                                                  gcCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCt                                 <  1:928361/150‑1 (MQ=255)

                                                                                                                   ccAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTg                                >  1:5307895/1‑150 (MQ=255)

                                                                                                                     aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                              >  2:2385522/1‑150 (MQ=255)

                                                                                                                     aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                              >  1:3673073/1‑150 (MQ=255)

                                                                                                                     aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                              >  2:3989962/1‑150 (MQ=255)

                                                                                                                     aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                              >  2:712495/1‑150 (MQ=255)

                                                                                                                     aGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTGTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTggg                              >  2:4188119/1‑150 (MQ=255)

                                                                                                                      ggCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGt                             <  2:2719485/150‑1 (MQ=255)

                                                                                                                      ggCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGt                             <  1:482786/150‑1 (MQ=255)

                                                                                                                         ttCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGGGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTg                          >  1:1985078/1‑150 (MQ=255)

                                                                                                                         ttCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTg                          >  1:1885490/1‑150 (MQ=255)

                                                                                                                            ggCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTg                                                     <  2:3213682/120‑1 (MQ=255)

                                                                                                                            ggCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTg                                                     >  1:3213682/1‑120 (MQ=255)

                                                                                                                             gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                      <  1:1648595/150‑1 (MQ=255)

                                                                                                                             gcgcTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGGTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCgg                      >  2:3332236/1‑150 (MQ=255)

                                                                                                                                 tGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCtt                  <  2:884913/150‑1 (MQ=255)

                                                                                                                                   gCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgc                <  1:3628450/150‑1 (MQ=255)

                                                                                                                                      tGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTtgctgc             <  2:4465892/150‑1 (MQ=255)

                                                                                                                                       gTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct            >  1:4479351/1‑150 (MQ=255)

                                                                                                                                       gTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct            >  1:1098756/1‑150 (MQ=255)

                                                                                                                                       actCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTgctgct            >  1:5391542/3‑150 (MQ=255)

                                                                                                                                         tCGCCGAACTGCTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGctct          >  2:2299527/1‑150 (MQ=255)

                                                                                                                                           gCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTc        <  2:2178622/150‑1 (MQ=255)

                                                                                                                                           gCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTc        >  2:1126926/1‑150 (MQ=255)

                                                                                                                                            ccGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTcc       <  2:3292576/150‑1 (MQ=255)

                                                                                                                                             cGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCa      <  1:4950174/150‑1 (MQ=255)

                                                                                                                                             cGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACACAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCa      <  1:761182/150‑1 (MQ=255)

                                                                                                                                               aGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGt    <  1:1762536/150‑1 (MQ=255)

                                                                                                                                                 cAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTc  <  2:5307991/150‑1 (MQ=255)

                                                                                                                                                 cAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTc  <  1:4753918/150‑1 (MQ=255)

                                                                                                                                                 cAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTc  <  2:867989/150‑1 (MQ=255)

TGGCCTTGAGCTCGCTGGCCAGTTGCACGCCTTGGGTGCGCAGATCTTGCAGGCGCTGGGCGGTGTCGGCGTCAAGGGTGCGGGCGTGCTCCAGCGCGGTGTCCATGTGCTGCTGCCAGGCTTCGGCGCTGGCATGTTCGCCGAGCAGTTCGGCGAGTGAGGCCAGGGCTTGCTGGCGTTGCTCGTCGAGGGCGGCCAGTTTCTGCTGTAACTGCTGCTGGCCTTGCACGCGCGCCTGCTGTTGGTCACGCAGCTTGTCCAGCTGGGTTTGGCGGGCTTGCTGCTCTTCCAGTTC  >  NC_002947/2299356‑2299650

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