BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,330,508	T→A	5.5%	noncoding (133/149 nt)	PP_mr52 →	PrrF ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,330,508	0	T	A	5.5%	1163.6 / 36.6	346	noncoding (133/149 nt)	PP_mr52	PrrF ncRNA
Reads supporting (aligned to +/- strand): ref base T (149/177); new base A (10/9); total (159/187)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 350 total aligned reads displayed.

AATAATCGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACC  >  NC_002947/5330359‑5330655

aataatCGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAt                                                                                                                                                     <  2:6783797/150‑1 (MQ=255)

    accGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGc                                                                                                                                                 <  2:5964134/148‑1 (MQ=255)

     tCGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGcc                                                                                                                                                <  2:6995779/150‑1 (MQ=255)

      cGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCg                                                                                                                                               <  2:3631646/150‑1 (MQ=255)

        gggCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCggg                                                                                                                                             <  1:8291905/150‑1 (MQ=255)

          gCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCGGGTc                                                                                                                                           <  1:8587968/150‑1 (MQ=255)

          gCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAAATTGCCGGGTc                                                                                                                                           <  1:2255267/150‑1 (MQ=255)

                 gATTGACGGCTAAATGAGAATTGTTATTCTGCTCGCACCTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTTTCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCttttttt                                                                                                                                    <  2:4466211/150‑1 (MQ=38)

                  aTTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTg                                                                                                                                   <  2:249369/150‑1 (MQ=255)

                   ttGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGc                                                                                                                                  <  1:2630366/150‑1 (MQ=255)

                       cGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCtgt                                                                                                                              <  1:1927293/150‑1 (MQ=255)

                       cGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCtgt                                                                                                                              <  2:2935369/150‑1 (MQ=255)

                        ggCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCtgtg                                                                                                                             >  1:1713440/1‑150 (MQ=255)

                         gCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGc                                                                                                                            <  1:6817343/150‑1 (MQ=255)

                         gCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGc                                                                                                                            <  1:5665435/150‑1 (MQ=255)

                               tGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTc                                                                                                                      <  2:8662112/150‑1 (MQ=255)

                                gagaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                     <  2:5751894/150‑1 (MQ=255)

                                gagaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                     <  2:529905/150‑1 (MQ=255)

                                gagaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                     <  1:5936779/150‑1 (MQ=255)

                                gagaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                     <  1:8818268/150‑1 (MQ=255)

                                gagaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                     <  2:3505027/150‑1 (MQ=255)

                                gagaATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCGCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCa                                                                                                                     <  1:5204432/150‑1 (MQ=255)

                                    aTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGAt                                                                                                                 >  1:1142151/1‑150 (MQ=255)

                                    aTTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGAt                                                                                                                 <  1:211045/150‑1 (MQ=255)

                                     ttGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATc                                                                                                                >  1:7098743/1‑150 (MQ=255)

                                      tGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCa                                                                                                               >  1:7391876/1‑150 (MQ=255)

                                      tGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCa                                                                                                               >  1:3559365/1‑150 (MQ=255)

                                       gTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGAACCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAg                                                                                                              >  2:7780036/1‑150 (MQ=255)

                                           ttattatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgcc                                                                                                          >  1:4395432/1‑150 (MQ=255)

                                           ttattatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgcc                                                                                                          <  2:3882219/150‑1 (MQ=255)

                                           ttattatCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGcgcc                                                                                                          >  2:7245047/1‑150 (MQ=255)

                                                 tCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCt                                                                                                    <  2:6504985/150‑1 (MQ=255)

                                                   gCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTgc                                                                                                  <  2:1494617/150‑1 (MQ=255)

                                                   gCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTgc                                                                                                  <  1:3171778/150‑1 (MQ=255)

                                                     aaCTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTgcgc                                                                                                >  1:3620542/1‑150 (MQ=255)

                                                      aCTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTACTCACGGTTATTGACCCGGCAATTTGCCGGGTCTATTTTTGCCTGTGCCTTTTCAGGATCAGTGCCGCAGCTGCGca                                                                                               <  2:5827709/150‑1 (MQ=255)

                                                       cTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGcac                                                                                              >  1:415926/1‑150 (MQ=255)

                                                            cgcgAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACaataa                                                                                         <  2:3683567/150‑1 (MQ=255)

                                                              cgAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATc                                                                                       <  2:1905901/150‑1 (MQ=255)

                                                              cgAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATc                                                                                       <  2:2158055/150‑1 (MQ=255)

                                                              cgAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATc                                                                                       <  2:8781477/150‑1 (MQ=255)

                                                               gagaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATcc                                                                                      <  1:3365043/150‑1 (MQ=255)

                                                                agaCTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGAACAATAATCCt                                                                                     <  2:3389975/150‑1 (MQ=255)

                                                                    tGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTc                                                                                 >  2:477485/1‑150 (MQ=255)

                                                                      gTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCg                                                                               >  1:549937/1‑150 (MQ=255)

                                                                         ggATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCa                                                                            <  1:6195045/150‑1 (MQ=255)

                                                                           aTAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATc                                                                          <  2:8208931/150‑1 (MQ=255)

                                                                           aTAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATc                                                                          >  1:6114352/1‑150 (MQ=255)

                                                                            tAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATcg                                                                         >  2:8341853/1‑150 (MQ=255)

                                                                               gCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATcgccg                                                                      >  2:7088054/1‑150 (MQ=255)

                                                                               gCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATcgccg                                                                      >  2:1515164/1‑150 (MQ=255)

                                                                               gCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGGTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCGGCGCACAATAATCCTGCTCCGGCATcgccg                                                                      >  1:1265382/1‑150 (MQ=255)

                                                                                cTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTCTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATcgccgc                                                                     >  1:2871523/1‑150 (MQ=255)

                                                                                     agCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgtgt                                                                <  1:8457837/150‑1 (MQ=255)

                                                                                     agCCCCTGGTTCGGACTCTCAGATGATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCgtgt                                                                <  2:1974725/150‑1 (MQ=255)

                                                                                       ccccGGGTTCGGACTGTCATCTTATCTCCTCATCAGGCTTTTGACGGTTTTGGACCCGGCTTTTTGCCGGGTCTTTTTTTGCCCGTGCCATTTCAGGATCAGCGCCGCAGCTGCGCACAATTATGGTGCTCCGGCATCGCCGCCGTGTag                                                              <  2:273301/150‑2 (MQ=255)

                                                                                        cccTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTAcc                                                             <  1:2349946/150‑1 (MQ=255)

                                                                                          cTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACcac                                                           >  2:6247835/1‑150 (MQ=255)

                                                                                           tGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACcaca                                                          <  2:4260559/150‑1 (MQ=255)

                                                                                           tGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACcaca                                                          <  1:3289264/150‑1 (MQ=255)

                                                                                           tGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACcaca                                                          <  1:3920754/150‑1 (MQ=255)

                                                                                                   tctctCAGATTATCTCCTCATCTGGCTAATCTCGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCa                                                  <  2:5004254/149‑1 (MQ=255)

                                                                                                       tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCt                                              >  2:4085378/1‑150 (MQ=255)

                                                                                                       tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCt                                              >  2:7809025/1‑150 (MQ=255)

                                                                                                       tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCt                                              >  1:4743171/1‑150 (MQ=255)

                                                                                                       tcAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCt                                              >  1:4272536/1‑150 (MQ=255)

                                                                                                        cAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTc                                             <  2:2543750/150‑1 (MQ=255)

                                                                                                           aTTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGcc                                          <  1:3811388/150‑1 (MQ=255)

                                                                                                            ttATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCa                                         >  1:665418/1‑150 (MQ=255)

                                                                                                             tATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAg                                        <  1:1578510/150‑1 (MQ=255)

                                                                                                              aTCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAgg                                       <  1:3091310/150‑1 (MQ=255)

                                                                                                                 tcctcATCAGGCTAATCACGGTTATTGACCCGGCACTTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTc                                    >  2:1929936/1‑150 (MQ=255)

                                                                                                                 tcctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTc                                    >  1:8635766/1‑150 (MQ=255)

                                                                                                                  cctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTATTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCg                                   <  1:8414724/150‑1 (MQ=255)

                                                                                                                   ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATc                                                                          <  1:5034421/110‑1 (MQ=255)

                                                                                                                   ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCgg                                  <  2:6586937/150‑1 (MQ=255)

                                                                                                                   ctcATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCgg                                  >  2:1677972/1‑150 (MQ=255)

                                                                                                                    tcatcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGt                                 >  2:7312068/1‑150 (MQ=255)

                                                                                                                      atcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtc                               >  1:7218305/1‑150 (MQ=255)

                                                                                                                      atcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtc                               >  1:7995882/1‑150 (MQ=255)

                                                                                                                       tcaGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcg                              >  2:463705/1‑150 (MQ=255)

                                                                                                                         aGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgct                            >  1:4556501/1‑150 (MQ=255)

                                                                                                                         aGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgct                            <  1:4210297/150‑1 (MQ=255)

                                                                                                                           gCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctcg                          >  1:2699706/1‑150 (MQ=255)

                                                                                                                            cTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTtcgctcgc                         >  1:627345/1‑150 (MQ=255)

                                                                                                                             tAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGcc                        <  1:7429387/150‑1 (MQ=255)

                                                                                                                             tAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGcc                        <  1:6469876/150‑1 (MQ=255)

                                                                                                                                  aCGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCt                   >  1:7773200/1‑150 (MQ=255)

                                                                                                                                    ggTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTca                 <  2:7550330/150‑1 (MQ=255)

                                                                                                                                    aaaaaaagACCCGGCAAATTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTca                 >  2:2315090/8‑150 (MQ=255)

                                                                                                                                     gTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTcag                >  1:7067185/1‑150 (MQ=255)

                                                                                                                                     gTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTcag                <  2:1324490/150‑1 (MQ=255)

                                                                                                                                       catggACCCGGCACTGTGCCGCGTCTTATCTTGCCTGTGCCATCGCAGGATCAGCGCCGCAGCGGCGCACACTAATCCGGCTCCGGCATCGCCGCCGGGTACCCCACAGACTCAGCCTCGCCCGGCGCGGTTCGCTCGCCCACCTcagcc              <  2:1708551/146‑1 (MQ=255)

                                                                                                                                          tGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTcagccagc           <  1:1497140/150‑1 (MQ=255)

                                                                                                                                           gACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCcagca          <  2:5131829/150‑1 (MQ=255)

                                                                                                                                            aCCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCcagcag         >  1:686221/1‑150 (MQ=255)

                                                                                                                                            aCCCGGCAAATTGCCGGGTGTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCAGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCcagcag         >  2:4960121/1‑150 (MQ=255)

                                                                                                                                              ccGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAgcagca       >  1:2812311/1‑150 (MQ=255)

                                                                                                                                               cGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCaa      >  2:3139832/1‑150 (MQ=255)

                                                                                                                                                 gCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATa    <  1:2508711/150‑1 (MQ=255)

AATAATCGGGGCAAAGGGATTGACGGCTAAATGAGAATTGTTATTATTATCGCAACTGGTCGCGAGACTGGTTGGATAAGCTAAGAGCCCCTGGTTCGGACTCTCAGATTATCTCCTCATCAGGCTAATCACGGTTATTGACCCGGCAATTTGCCGGGTCTTTTTTTGCCTGTGCCTTTTCAGGATCAGCGCCGCAGCTGCGCACAATAATCCTGCTCCGGCATCGCCGCCGTGTACCACACATAATCAGCCTCGCCAGGCTCGGTTCGCTCGCCCACCTCAGCCAGCAGCAATACC  >  NC_002947/5330359‑5330655

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