BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,786,987	T→C	5.8%	S134P (TCG→CCG) ‡	trxC →	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,786,987	0	T	C	5.8%	302.4 / 5.3	103	S134P (TCG→CCG) ‡	trxC	thioredoxin
Reads supporting (aligned to +/- strand): ref base T (51/46); new base C (3/3); total (54/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 104 total aligned reads displayed.

GCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCAC  >  NC_002947/5786839‑5787132

gCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGt                                                                                                                                                   >  6:122723/1‑149 (MQ=255)

  ttGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCg                                                                                                                                                 >  7:119582/1‑149 (MQ=255)

  ttGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCg                                                                                                                                                 >  3:345038/1‑149 (MQ=255)

  ttGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCg                                                                                                                                                 >  2:19080/1‑149 (MQ=255)

      tGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCgctgc                                                                                                                                             >  7:182522/1‑149 (MQ=255)

        gCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGctgctg                                                                                                                                           <  7:399795/149‑1 (MQ=255)

           gCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCTAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAg                                                                                                                                        <  7:90267/149‑1 (MQ=255)

               ccgccTGGCCAAGCTCGACAGCGAAGCCAACCGTACCCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGTgga                                                                                                                                    <  5:287137/149‑2 (MQ=255)

                cgccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGTGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCt                                                                                                                                   <  2:302088/149‑1 (MQ=255)

                 gccTGGCCAAGCTCGACAGTGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTg                                                                                                                                  >  7:203933/1‑149 (MQ=255)

                 gccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGt                                                                                                                                       >  2:46682/1‑144 (MQ=255)

                 gccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGt                                                                                                                                       <  1:46682/144‑1 (MQ=255)

                   cTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCg                                                                                                                                >  2:232511/1‑149 (MQ=255)

                        cAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGcc                                                                                                                           <  4:345038/149‑1 (MQ=255)

                        cAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGcc                                                                                                                           <  6:96820/149‑1 (MQ=255)

                           gCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAg                                                                                                                        <  1:232511/149‑1 (MQ=255)

                            cTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAgg                                                                                                                       >  3:52518/1‑149 (MQ=255)

                              cGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAgggg                                                                                                                     >  7:176642/1‑149 (MQ=255)

                               gACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCATGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGt                                                                                                                    <  3:64661/149‑1 (MQ=255)

                               gACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGc                                                                                                                                 <  5:259297/136‑1 (MQ=255)

                               gACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGc                                                                                                                                 >  6:259297/1‑136 (MQ=255)

                               gACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAggg                                                                                                                      <  2:421028/147‑1 (MQ=255)

                               gACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGt                                                                                                                    >  6:234089/1‑149 (MQ=255)

                               gACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCACGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGt                                                                                                                    >  1:310261/1‑149 (MQ=255)

                                 cAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGtt                                                                                                                   >  6:247654/1‑148 (MQ=255)

                                  aGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGttt                                                                                                                  >  1:425653/1‑148 (MQ=255)

                                        gCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGtt                                                                                                           >  8:84321/1‑149 (MQ=255)

                                        gCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCACCCGCCACGCCGGGGCCTTTCCACTCCACTCGCTGCTCGACTGGCTCCGTAGCAAAGGCGTTTAAGAGtg                                                                                                           >  4:348661/1‑148 (MQ=255)

                                        gCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCCGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGtt                                                                                                           >  7:35927/1‑149 (MQ=255)

                                         ccAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGttt                                                                                                          >  2:280668/1‑149 (MQ=255)

                                         ccAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGTgg                                                                                                                                     >  5:310971/1‑122 (MQ=255)

                                         ccAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGTgg                                                                                                                                     <  6:310971/122‑1 (MQ=255)

                                           aaCCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttc                                                                                                        >  3:140501/1‑149 (MQ=255)

                                             ccGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttctt                                                                                                      >  8:75302/1‑149 (MQ=255)

                                              cGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttcttc                                                                                                     <  7:75302/149‑1 (MQ=255)

                                              cGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttctt                                                                                                      >  6:135171/1‑148 (MQ=255)

                                               gTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCa                                                                                                    <  2:310261/149‑1 (MQ=255)

                                                tAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCaa                                                                                                   <  2:289668/149‑1 (MQ=255)

                                                tAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCaa                                                                                                   <  2:251070/149‑1 (MQ=255)

                                                  aCCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAag                                                                                                  <  4:140501/148‑1 (MQ=255)

                                                   ccTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAag                                                                                                  >  8:251209/1‑147 (MQ=255)

                                                         aGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAAt                                                                                          >  6:197658/1‑149 (MQ=255)

                                                         aGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAAt                                                                                          <  5:234089/149‑1 (MQ=255)

                                                          gggCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATg                                                                                         <  5:212451/149‑1 (MQ=255)

                                                          gggCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTAAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATg                                                                                         <  4:67318/149‑1 (MQ=255)

                                                             cAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAg                                                                                      <  7:77628/149‑1 (MQ=255)

                                                              aaCTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGc                                                                                     >  7:113931/1‑149 (MQ=255)

                                                                cTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTc                                                                                   <  8:385386/149‑1 (MQ=255)

                                                                   ggCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGa                                                                                 >  5:434593/1‑148 (MQ=255)

                                                                   ggCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGAc                                                                                <  8:203933/149‑1 (MQ=255)

                                                                   ggCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGAc                                                                                <  7:39030/149‑1 (MQ=255)

                                                                   ggCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGAc                                                                                <  6:434593/149‑1 (MQ=255)

                                                                      aTACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACaaa                                                                             >  8:109497/1‑149 (MQ=255)

                                                                       tACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAAc                                                                            >  2:224372/1‑149 (MQ=255)

                                                                       tACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAAc                                                                            >  4:395933/1‑149 (MQ=255)

                                                                         cGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActg                                                                          <  4:311653/149‑1 (MQ=255)

                                                                         cGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAAct                                                                           <  7:251209/148‑1 (MQ=255)

                                                                          gCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgc                                                                         >  8:212567/1‑149 (MQ=255)

                                                                          gCTCGATCACCAGCCTGCTGCTGTTCAAGAACGGACGAGAAATCACCCGACAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgc                                                                         >  1:11402/1‑149 (MQ=255)

                                                                           cTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgct                                                                        >  1:68182/1‑149 (MQ=255)

                                                                            tCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCTGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgct                                                                        >  2:199488/1‑148 (MQ=255)

                                                                              gATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGa                                                                      <  8:26619/148‑1 (MQ=255)

                                                                              gATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGa                                                                      <  4:52518/148‑1 (MQ=255)

                                                                                 ccccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGt                                                                    >  1:25808/1‑147 (MQ=255)

                                                                                 ccccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGt                                                                    >  5:363653/1‑147 (MQ=255)

                                                                                  cccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGt                                                                                                                                       <  7:77024/79‑1 (MQ=255)

                                                                                  cccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGt                                                                                                                                       >  8:77024/1‑79 (MQ=255)

                                                                                   ccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGc                                                                <  8:176642/149‑1 (MQ=255)

                                                                                   ccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGc                                                                >  2:421136/1‑149 (MQ=255)

                                                                                       ccTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTtgt                                                            <  8:113931/149‑1 (MQ=255)

                                                                                                  tCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActg                                                                          >  7:93807/1‑124 (MQ=255)

                                                                                                    aaGAACGGCCGCGAAATCAGCCGCCGGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAg                                               <  6:432286/149‑1 (MQ=255)

                                                                                                    aaGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAg                                               <  1:224372/149‑1 (MQ=255)

                                                                                                    aaGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAg                                               >  3:367608/1‑149 (MQ=255)

                                                                                                       aaCGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAg                                                                 >  8:55091/1‑128 (MQ=255)

                                                                                                       aaCGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAg                                                                 <  7:55091/128‑1 (MQ=255)

                                                                                                         cGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGcc                                                         <  6:164133/134‑1 (MQ=255)

                                                                                                         cGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGcc                                                         >  5:164133/1‑134 (MQ=255)

                                                                                                         cGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGaa                                          <  5:135171/149‑1 (MQ=255)

                                                                                                            ccGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAAt                                         >  6:13750/1‑147 (MQ=255)

                                                                                                             cgcgAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAg                                      >  1:390262/1‑149 (MQ=255)

                                                                                                                gAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                   >  3:35407/1‑149 (MQ=255)

                                                                                                                 aaaTCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGc                                  <  3:395933/149‑1 (MQ=255)

                                                                                                                  aaTCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                   >  7:101195/1‑147 (MQ=255)

                                                                                                                    tCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                   >  6:358432/1‑145 (MQ=255)

                                                                                                                    tCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                   <  5:358432/145‑1 (MQ=255)

                                                                                                                         cgccAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGc                          <  8:35927/149‑1 (MQ=255)

                                                                                                                            cAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtg                       <  1:280668/149‑1 (MQ=255)

                                                                                                                               gCCGGGGCGTTTCCACTCCAGTCTCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtggtg                    <  8:119582/149‑1 (MQ=255)

                                                                                                                               gCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtggtg                    <  2:390262/149‑1 (MQ=255)

                                                                                                                                ctggtgCGTTTCCACTCCAGTCGCTGCTGTATTGGCTGCTTATCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtggtg                    <  5:122723/143‑1 (MQ=255)

                                                                                                                                    ggCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGAtt               <  7:109497/149‑1 (MQ=255)

                                                                                                                                       gTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgc            <  8:247381/149‑1 (MQ=255)

                                                                                                                                        tttCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgcg           >  2:280030/1‑149 (MQ=255)

                                                                                                                                          gccACTCCAGCAGCGACTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCa         >  5:248139/2‑149 (MQ=255)

                                                                                                                                               tCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtttt     >  7:238450/1‑148 (MQ=255)

                                                                                                                                               tCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtttt     >  7:20658/1‑148 (MQ=255)

                                                                                                                                               tCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtttt     >  7:14292/1‑148 (MQ=255)

                                                                                                                                               tCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtttt     >  6:441852/1‑148 (MQ=255)

                                                                                                                                                 cAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCAc  <  1:199488/149‑1 (MQ=255)

GCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCAC  >  NC_002947/5786839‑5787132

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