BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,791,997	G→A	5.9%	L135L (CTG→CTA) ‡	trxC →	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,791,997	0	G	A	5.9%	1184.6 / 37.5	374	L135L (CTG→CTA) ‡	trxC	thioredoxin
Reads supporting (aligned to +/- strand): ref base G (160/192); new base A (13/9); total (173/201)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 385 total aligned reads displayed.

CTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGCTGC  >  NC_002947/5791849‑5792145

 tGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCgctgct                                                                                                                                                    >  2:1616595/1‑150 (MQ=255)

 tGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCgctgct                                                                                                                                                    >  1:6180107/1‑150 (MQ=255)

  ggCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGctgctg                                                                                                                                                   <  2:5997200/150‑1 (MQ=255)

     cGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAg                                                                                                                                                <  2:4009023/150‑1 (MQ=255)

      gCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGt                                                                                                                                               <  1:7047887/150‑1 (MQ=255)

          ccgccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCt                                                                                                                                           >  2:1237219/1‑150 (MQ=255)

          ccgccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCACGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGACCGGCCGCGAAAGCAGCCGCCCGGCCGGGGCGTTTCCACGCCCGGCGCTGCGGGAGTGGCt                                                                                                                                           >  2:5646752/1‑150 (MQ=255)

            gccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGc                                                                                                                                         >  1:4042116/1‑150 (MQ=255)

            gccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGc                                                                                                                                         >  1:3145844/1‑150 (MQ=255)

             ccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCg                                                                                                                                        >  2:4419949/1‑150 (MQ=255)

              cTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTGCCACTCCAGCAGCGACTGGAGTggaaacgc                                                                                                                                       <  2:5804528/150‑7 (MQ=255)

               tGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTa                                                                                                                                      <  1:7527287/150‑1 (MQ=255)

                  ccAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGcc                                                                                                                                   >  2:4911032/1‑150 (MQ=255)

                  ccAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTGTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGcc                                                                                                                                   >  2:6469054/1‑150 (MQ=255)

                    aaGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCaa                                                                                                                                 <  1:3317122/150‑1 (MQ=255)

                      gCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAgg                                                                                                                               <  1:5844822/150‑1 (MQ=255)

                      gCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAgg                                                                                                                               <  2:2199479/150‑1 (MQ=255)

                       cTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAggg                                                                                                                              <  1:7240962/150‑1 (MQ=255)

                         cGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGt                                                                                                                            <  2:4252884/150‑1 (MQ=255)

                          gACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGtt                                                                                                                           >  1:4873582/1‑150 (MQ=255)

                           aCAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGACCGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGTGTGGCTGCGTAGCCAAGGGGttt                                                                                                                          <  2:8210409/150‑1 (MQ=255)

                             aGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTaa                                                                                                                        <  2:1522142/150‑1 (MQ=255)

                             aGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTaa                                                                                                                        <  1:5324832/150‑1 (MQ=255)

                               cGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGc                                                                                                                      <  2:6196672/150‑1 (MQ=255)

                                gAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCg                                                                                                                     >  2:3177962/1‑150 (MQ=255)

                                 aaGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCTAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGt                                                                                                                    >  2:7348071/1‑150 (MQ=255)

                                 aaGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGt                                                                                                                    <  1:2039026/150‑1 (MQ=255)

                                  aGCCAACCGTAACCTCGCTGGGCAACTGGTCATACGCTCGATCCCCAGCCTGCTGCTGTTCAATAACGTCCGCGTAATCTGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGAGTAGCCAAGGGGTTTAAGCGtt                                                                                                                   <  1:1825224/150‑1 (MQ=255)

                                  aGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCATGGGGTTTAAGCGtt                                                                                                                   >  1:8292531/1‑150 (MQ=255)

                                    ccAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACAGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGtttt                                                                                                                 >  2:1758975/1‑150 (MQ=255)

                                     cAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGCACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttc                                                                                                                <  1:2538320/150‑1 (MQ=255)

                                      aaCCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttct                                                                                                               >  1:3834171/1‑150 (MQ=255)

                                      aaCCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttct                                                                                                               >  1:798743/1‑150 (MQ=255)

                                      aaCCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttct                                                                                                               >  2:7490782/1‑150 (MQ=255)

                                      aaCCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttct                                                                                                               >  1:2634099/1‑150 (MQ=255)

                                        ccGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttcttc                                                                                                             <  2:2714305/150‑1 (MQ=255)

                                            aaCCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAagc                                                                                                         <  1:7122390/150‑1 (MQ=255)

                                                         aaCTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCt                                                                                            <  1:2140733/150‑1 (MQ=255)

                                                         aaCTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCt                                                                                            >  1:5486777/1‑150 (MQ=255)

                                                               gCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACaa                                                                                      >  1:8740877/1‑150 (MQ=255)

                                                                 aTACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAAc                                                                                    <  2:3890906/150‑1 (MQ=255)

                                                                 aTACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAAc                                                                                    >  1:4451157/1‑150 (MQ=255)

                                                                  tACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCACGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACACAct                                                                                   >  2:1985960/1‑150 (MQ=255)

                                                                  tACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAg                                                                                                                                >  1:1954100/1‑105 (MQ=255)

                                                                  tACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCgtcgccac                                                                                   >  1:109259/1‑143 (MQ=255)

                                                                  tACGCACGATCCCCAGCCTGCAGCAGTTCAAGAACCGCCGCGAAATCAGCCGCCAGGCCGGGGCGGTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAAct                                                                                   >  2:5495226/1‑150 (MQ=255)

                                                                   aCGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActg                                                                                  >  2:177841/1‑150 (MQ=255)

                                                                    cGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgc                                                                                 <  1:8231527/150‑1 (MQ=255)

                                                                    cGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgc                                                                                 <  1:6090651/150‑1 (MQ=255)

                                                                     gCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgct                                                                                <  1:6571333/150‑1 (MQ=255)

                                                                       tCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGa                                                                              <  2:8675971/150‑1 (MQ=255)

                                                                       tCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGa                                                                              <  1:3083082/150‑1 (MQ=255)

                                                                       tCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGa                                                                              <  1:5748855/150‑1 (MQ=255)

                                                                         gATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGt                                                                            <  2:7154602/150‑1 (MQ=255)

                                                                         gATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGt                                                                            <  2:6274508/150‑1 (MQ=255)

                                                                             cccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCCAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGc                                                                        >  2:6323356/1‑150 (MQ=255)

                                                                              ccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCt                                                                       >  1:879016/1‑150 (MQ=255)

                                                                               cAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCtt                                                                      <  1:3588150/150‑1 (MQ=255)

                                                                               cAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCtt                                                                      >  1:3417033/1‑150 (MQ=255)

                                                                                     gctgctgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCcg                                                                >  2:5973376/1‑150 (MQ=255)

                                                                                        gctgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCgg                                                             <  1:6523790/150‑1 (MQ=255)

                                                                                         ctgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGc                                                            <  2:2931114/150‑1 (MQ=255)

                                                                                         ctgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGc                                                            <  2:3725474/150‑1 (MQ=255)

                                                                                          tgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCt                                                           >  2:7814834/1‑150 (MQ=255)

                                                                                            ttCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTcc                                                         >  2:6630559/1‑150 (MQ=255)

                                                                                               aaGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAgg                                                      <  2:8089896/150‑1 (MQ=255)

                                                                                                 gAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTg                                                    >  2:5318004/1‑150 (MQ=255)

                                                                                                  aaCGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGa                                                   >  1:755385/1‑150 (MQ=255)

                                                                                                  aaCGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGGGTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGa                                                   >  2:6955613/1‑150 (MQ=255)

                                                                                                     ggCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATc                                                >  2:7532065/1‑150 (MQ=255)

                                                                                                       ccGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGAAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAg                                              >  2:5389064/1‑150 (MQ=255)

                                                                                                        cgcgAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGc                                             >  1:3727203/1‑150 (MQ=255)

                                                                                                        cgcgAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGc                                             >  2:1093948/1‑150 (MQ=255)

                                                                                                         gcgAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCg                                            <  2:6066328/150‑1 (MQ=255)

                                                                                                          cgAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                           <  2:2563082/150‑1 (MQ=255)

                                                                                                          cgAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                           >  1:7388392/1‑150 (MQ=255)

                                                                                                            aaaTCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCa                                         <  1:1459372/150‑1 (MQ=255)

                                                                                                                  gccgccAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTgg                                   >  1:2680459/1‑150 (MQ=255)

                                                                                                                     gccAGGCCGGGGCGTTCCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgt                                <  2:233310/150‑1 (MQ=255)

                                                                                                                       cAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtgg                              >  2:5763091/1‑150 (MQ=255)

                                                                                                                        aGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtggt                             >  1:5049360/1‑150 (MQ=255)

                                                                                                                         ggCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtggtg                            <  1:7812061/150‑1 (MQ=255)

                                                                                                                           ccGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCg                          <  1:93572/150‑1 (MQ=255)

                                                                                                                            cGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGa                         <  1:4001316/150‑1 (MQ=255)

                                                                                                                             ggggCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGAt                        <  2:7934852/150‑1 (MQ=255)

                                                                                                                               ggCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTc                      <  1:4652220/150‑1 (MQ=255)

                                                                                                                                 cGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgc                    <  2:7089102/150‑1 (MQ=255)

                                                                                                                                 cGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgc                    <  1:6698005/150‑1 (MQ=255)

                                                                                                                                   tttCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCgcgc                  >  2:2812833/1‑150 (MQ=255)

                                                                                                                                        aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTACGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCGGAGTCAGCTTGTGCCGCGGCTCCAGGTGACTCAGCGGCAGGCTGGCGGGGTGCGATTGGCGCAtttt             >  2:8457672/1‑150 (MQ=255)

                                                                                                                                        aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtttt             >  1:3169285/1‑150 (MQ=255)

                                                                                                                                        aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtttt             >  1:1724163/1‑150 (MQ=255)

                                                                                                                                        aCTCCAGCAGCGACTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtttt             >  2:3479672/1‑150 (MQ=255)

                                                                                                                                           ccAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCAc          <  1:1361643/150‑1 (MQ=255)

                                                                                                                                           ccAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCAc          >  2:3640606/1‑150 (MQ=255)

                                                                                                                                            cAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCAcc         >  2:5566504/1‑150 (MQ=255)

                                                                                                                                               tCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAg      <  1:4333253/150‑1 (MQ=255)

                                                                                                                                                cGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGc     <  2:3323643/150‑1 (MQ=255)

                                                                                                                                                 gctgctGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGCt    >  2:4664577/1‑150 (MQ=255)

                                                                                                                                                   tgctgGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGCTGc  <  1:2244551/150‑1 (MQ=255)

CTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGCTGC  >  NC_002947/5791849‑5792145

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