BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,765,285	T→G	5.7%	S136S (TCT→TCG)	PP_1574 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,765,285	0	T	G	5.7%	708.0 / 21.8	204	S136S (TCT→TCG)	PP_1574	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (108/84); new base G (7/5); total (115/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

CGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGC  >  NC_002947/1765137‑1765434

cGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCACGCGCCTGGGATCTCTTGATGGAGACCTGTCCGTCGCGATGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGctctt                                                                                                                                                      <  1:4051963/150‑2 (MQ=255)

cGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTg                                                                                                                                                      >  1:2764305/1‑150 (MQ=255)

  cTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGc                                                                                                                                                    >  1:3825657/1‑150 (MQ=255)

     gACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACg                                                                                                                                                 <  1:1955052/150‑1 (MQ=255)

      aCCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGa                                                                                                                                                <  1:3622181/150‑1 (MQ=255)

       cccGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAg                                                                                                                                               <  2:4001667/150‑1 (MQ=255)

          ggTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTc                                                                                                                                            <  2:2985260/150‑1 (MQ=255)

             tGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTc                                                                                                                                         <  2:4450094/150‑1 (MQ=255)

               aaCTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCGTGCCAGAGCTCGTcga                                                                                                                                       <  1:2829296/150‑3 (MQ=255)

                 cTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCctct                                                                                                                                     <  2:3006200/150‑1 (MQ=255)

                   ggAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTcc                                                                                                                                   >  1:4068712/1‑150 (MQ=255)

                         ggACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATc                                                                                                                             >  2:1330786/1‑150 (MQ=255)

                         ggACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCGATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATc                                                                                                                             >  2:4301422/1‑150 (MQ=255)

                              cGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCGGGGACAGACGAAGAGCTCACCGAGCTCGCCGAGCTCTGGCACGAGCTCGTGCTCTCCGGCCTcggcgg                                                                                                                        >  2:4037065/1‑150 (MQ=255)

                                   agCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcga                                                                                                                   <  1:2894465/150‑1 (MQ=255)

                                   agCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcga                                                                                                                   <  2:3795247/150‑1 (MQ=255)

                                     cACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcgacg                                                                                                                 >  1:2991687/1‑150 (MQ=255)

                                        cAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATc                                                                                                              >  2:3045032/1‑150 (MQ=255)

                                        cAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATc                                                                                                              >  1:4034229/1‑150 (MQ=255)

                                        cAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATc                                                                                                              <  1:3836659/150‑1 (MQ=255)

                                         aaGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATcg                                                                                                             >  1:1555827/1‑150 (MQ=255)

                                             gcCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCa                                                                                                         <  1:870901/150‑1 (MQ=255)

                                              ccTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                        <  2:3876742/150‑1 (MQ=255)

                                              ccTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                        <  2:3601619/150‑1 (MQ=255)

                                              ccTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                        <  1:771412/150‑1 (MQ=255)

                                              ccTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                        >  2:214307/1‑150 (MQ=255)

                                               cTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAgg                                                                                                       >  1:4369375/1‑150 (MQ=255)

                                                tGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGc                                                                                                      >  2:1078173/1‑150 (MQ=255)

                                                   gATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCaa                                                                                                   <  2:2310081/150‑1 (MQ=255)

                                                    aTCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCaaa                                                                                                  <  2:2585065/150‑1 (MQ=255)

                                                    aTCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCaaa                                                                                                  <  2:173395/150‑1 (MQ=255)

                                                    aTCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCaaa                                                                                                  <  1:1712754/150‑1 (MQ=255)

                                                    aTCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCaaa                                                                                                  <  2:3338879/150‑1 (MQ=255)

                                                     tctcTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAg                                                                                                 <  2:2861136/150‑1 (MQ=255)

                                                        ctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGTCGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaaa                                                                                              >  2:3020115/1‑150 (MQ=255)

                                                         ttGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGaaaa                                                                                             >  1:601419/1‑150 (MQ=255)

                                                          tGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAAc                                                                                            >  1:2908442/1‑150 (MQ=255)

                                                             tGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTg                                                                                         >  1:342100/1‑150 (MQ=255)

                                                              ggcgtgcTGTCCGTCGCGCTGCTGGTAGGTATTCGTGAGGTGATCGAGCTGGGTTTGACGAAGAGCTCTCCGTGGTCGACGAGCTCTGGCTCGAGCTGGTCCTGTCCGGAAGCGGCGGGTCGACGATGGGGCTGGGCAAAGTAAAGCTGa                                                                                        <  2:3511111/144‑1 (MQ=255)

                                                              ggCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGTGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                        <  2:2146677/150‑1 (MQ=255)

                                                               gCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAg                                                                                       >  2:3731981/1‑150 (MQ=255)

                                                                cGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCG‑‑‑‑‑‑‑‑TGCCA‑GAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAg                                                                             <  1:3939120/150‑1 (MQ=255)

                                                                     tGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCg                                                                                 <  1:2122805/150‑1 (MQ=255)

                                                                     tGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCg                                                                                 >  2:14628/1‑150 (MQ=255)

                                                                     tGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCg                                                                                 >  1:1494294/1‑150 (MQ=255)

                                                                       tCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTg                                                                               >  2:1647809/1‑150 (MQ=255)

                                                                       tCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTg                                                                               >  2:4086545/1‑150 (MQ=255)

                                                                           tCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                        <  2:1824339/121‑1 (MQ=255)

                                                                           tCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAg                                                                                                        >  1:1824339/1‑121 (MQ=255)

                                                                           tCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGtt                                                                           >  2:1115538/1‑150 (MQ=255)

                                                                           tCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCAGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGTGCCTGCGTGAGtt                                                                           >  2:315522/1‑150 (MQ=255)

                                                                               gctgctGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAgg                                                                       >  1:416771/1‑150 (MQ=255)

                                                                                 tgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTc                                                                     >  2:2114078/1‑150 (MQ=255)

                                                                                 tgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTc                                                                     >  1:1016748/1‑150 (MQ=255)

                                                                                    tgGTAGCTCTTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTg                                                                  <  1:983237/150‑1 (MQ=255)

                                                                                     ggTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTgg                                                                 <  1:3811017/150‑1 (MQ=255)

                                                                                      gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCTAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                >  2:1039785/1‑150 (MQ=255)

                                                                                      gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                >  2:1039731/1‑150 (MQ=255)

                                                                                      gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTggg                                                                <  1:1916237/150‑1 (MQ=255)

                                                                                       tAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGga                                                               <  2:2649267/150‑1 (MQ=255)

                                                                                         gCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGgaga                                                             >  2:779888/1‑150 (MQ=255)

                                                                                         gCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGgaga                                                             >  1:3161873/1‑150 (MQ=255)

                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                        <  1:2596694/122‑1 (MQ=255)

                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                            <  1:134753/150‑1 (MQ=255)

                                                                                          cTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                            <  1:2253566/150‑1 (MQ=255)

                                                                                           tATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAg                                                           >  2:1373204/1‑150 (MQ=255)

                                                                                           tATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAg                                                           >  2:2618034/1‑150 (MQ=255)

                                                                                           tATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTGGGGAGAAg                                                           >  2:2339769/1‑150 (MQ=255)

                                                                                            aTTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGt                                                          >  2:1365271/1‑150 (MQ=255)

                                                                                             ttCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTa                                                         >  2:455295/1‑150 (MQ=255)

                                                                                             ttCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTa                                                         >  1:2554953/1‑150 (MQ=255)

                                                                                                gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                      <  2:2004894/150‑1 (MQ=255)

                                                                                                gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                      >  1:1263837/1‑150 (MQ=255)

                                                                                                gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAg                                                      <  1:2961237/150‑1 (MQ=255)

                                                                                                 ggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAgg                                                     <  2:601419/150‑1 (MQ=255)

                                                                                                 ggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAgg                                                     >  2:2474761/1‑150 (MQ=255)

                                                                                                      tCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCt                                                <  2:2628052/150‑1 (MQ=255)

                                                                                                        gaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGGCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAGGTACAGGCAGCTcc                                              >  1:4303789/2‑150 (MQ=255)

                                                                                                        aaCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTcc                                              >  2:3440134/1‑150 (MQ=255)

                                                                                                               tGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTc                                       <  2:1603272/150‑1 (MQ=255)

                                                                                                                  gAAAGACGAAGAGCTCACCGAGCTCGTCGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCt                                    <  1:3020503/150‑1 (MQ=255)

                                                                                                                  gAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCt                                    >  1:2377524/1‑150 (MQ=255)

                                                                                                                  gAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCt                                    >  2:3952302/1‑150 (MQ=255)

                                                                                                                      gACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCAt                                >  2:1030484/1‑150 (MQ=255)

                                                                                                                       aCGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATg                               >  1:4418171/1‑150 (MQ=255)

                                                                                                                         gAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTc                             >  1:3008949/1‑150 (MQ=255)

                                                                                                                           agagCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCgg                           >  1:949211/1‑150 (MQ=255)

                                                                                                                                    cGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGt                  >  1:803089/1‑150 (MQ=255)

                                                                                                                                    cGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGt                  <  2:1523113/150‑1 (MQ=255)

                                                                                                                                    cGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGt                  >  1:162878/1‑150 (MQ=255)

                                                                                                                                    cGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGt                  <  2:342100/150‑1 (MQ=255)

                                                                                                                                         tCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgc             >  2:522039/1‑150 (MQ=255)

                                                                                                                                           gacgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc           >  2:2259114/1‑150 (MQ=255)

                                                                                                                                           gacgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc           >  1:500171/1‑150 (MQ=255)

                                                                                                                                           gacgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc           >  1:1160359/1‑150 (MQ=255)

                                                                                                                                           gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgc           >  2:573153/1‑150 (MQ=255)

                                                                                                                                            acgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgcg          >  1:738084/1‑150 (MQ=255)

                                                                                                                                             cgaGCTCTGGCACGCGCTCGTCCTCTCCGGCATCGGCGGCGCGACGATCGCGCAGGCCAAAGAACCCCTGAGCCTGCGTGAGTTCAGGTCGTGGGAGAAGTACAGGCAGCTCCGGGGCTCGCTGCATGTCGGGATGCGCGTGGAgcgcgc         <  2:1241403/150‑1 (MQ=255)

                                                                                                                                                 ctctGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACg     >  2:2829270/1‑150 (MQ=255)

                                                                                                                                                    tGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGc  <  1:1078173/150‑1 (MQ=255)

CGCTTGACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGC  >  NC_002947/1765137‑1765434

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