BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,545,558	C→G	5.4%	S9S (TCC→TCG)	rnz →	ribonuclease Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,545,558	0	C	G	5.4%	829.3 / 27.8	240	S9S (TCC→TCG)	rnz	ribonuclease Z
Reads supporting (aligned to +/- strand): ref base C (106/121); new base G (7/6); total (113/127)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 246 total aligned reads displayed.

GACACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCT  >  NC_002947/4545409‑4545706

gACACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTcc                                                                                                                                                      <  1:889773/150‑1 (MQ=255)

      aGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCggg                                                                                                                                                <  1:3151137/150‑1 (MQ=255)

       gCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCgggg                                                                                                                                               >  2:4238388/1‑150 (MQ=255)

        cAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGt                                                                                                                                              >  1:2666650/1‑150 (MQ=255)

         aaGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTg                                                                                                                                             >  2:1887675/1‑150 (MQ=255)

           gCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGcc                                                                                                                                           >  1:1539483/1‑150 (MQ=255)

            cTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGccc                                                                                                                                          <  2:681447/150‑1 (MQ=255)

            cTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGccc                                                                                                                                          <  2:4433747/150‑1 (MQ=255)

             tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCcca                                                                                                                                         >  1:3552300/1‑150 (MQ=255)

             tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                         <  1:2502896/150‑3 (MQ=255)

             tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                         <  1:2264661/150‑3 (MQ=255)

             tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                         <  1:2993907/150‑3 (MQ=255)

              cATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCccac                                                                                                                                        <  2:3094656/150‑1 (MQ=255)

               aTCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCccacc                                                                                                                                       <  1:947194/150‑1 (MQ=255)

                    gCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGc                                                                                                                                  <  1:2481673/150‑1 (MQ=255)

                     cGACAGGATAACGAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCa                                                                                                                                 <  1:437178/150‑1 (MQ=255)

                       aCAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACg                                                                                                                               >  2:4180268/1‑150 (MQ=255)

                           gATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAAt                                                                                                                           >  2:1240799/1‑150 (MQ=255)

                           gATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAAt                                                                                                                           <  2:1035768/150‑1 (MQ=255)

                               aCTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGTGTGCCCACCAAGGCACGCAATGTCa                                                                                                                       <  1:3837259/150‑1 (MQ=255)

                                cTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAg                                                                                                                      >  1:3756971/1‑150 (MQ=255)

                                 tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgc                                                                                                                     <  2:4014993/150‑1 (MQ=255)

                                   gcgggcggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgcgc                                                                                                                   <  1:3998875/150‑1 (MQ=255)

                                        cggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAg                                                                                                              >  2:4357064/1‑150 (MQ=255)

                                         ggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGc                                                                                                             >  2:1078250/1‑150 (MQ=255)

                                           tCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCg                                                                                                           >  1:2772730/1‑150 (MQ=255)

                                            cGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGt                                                                                                          <  2:3573740/150‑1 (MQ=255)

                                              aCGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGa                                                                                                        <  1:2862897/150‑1 (MQ=255)

                                              aCGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGa                                                                                                        <  2:1506921/150‑1 (MQ=255)

                                              aCGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGa                                                                                                        <  2:1729094/150‑1 (MQ=255)

                                               cggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                       >  1:1825138/1‑150 (MQ=255)

                                               cggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                       <  2:3550018/150‑1 (MQ=255)

                                               cggcggGTGGCAACCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                       <  2:3899859/150‑1 (MQ=255)

                                                ggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATc                                                                                                      <  2:1984414/150‑1 (MQ=255)

                                                ggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATc                                                                                                      <  2:3205440/150‑1 (MQ=255)

                                                ggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATc                                                                                                      >  1:3121176/1‑150 (MQ=255)

                                                   gggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                   <  1:109274/150‑1 (MQ=255)

                                                   gggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                   <  1:3498576/150‑1 (MQ=255)

                                                   gggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                   <  1:1951904/150‑1 (MQ=255)

                                                   gggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                   <  2:3050438/150‑1 (MQ=255)

                                                     gTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGc                                                                                                 >  2:1411356/1‑150 (MQ=255)

                                                        gCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAg                                                                                              <  2:1114623/150‑1 (MQ=255)

                                                          aaTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCg                                                                                            >  1:659358/1‑150 (MQ=255)

                                                             ccGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCa                                                                                         <  2:3413702/150‑1 (MQ=255)

                                                               ggCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGc                                                                                       <  1:4108417/150‑1 (MQ=255)

                                                               ggCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGc                                                                                       <  2:2139782/150‑1 (MQ=255)

                                                                    aGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTg                                                                                  <  2:2440027/150‑1 (MQ=255)

                                                                      gtgtGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGt                                                                                >  2:3045958/1‑150 (MQ=255)

                                                                      gtgtGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGt                                                                                >  1:1401846/1‑150 (MQ=255)

                                                                       tgtgCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTa                                                                               >  2:1403417/1‑150 (MQ=255)

                                                                         tgCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCCAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTg                                                                                  <  2:1628743/145‑1 (MQ=255)

                                                                         tgCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATc                                                                             <  1:4439005/150‑1 (MQ=255)

                                                                             ccccGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGt                                                                         <  1:1881473/150‑1 (MQ=255)

                                                                              cccGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTc                                                                        <  1:197438/150‑1 (MQ=255)

                                                                              cccGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTc                                                                        >  2:225641/1‑150 (MQ=255)

                                                                                   aTCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTg                                                                   <  1:2585760/150‑1 (MQ=255)

                                                                                    tCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgc                                                                  >  2:1581272/1‑150 (MQ=255)

                                                                                     cACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcg                                                                 <  2:3074744/150‑1 (MQ=255)

                                                                                     cACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcg                                                                 >  1:761815/1‑150 (MQ=255)

                                                                                      aCCTTGGCCGCTTCTGGCCTGAAACCGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcgg                                                                <  2:591785/150‑1 (MQ=255)

                                                                                      aCCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcgc                                                                <  2:1727923/150‑2 (MQ=255)

                                                                                       ccTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcggc                                                               >  1:2623979/1‑150 (MQ=255)

                                                                                          tGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                            <  1:4146036/150‑1 (MQ=255)

                                                                                          tGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                            <  1:1644410/150‑1 (MQ=255)

                                                                                          tGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                            <  1:3985436/150‑1 (MQ=255)

                                                                                           ggCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTCTCTGGTCGATTGCGGCGAAg                                                           >  2:932643/1‑150 (MQ=255)

                                                                                            gCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAgg                                                          <  2:4500697/150‑1 (MQ=255)

                                                                                               gCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAc                                                       >  1:1373263/1‑150 (MQ=255)

                                                                                                  tCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCa                                                    <  1:1936297/150‑1 (MQ=255)

                                                                                                   cTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCACTGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAccccg                                                   >  2:3347910/1‑148 (MQ=255)

                                                                                                     ggCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCa                                                 <  2:1519246/150‑1 (MQ=255)

                                                                                                      gCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCAc                                                >  2:3463114/1‑150 (MQ=255)

                                                                                                        cTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCg                                              >  1:847454/1‑150 (MQ=255)

                                                                                                                      tGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACAccc                                >  2:137729/1‑150 (MQ=255)

                                                                                                                      tGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACAccc                                >  1:2888801/1‑150 (MQ=255)

                                                                                                                             ggACCGGCTGTTCCGTGGCACCTGCGCCGGGGTGCGCACCAAGGCACGCAATGTCTGCGCGACAGCCGTGATCGAAGCCAGCGGCAGCCTCTGGTATCTGGTCGATTGCGGCGTAGGCACCCTGCTCCGCCTGCTGCACACCCCGCTTTc                         <  2:4376334/150‑1 (MQ=255)

                                                                                                                                   gctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcg                   >  1:1672051/1‑150 (MQ=255)

                                                                                                                                   gctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcg                   <  1:281076/150‑1 (MQ=255)

                                                                                                                                    ctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgc                  <  1:646304/150‑1 (MQ=255)

                                                                                                                                     tgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgcg                 >  1:1008405/1‑150 (MQ=255)

                                                                                                                                     tgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTGTCGATCcgcg                 >  2:1617579/1‑150 (MQ=255)

                                                                                                                                       ttggtgggCACCCCGGCGGAGGTG‑CCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGAcc              >  2:871189/7‑150 (MQ=255)

                                                                                                                                       ttCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGAc               >  2:3988418/1‑150 (MQ=255)

                                                                                                                                         ccTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCCGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCt             >  2:2320081/1‑150 (MQ=255)

                                                                                                                                         ccTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCt             >  1:3586359/1‑150 (MQ=255)

                                                                                                                                          cTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTg            >  1:3338551/1‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCATCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc           >  1:2356273/1‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc           >  2:3911111/1‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc           >  1:3991876/1‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc           <  2:4311158/150‑1 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc           >  2:1005672/1‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc           >  1:1257821/1‑150 (MQ=255)

                                                                                                                                           tgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc           >  2:2944956/3‑150 (MQ=255)

                                                                                                                                              gCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCCGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTATCGATCCGCGACCTgcgcg        >  2:2566111/1‑150 (MQ=255)

                                                                                                                                              gCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgcg        <  2:2756431/150‑1 (MQ=255)

                                                                                                                                              gCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgcg        >  1:3726282/1‑150 (MQ=255)

                                                                                                                                               cACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgcgc       <  1:3838631/150‑1 (MQ=255)

                                                                                                                                                aCCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGcgcgcg      >  1:159500/1‑150 (MQ=255)

                                                                                                                                                  cTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGAt    <  2:1367711/150‑1 (MQ=255)

                                                                                                                                                   tCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATc   >  1:2558140/1‑150 (MQ=255)

GACACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCT  >  NC_002947/4545409‑4545706

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